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ABSTRACT 



A procedure to approximate data given at arbitrary intervals in 
two-, three- and four dimensions using polynominal expressions is 
described. The programming of the problem is explained in each case 
and a user manual is given for software implemented on the TPL 
Hewlett-Packard computer system. The method, which is general, xvas 
derived and has been applied to represent the calibration of flow 
probes which have multiple sensors o 
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1 . INTRODUCTION 



This report describes algorithms to approximate data patterns such 
as those which occur, for instance, when pneumatic-velocity probes are 
calibrated. The functional value (Y) can depend on either one (XI), 
two (XI and X2) or even three (XI, X2 and X3) parameters. The approxi- 
mation Y = f(Xl) will be referred to as the two dimensional, the approxi- 
mation Y = f(Xl, X2) as the three dimensional and the approximation 
Y = f(Xl, X2, X3) as the four dimensional approximation. These three 
options meet the requirements for probes used in the Turbopropulsion 
Laboratory. 

Using the least-squares criterion to obtain the coefficients in 
assumed polynominals leads to a system of linear equations, which are, in 
principle easily solved. Numerical problems may arise however, since 
the software, as described herein and as implemented in the NFS 
Turbopropulsion Laboratory Hewlett Packard 21MX computer, uses 32-bit 
real constants. 

The author wishes to express his thanks to Professor Ray Shreeve, 
whose constant and critical interest was a crucial help to solve this 
mathematical problem and link it to engineering application. 
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2. TOO DIMENSIONAL APPROXIMATION 



2.1. Problem: 




' 1 > 

Xi Xi 



r 

A data set of NPNTSl data points is given, where Y depends on 
parameter XI. The data pattern is to be approximated by a function 
Y = f(Xl), so that the error between data points and analytically 
determined points is lowest. 
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2.2, Approach : 



Polynorainals are commonly used in order to approximate data 
patterns. Y = f (XI) is a function which approximates the data value, 
V, at each value of the variable XI, we look for our expression of 
the form 

2 ( L. i ) 

Y = C + C *Xi + C JtcXl** + ... + C *Xi 
12 3 L 



or 



Y 



L 

= S c *xi 

i = l 1 



<i-l > 



( 2 . 1 .) 



in which the coefficients C. are to be determined by the method of 

1 

least squares. As Ref. 1 shows Che least squares criterion leads to 
a linear equation system that can easily be solved. We define the error 



NPNTSi 2 

R = E [f(Xi )-Y ] 

r=l r r 



( 2 . 2 .) 



where the index r denotes the individual data point. Using equation 
(2.1.) R becomes 



NPNTSi 2 <L-i) 

S tC + C 4cxi + C *Xi + . . . + C «Xi - Y 

= i i 2 r 2 r L > 
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2.3. Solution: 



R depends on the selection of the coefficients 

In order to minimize the error, R is differentiated with respect 

to C , ...,C and the partial derivatives are set to zero. Thus 
1 u 



3 R ! 

3 C ~ ° 
i 



i = i, . . . ,L 



or 

3R NPNTSi 2 (L-i) 

= E 2*LC +C #Xi +C *Xi + . . . + C Xi -Y 

5C r = l 12r3r Ir r 

i 

^ 2 (L-i) 

^ [C +C *Xi +C «Xi + . . . + C XI -Y 1 

12r3r Lr r 

t ^ = 0 



^ c 



I 



Assuming that the summations extend over all the data points, Z 

NPNTSi 

should be understood to mean E • Performing the differen- 

r = i 

tiations and rearranging the equations, we get 



NPNTSi*C 

t 


+ Zxi 

r 2 


2 

+ Zxi »C + , 
r 3 


V. (L-i) 
. . + EXi 


*c = Cy 

L r 


Z^Xi :*:C 


+ Z7xi‘'*c 


+ Z!xi^*c + . 


. . + 2Hxi^ 


*c = Ey *Xi 


r i 


r 2 


r 3 


r 


1 r 


r-. 2 

CXi *C 
r i 


+ ZIXi^»C 
r 2 


+ Z’xi'^*c + . 

r 3 


m <L+1) 

. . +Z7Xi 

r 


*c =Cy *xi 

L r 



(L-i) ml- m (L + i) ^ (2L-2) I -i 

27xi *c +I7xi +Exi *c +.. + i:xi »c =i;y »xi 



r r 
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In matrix notation 



( 




Q 

c. 


L-i > 




/ \ 


( 


NPNTSl 


Exi 


Exi 


. . . Exi 




c 


\ 




r 


r 


r 




1 






2 

EXi 


3 


L 








CXi 


EXi 


. . . EXi 




c 


\ 


r 


r 


r 


r 




2 




2 


3 


4 


L + i 


* 




r= 


EXi 


EXi 


EXi 


. . . EXi 




C 




r 


r 


r 


r 




3 




L-1 


L. 


L+1 


2L-2 








EXl 


EXi 


EXi 

r* 


. . . EXi J 




C 

k I 1 


1 



Cy 



r 

o 



EY *Xi 






% C 



= !■:< (2 



... Systen Matrix 
C' ... Coefficients vector 

... Right hand side vector 
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2.4, 



Equation System 



2,4,1. System Matrix A 

Since all Matrix elements on diagonal lines running up from 
left to right (2. order, as defined in Appendix A) are identical, the 
elements are renamed for simplification as follows: 



















... 








LA ^ 




= ^ ^Z,U 



so that the set of elements, 






can be written as 



HPf^rz! 

/ .. 2L-/ 

tz.\ ^ f f 



(2.4.) 



This operation is programmed in REAL FUNCTION S2 (NPNTSl , IPOWRl , lY) . 
The data XI and Y are known to this function through a common block 
named DTA2. So only NPNTSl, IPOWRl (= k-1) and lY (=0) have to be 
passed to the function, and the value of a^ is returned through the 
function name S2. 
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SUBROUTINE MAT2 presets the system matrix in the following way: 



i) Preset edge section elements (using REAL FUNCTION S2) 



I 

^L*Z 

a, . 



lop 

^ 2l-I 

Aa^e/ s<^h6^ 



ii) Copy defined elements diagonally 




O, 



*!Z 



a, 







Oil 



a 



SL 



^LL 
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a. 



a. 



'Zl 



<7,7 



% ‘^23 



a 



IL 



■‘ZL 



y y y 









Gl_l 


O'/, 




0,Z 


• • • 


Oz, 


On 




♦ - . ^Zl 


^3l 


^3Z 


Ojs 


03 L 



m • • ••• «•# >v*** 

2 ::/^ ^ 







^>2 




• • • 


0 ,, 


j 






Ou ^ 2 Z 


O 23 








A 


— 


O 31 


<=^Z3 


• • • 










* • • • • • 






• • • 








^LZ 






OCL 




2.4.2. 


Right hand side vector B 










The elements of the 


right 


hand side 


vector B can be written as 




NTWI3I 

UX-K) 

\ 


k = 


h ••• 




(2.5.) 


To calculate 


b, , again REAL 


FUNCTION S2 is 


used . 


Data XI and Y are 



available through COMMON block DTA2. NPNTSl, IPOWRl (=k-l) and lY (=1) 



have to be passed to the function, and the value of b is returned 

IX 

through the function name S2. 
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2.5. Software: 



The software to compute the coefficients for a two dimensional 
approximation is described in Appendix B and is implemented in the 
TPL HP-21MX computer system 

To work correctly with these program modules, the user must 



conform to the following conventions: 



i) Provide the data in two arrays (Type; REAL) of 256 
elements through a COMMON block, named DTA2. 

COMMON / DTA2 / X1,Y 
REAL X1(256),Y(256) 

ii) Dimension an array (Type: REAL) of 7 (seven) elements 
to contain the coefficients. 



REAL C0EE(7) 



iii) Define the parameters NPNTSl, L and IPRINT (Type: all 
INTEGER) 

NPNTSl ... # of data points 

L ^ NPNTSl ^ 256 

L ... (desired order of polynominal) + 1 

1 ^ L ^ 7 



IPRINT 



controls 

2 



1 

<0 . . . 

>0 . . . 



quantity of print out 
Print system matrix and right 
hand side vector before and 
after Gauss Jordan Elimination 
Print equation system after 
Gauss Jordan Elimination 
No print out 

Print equation (1,1) with the 
actual parameters 



iv) When loading a program, that uses the subroutine MAT2, 

the binary library file has to be searched for externals. 



Suppose, the source file &USER::26 contains a user 
program named USER. This program calls MAT2. After 
the compilation, the relocatable binary file for this 
user program is %USER::26. To load the program, the 
following procedure is recommended: 
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Type (from FMGR) 

:RU,LOADR 

and the loader program will respond 
/LOADR: RE,%USER; :26 

where the underlined information already is the user's 
input. This causes the loader to load all program 
modules of the user program; a load map is listed on 
the terminal. Upon completion the loader prompts 

/LOADR: MS,%TPLBL: :26 

where the underlined information already is the user's 
input. Now the loader conducts a search for all unsatis- 
fied externals of the user program in the binary lib- 
rary file %TPLBL::26. Since some library programs 
have externals themselves the search has to be repeated 
(MS . . . multiple search) as many time as is necessary 
for all externals to be satisfied. The loader prompts 

/LOADR: END 



where the underlined information is the user's input. 
The loader now loads all system programs, outputs a 
load map and generates the program. Upon completion, 
the loader outputs a ready message 

/LOADR:USER READY AT 1:28 PM WED., 10 SEPT, 1980 
/LOADR: END 

Now the program can be run 
: RU,USER 

Again, the underlined information is the user's input. 

If all these requirements are met, the correct call for the sub- 
routine is: 



CALL MAT2 (NPNTSl ,L, COEF , IPRINT) 

Upon completion, the array COEF contains the coefficients (C0EF(1)=C^ 
, . . . ,C0EF (L)=C^) . Externals used by MAT2 are: AB2, DTA2, S2 under 

no circumstances may the user use any of these names for modules of 
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his own user program. In some cases the program may not be able to 
perform a Gauss Jordan Elimination to the system matrix and the right 
hand side vector and thus cause the program to stop. If this happens, 
an error message is displayed. 

It is highly encouraged, to use the system function FNP to 
calculate the value of a polynominal at one specified XI rather than a 
loop such as the following: 

S=.0 

DO 08 1=1, L,1 
08 S=S+X’*f*(I-l)*COEF(I) 

YCALC=S 

To call FNP is easier and faster, because FNP is a programmed Horner 
Scheme. The recommended call is then; 

L1=L"1 

YCALC=FNP ( COEF , X , LI ) 
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2.6. Sample User Program 
2.6.1. Listing 

0001 FTN. 12:40 PM THU., 18 SEP., 1980 



FTN4,L 

C 
C 
C 
C 
C 
C 
C 
C 



PROGRAM DEM02 (3,99) 



THIS IS A DEMONSTATION PROGRAM AND IT SHOWS THE CORRECT USE 
OF THE TPL BINARY LIBRARY (ON TYPE 6 FILE %TPLBL) USING 
EXAMPLES. FOR FURTHER QUESTIONS READ THE TPL-LIBRARY-BINDER 
OR CONSULT THE SOFTWARE MANAGER . 



*, DEMONSTRATE THE USE OF THE TPL BINARY LIBRARY. 

COMMON / AFLD / A 




INTEGER N0LF,N0CR(2) ,ICLR(3) 



DATA NOLF /006537B/ 

DATA NOCR /OO 0033B, 040433B/ 

DATA ICLR /015524B, 015515B,006537B/ 

DATA PI /3. 141593/ 

r-" a R M *T' i> ih: m a r a r-i t 

101 FORMAT (///" HELLO! THIS IS PROGRAM DEMO AND WE»LL 
»OW TO USE THE Mi^RUELLOUS TPL BINARY"A2/“ LIBRARY! 
*WILL BE GREATLY AWARDED BY EASIER PROGRAMMING."///) 

102 FORMAT 



SEE H 

YOUR INTERES1 



I) ") 



FORMAT (" GENERATING DATA POINTS "/5X" I "5X"X( I ) "5X " Y ( 

103 FORMAT ( 2X , 14 , 2F9 . 3 , A2 ) 

104 FORMAT (" DATA POINTS GENERATED!"/) 

105 FORMAT (" INITIALIZING THE PLOTTER") 

106 FORMAT (" PLOTTER INITIALIZED!"/) 

107 FORMAT (" DEFINE PLOTTER AND USER. AREAS AND DRAW AXES") 

108 FORMAT (9X""6X"XMIN "6X“XMAX "6X"YMIN 
*"6X"YMAX"/9X,4( ""lOX" ")/ 

*" HP9872".4(" "FIO .3" " ) . ""/ 

* ) 

P LO T T t'R Al^’ffS^EF I NED AND AXES DRAWN!"/) 

DRAWING DATA POINTS INTO COORDINATES SYSTEM") 
DATA POINTS DRAWN!"/) 

CALCULATING A CURUE FIT THROUGH THE DATA POINTS 



*" 05ER 

109 FORH75rr~ 

110 FORMAT 

111 FORMAT 

112 FORMAT 
«") 

113 FORMAT 
*UGH THE 

114 FORMAT 



( " ENTER DEGREE 
POINTS "2A2) 



OF POLYNOMIAL (NORDER) TO FIT THRO 



(" ENTER IPRINT "2A2) 

115 FORMAT (/" CURVE FIT DONE!"/) 

116 FORMAT (" CURVE FIT DONE!"/) 

117 FORMAT (" PLOTTING THE CURVE FIT") 

118 FORMAT (" CURVE FIT DRAWN!"/) 

149 FORMAT (""((3A2) ) ) 

I" O R M T* F:5 I> In. M O 5!> T' Q \-> 



C 

C 

c 

c 

c 



GET THE LU OF THE TERMINAL. 



LI = LOGLU(I) 

WRITE (LI, 101) NOLF 



C 

C 

C 

C 

C 



GENERATE DATA POINTS SCATTERED AROUND A POLYNOMIAL. 



CALL INITG (13) 
WRITE (LI, 102) 
XSTART = -6. 
XSTOP = +2. 
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0002 DEM02 12:40 PM THU., 18 SEP., 



1980 



DX = (XST0P-XSTART)/<NPNTS1-1) 

NORDER = 4 
COEFd) = -1 . 

C0EF<2) = -0.500000 
C0EF(3) = -0.000340 
CDEF<4) = +0.300000 
C0EF<5) = +0.055000 
A = <2. 0*PI)*0 . 03125 

DO 01 I=1,NPNTS1,1 
XKI) = XSTART+(I-1)*DX 

Y<I) = FNP<CQEF,X1< I) , NORDER) + 0 . 25*SIN< XI < I ) /A ) 

01 WRITE <LI, 103) I,XKI),Y(I),NOLF 
ISTOPl = NPNTSl-1 

DO 02 1=1, ISTOPl, 2 
DUMMY = Y<I+1) 

Y(I+1) = Y(I) 

02 Yd) = DUMMY 

WRITE <LI, 149) <ICLR,I=1,3,1) 

WRITE <LI, 104) 



C 

C 

C 

C 

C 



£ 

C 

C 

c 



INITIALIZE THE PLOTTER 



WRITE <LI, 105) 
LP = 13 

CALL INITG <13) 
WRITE (LI, 149) 
WRITE <LI, 106) 



ICLR 



DEFINE PLOTTER AND USER AREAS) DRAW AXES. 



WRITE 

XPMIN 

XPMAX 

YPMIN 

YPMAX 

XUMIN 

XUMAX 

YUMIN 

YUMAX 

WRITE 

ALPHAX 

XA 

XB 

XL 

ALPHAY 
YA 
YB 
YL 

CALL 
CALL 



(LI, 107) 

= 2 . 

= 12 . 

= 12 . 

= 22 . 

= - 6 . 

= + 2 . 

= -5. 

10 . 

, 108) 
0.00 
(XPMAX 



(LI 



CALL AXIS 
WRITE (LI, 
WRITE (LI, 



XPMIN, XPMAX, YPMIN, YPMAX, XUMIN, XUMAX, YUMIN, YUM 

-XPMIN )/(XUMAX-XUMIN) 

( XPMIN*XUMAX-XPMAX«XUMIN ) / ( XUMAX-XUMIN ) 

(XPMAX -XPMIN ) 

90 . 00 

= (YPMAX -YPMIN )/(YUMAX-YUMIN) 

= (YPMIN*YUMAX-YPMAX*YUMIN)/(YUMAX-YUMIN) 

= (YPMAX -YPMIN ) 

SETSM (113,1. ) 

AXIS ( XP MIN , YPMIN , XL , ALPHAX , 2HX 1 , +2 , XUMIN , XUMAX , 4HF4 . 1 , 4 



( XPMIN, YPMIN, YL,ALPHAY,2H Y , -2 , YUMIN , YUMAX ; 2HI2,2 , 3 
149) (ICLR, 1=1, 6,1) 

109) 



C 

C 

C 

C 

C 



PLOT DATA POINTS. 



WRITE (LI, 110) 

CALL SETSM (113,3. ) 
DO 03 I=l,NPNTSi,l 
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0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 
‘'173 

174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 
199 

^200 

0201 

0202 

0203 

0204 

0205 



C 

c 

c 

c 

c 



c 

c 

c 

c 

c 



XPLOT = X1(I)*XA+XB 
YPLOT = Y(I)»YA+YB 
CALL PLOT (XPLOT, YPLOT, 2) 
03 CALL SYMBL (1) 

WRITE (LI, 149) ICLR 
WRITE (LI, ill) 



CALCULATE CURVE FIT THROUGH DATA POINTS. 



WRITE (LI, 112) 

04 WRITE (Ll' 113) NOCR 

READ (LI, NORDER 

WRITE (LI, 149) ICLR 

IF ( NORDER .LT. 0 
WRITE (LI, 114) NOCR 



OR 



NORDER .GT. 6 ) GO TO 04 



WRITE (LI, 149) (ICLR, 1=1, 2,1) 



IFn?RINT^GE?0?*SR?TE*Tf??*i?5^ 

IF (IPRINT .LT . 0) WRITE (LI, 116) 



PLOT POLYNOMIAL FIT. 



09 



WRITE (LI, 117) 

CALL SETSM (113.2. ) 

DO 09 I=l,NPNTSi,i 
XPLOT = Xi(I)»XA+XB 

YPLOT = FNP(COEF,Xi(I) .NORDER)*YA+YB 
IF (I.EQ.l) CALL ^LOT ()i:PLOT , YPLOT , 1 ) 
IF (I.GT.i) CALL PLOT (XPLOT , YPLOT ,3 ) 
CALL SETSM (115, DUMMY) 

WRITE (LI, 149) ICLR 
WRITE (LI, 118) 



TERMINATE GRAPHICS. 



CALL STOPG 
STOP 7777 
END 



FTN4 COMPILER; HP92060-16092 REV. 1926 (790430) 

** NO WARNINGS tt NO ERRORS ** PROGRAM = 01409 COMMON = 00000 
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2.6.2. Load map 



0EMO2 

FNP 



10042 

12643 



12642 

12746 



DEMONSTRATE THE 



USE OF THE TPL BINARY LIBRARY 
Calculate n-th order polyn 



to load these programs, type (from LOADR) ; MS,%TPLBL 



iwra 



12747 

14747 

17000 

17132 



14746 

16777 

17131 

17311 



2D-Appr 0 X inat i on / DTA2 
2D-Appr 0 X inat ion arrange 
2D-Appr 0 X inat i on 
2D-Approxiwation / AB2 



/ (DATA2 DTA2) 

system matrix and vector 
compute summationS| 



AFLD 

AXIS 

INITG 

PLOT 

SETSM 

STOPG 

SYMBL 

FCTR 

FACTR 

LOGLU 

READF 

OPEN 

CLOSE 

CLRIO 

OVRD. 

$SMOE 

LURQ 

.DADS 

.DMP 

.DDI 

SESSN 

R/W* 

P .PAS 
.DNG 
PAUSE 
$ALRN 
FMTIO 
ERRO 
TAN 
ABS 
.SNCS 
.DDE 
.DIN 
.RTOI 
.FPUR 
.SBT 
. .FCM 
PAU .E 
ERO.E 
.CMRS 
$OPEN 
RU$UB 
RUNDS 
FRMTR 
FMT. E 
IB4A2 
REID 
RMPAR 
PNAME 
LUTRU 
IPUT 
SSETP 
.CFER 
.LBT 



17312 

20312 

22632 

2364S 

24376 

25337 

25402 

26137 

26143 

26163 

26241 

27203 

27501 

27711 

27720 

27721 
30010 
30373 
30503 
30651 
31152 
31170 
31327 
31356 
31366 
31467 
31605 
33104 
33174 
33301 
33310 
33452 
33464 
33472 
33566 
33630 
33671 

33706 

33707 
33710 
33774 
34151 
34517 
34642 

40300 

40301 
41602 
41727 
41772 
42040 
42147 
42170 
42215 
42273 



20311 

22631 

23644 

24375 

25336 

25401 

26136 

26142 

26162 

26240 
27202 
2750 0 
27710 
27717 
27720 
30007 
30372 
30502 
30650 
31151 
31167 
31326 
31355 
31365 
31466 
31604 
33103 
33173 
3330 0 
33307 
33451 
33463 
33471 
33565 
33627 
33670 

33705 

33706 

33707 
33773 
34150 
34516 
34641 
40277 
40300 
41601 
41726 
41771 
42037 
42146 
42167 
42214 
42272 
42323 



GSP 

GSP 

GSP 

GSP 

GSP 

GSP 

GSP 

GSP 

GSP 



/ AFLD / A(256) 
read control 



array 

move 



/ FCTR / 



print 

FX,FY 



Drew and label axes 
Aj initialize plotter 
pen to a defined point 
change plotter modes 
terminate grafics 
symbol at location of pen 

vary size of the plot 



92067-16268 REV. 1903 
92067-16125 REV. 1940 
92067-16125 REV. 1903 
92067-16125 REV. 1903 
750701 24998-16001 

92067-16125 REV. 1903 
92067-16268 REV. 1903 
92067-16268 REV. 1903 
780818 24998-16001 

780818 24998-16001 

781021 24998-16001 

92067-16125 REV. 1903 
92067-16125 REV. 1903 
92067-16125 REV. 1903 
780818 24998-16001 

771122 24998-16001 

92067-16268 REV. 1903 
24998-16002 REV. 192 
771122 24998-16001 

780424 24998-16001 

750701 24998-1600 

24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
REV. 1903 
REV. 1903 



790228 

790719 

790215 

781229 

780526 

790202 

790223 



780413 

781214 

740801 



770715 
6 790417 



780424 
780818 
780818 
780921 
781106 
770518 
750701 
750701 
750701 
780424 
92067-16125 
92067-16125 



92067-16125 REV. 1903 
24998-16002 REV. 192 
24998-16002 REV. 190 
59310-1X013 REV 1940 
92067-16268 REV. 1903 
781106 24990-16001 

771121 24998-16001 

92067-16268 REV. 1903 
92067-16125 REV. 1903 
781106 24998-16001 

750701 24998-16001 

770518 24998-16001 



790103 
781003 
780801 
6 790503 
1 781107 
790802 1153 
790316 



790223 

740801 



15 PAGES RELOCATED 15 PAGES REQ>D NO PAGES EMA 
LINKS:BP PROGRAMrLB LOAD : TE COMMON:NC 
/LOADR :DEM02 READY AT 12:44 PM THU., 18 SEPT, 1980 

/LOADR : SEND 



NO PAGES MSE 
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2 . 6.3 



Results 



Affect of IPRINT on quantity of print out 



i) printed output 

Systert Matrix A and vector B before Gauss Jordan EliMation "\ 



1 

2 - 

3 

4 - 

5 



SiOE+02 

102E+03 

487E+03 

2iiE+04 

104E+05 

i 



102E+03 
487E+03 
2i iE+04 
i04E+05 
S2SE+05 
2 



. 487E+03 
.211E+04 
. 104E+05 
. 525E+05 
. 274E+06 
3 



,2iiE+04 

104E+05 

525E+05 

274E+06 

146E+07 

4 



5 

. 104E+05 
. 52SE+05 
. 274E+06 
.146E+07 
. 787E+07 
5 



S78E+02 
.998E+02 
il3E+03 
-.944E+0 
. ii3E+0 



Equation syst 


eM after 1 


Sauss Jordan 


EliMinat 


i on 






i 


2 


3 


4 


5 




i 


. iOOE+Oi 


. OOOE-i-00 


. OOOE+00 


. OOOE+00 


. OOOE+00 


970E+00 


2 


. OOOE+00 


, iOOE+Oi 


. OOOE+00 


. OOOE-t-00 


. OOOE-i-00 


-.416E+00 


3 


. OOOE+00 


. OOOE+00 


. iOOE+Oi 


. OOOE+00 


. OOOE-i-00 


i62E-0i 


4 


. 000E-«-00 


. OOOE-t-00 


. OOOE+00 


. iOOE+Oi 


. OOOE+00 


.275E+00 


5 


. OOOEt-00 


. OOOE+00 


. OOOE-t-00 


. OOOE-t-00 


. iOOE+Oi 


. 5i4E-0i 




i 


2 


3 


4 


5 





> 



Coefficients COEF(I) 



I 

I 

I 

I 

I 



t 

2 

3 

4 

5 



970E+00 
. 4i6E+00 
, 162E-01 
. 275E+00 
,Si4E-0i 



J 



Equation syst 


eM after Gauss Jordan 


EliMinat 


ion 






i 


2 


3 


4 


S 




i 


.iOOE+Oi 


. OOOE-t-00 


, OOOE+00 


. OOOE-i-00 


. OOOE+00 


970E+00 


2 


. OOOE-i-00 


. iOOE+01 


. OOOE+00 


. OOOE+00 


. OOOE+00 


- . 416E+00 


3 


. OOOE+00 


. OOOE+00 


. iOOE+Oi 


. OOOE+00 


. OOOE-4-00 


i62E-0i 


4 


. OOOE+00 


. OOOE-i-00 


. OOOEt-00 


. lOOE+Oi 


. 000E-«-00 


.275E+00 


5 


. OOOE+OO 


. OOOE+00 


, OOOE+OO 


. OOOE+00 


. iOOE+Oi 


. 514E-01 




i 


2 


3 


4 


5 





Coefficients COEF ( I ) 



I 

I 

I 

I 

I 



1 

2 

3 

4 

5 



-.970E+00 
- .416E+00 
162E-0i 
.27SE+00 
.5i4E-0i 



J 



Coefficients COEF(I) 



I 

I 

I 

I 

I 



= i - 

= O - 

= 3 - 
= 4 
= S 



, 970E+00 
416E+00 
, i62E-0i 
27SE+00 
,Si4E-0i 



IPMNT « I 



IPRINT - 1 



IPRINT < I 
>2 
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ii) graphic output 




Grafic output from PROGRAM DEM02 
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3. THREE DIMENSIONAL APPROXIMATION 



3.1, Problem: 




A data set of NPNTS1*NPNTS2 data points is given, where Y 
depends on the parameters XI and X2, The data pattern is to be 
approximated by a function Y= f (XI, X2), so that the error between 
data points and analytically determined points is lowest. 
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3.2. Approach : 



As for the two dimensional approximation, a polynominal, now 
with two independent variables XI and X2, is used. 

Y- fCxl. XZ) 

-Cc„->c^xz 4 ... 4 C2^xz'^"’]-xi i- 

4 CCs, 4 4 Cgj X2^+ ... 4 Cj„x2 -X/ ^ 4 



4 [€^,4 C^XZ4 C^XZ^4 ... 4 c^^xz’^'^l-xi"' 



L 



M 



Y-E{ Z C,yXZ^'''’/-XJ 



(/-/) 



c= / j\. / y 



(3.1.) 



The data arrangement for the three dimensional approximation 
already yields to the application of this method: at each constant 

Xl (i.e.: r = constant, too) there is a data set for various X2. 

This is the case for example, when a pneumatic-velocity probe is 
calibrated. At each Mach number, the probe is balanced in yaw and the 
pitch angle is varied over the range of application. The probe out- 
put (each of two separate non-dimensional pressure differences) is a 
function of the Mach number and the pitch angle. 
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The least squares criterion becomes 



Z f ZC{ryi„,x 2 ^)- Y„J J 0 . 2 . 

/ £•- / 



R 



where the indices r and s denote the individual data points. Using 
equation (2.1.) R becomes 



hJPHTSf hiPHTSZ 



R 



/ 



J7 • f Z1 ^ ^ 






H-t 



- 7 V 



The term in the [ ] - bracket shall be called B. 
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3.3. Solution: 



To minimize the error, R is partially differentiated to the 

coefficients C.. and then the partial derivatives are set to zero 
ij 









Assuming, that the summations extend over all data points, ZZ 

NP MTS / 

should be understood to mean ZJ-[ u • Performing the differen- 

/ C «/ 

tiations and rearranging the equa^tions in matrix notation (which, 
while lengthy, is logically straight forward), we get, in matrix 
notation. 



3 ^ ! 

- G 



3C. 



/■- /...,4 • 
/ ) / 






^ . r-f 7 



"J 



!/ 
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/ 

-A. 



2 









/< 



/4/f 

# 

0 

t 

Z*M 



f(L-!)*-M 



L*M 



■ Ci-j}»H-il 




/*M 


(Z-i)fiM'tl 






NPHTS! f- HfNTSt 


ZZxz,^ 


SExC 


22 >/ 

FT 


HZxhXZ^ 


-•■- 


ZZxz 

. r ' 3 - 


ZExzl 


IZxzl, 


ZZxi.xz 


YZxhXZ^^ 


... 


I^EXZ'"' 

rs 


ZZxz: -- 


* # • 


* * § 

ZZxi^xz 


JUxhXzl 




LLxi^ 

ZZxi-xz^ 

* • « 1 


ZZxi^xi^ ... . 

ZZx^l ... 

\ 


22 y/X" 

ZZxi-xz^, 


- 


- T.Zxil:xz 


• • • 

I 


1 

ZZxi-xz''' ■ 

rs rs 




ZZx,^-xC 


P # « 

ZZxtxz'' 


EC^.y/ 


... 









2EyC 


EEy/">2 ... 


EC^'-yz'’" 


EE,;^ 


EEy/>2 


ZZxi^''xz'^ ... 


EE^Ui" 






Er;./>2" ... 


p o L-/ 2rf*2 

LliA'/,^*r2 


0 0 0 

Y,'Zxi}xz„ 

ts rs 






tt 



tz ^ 



HUxt-yz 

rs 



fx-f 
rs nr 



E1LX/-X2. 



2ff-Z 



«u 



ZZx>l-xz 
E r ^'\sxL 

rs rr 



n 

ft 



Y.'Zxi-xi 

(s rr 



:lM^ 



zt:x/J-x2^ 

ZZxt-xz, 



ZZxCjxz^ 



A 



L 




• • • 




0 0 W ' 




I 



0 0 0 



/ 

tixt: 

LZxC-xi^ 


eex/^Vaz^ .:: 
2ZxCW^ ... 




^ — \ 

ESxC.yz: • 

E 2 <.< 

1 




Cu 




\ 

EEV„ 

EE X,-xz^ 


— 


i^xr.xzi 


ZZx)^-xC ... 




2 z<v"; 




^Ifl 




2 E 





22 a/;, 

HZxt'Xi 

r 5 rj 

9 # * 


S 2 «>« 

... 


EEvz;; 

2 E)^i*y 2 j 




• 




22 Y„ .XI, 

EE Y„.y 2 ,.xi^ 


— 


E 2 « 


sExc^:. - .-.. 




_ _ 1 _ 2M1 

EEx 4 ^.a 2 ^ 


--- 


# 




4 

iz 


— 






i 


* 




- 






” 


i 




- 




0 

♦ 

# 


$ 


f 

• # 


— 


EEx^' 

2 Z<-a 3 . 

» • f 


ZLxi^'^-xz,, ... 

... 


1 

! 

1 

i 


21-2 fr-i 

ELX/j, 

s 2 x/r-A 2 : 




f 

0 

• 




2 EY, . <' 

E 2 Xx*X 2 „ yi, 

« 

_ 0 




EEx/:^yz:: 


EEC^: ... 


EE^'Wr 




^LM 




EEX-xz^-x.: 





* 




€ 



(3.3.) 



I 







I i 

J - 



A/ 
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3.4, Equation System : 
3,4.1. System Matrix A 



This system matrix needs some examination. A crucial step 
is the introduction of submatrices 1. class into the system matrix. 





K 


A ,i_ 




- 


/aI 4 


/A±l 






*•* 


A* 




The asterisk * denotes a 


submatrix 1. class where 




M>NTSI *MPntSZ 


TLn„ 


EE xC' 




EEx2„ 


llxil 






ILxi^ 


Ll^C 


' EXxz;"-" 






H^^Ts 


EZx/,.x2„'^' 


li 


EEx(„x2„ 


ll X/^X2* 


TLxi^yi^'" 






11 






Ely/' 


1 1 ^4 


ll^C<" 




IIKK 


2E<xJ‘ 


E2 




llK<" 


El 


inK^C 



etc , 
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Now it can seem that all submatrices on the same diagonal bands 



2. order 


• are 


identical 


. The submatrices 


II 


At 








A* = 


A, 


- 






A - 




' A 22 






If 


K 


- 4, 


' A, - 




= 


K, 


- <1. 


■K.z. 


r 


A * 

'^L*l 








- 


A* 








“ * Viiw- 



[ 




where 






can be written as 



A., subsequently are renamed 



-- A 


'-A, , 


Z L-f 


/ L 




/f* 

-Az , L 


'Kl., 


ll 



'At 



* 






X2^ 




K- 


V* t 

LLXI,,XZ^ 


EExiIX^ 




^ ^ l(*l 

EE/i,,xz„ 






.2L-< 





S E nX'n!; 



Like the system matrix for the two dimensional approximation this 
submatrix 1. class is not only symmetrical to its main diagonal line 
1, order, but also the elements on each diagonal line 2. order are 
identical* The elements submatrix therefore are 
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renamed to a , where k specifies the place of the submatrix in the 

lx ^ 

:k 

system matrix and i and j define the place of a . . in the submatrix A, . 

k > ij k 











^k;Zf 


" ^k;tZ 




^ k;?/ 


' ^k;zz 


* 


M/ 




















^ /3 

^ O'k; H-Z 3 ' — 

- z - ■■ ■ 

- . /V-T - 





' ^ k; tM 


^ K'l z h-< 


■ ^k;ZH 


O^k ; f- 






- 



where a, ^ can be written as 
k; L 



a 



Kji 



L 1 II 

• ECxC < ? 



L- 



(3.4.) 



Equation (3.4.) therefore is the fundamental step towards the pro- 
gramming of the system matrix and is programmed in REAL FUNCTION S3. 
Since the data XI, X2 and Y are already provided through a COMMON 
block named DATA3, only NPNTSl, NPNTS2, IPOWRl (= k-1) and IP0WR2 
(= L-1) have to be passed to the function that returns the value of 
a through the function name S3. 

K ^ Li 
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* 

SUBROUTINE MAT31 presets a submatrix A in the following way: 

iC 

i) Preset edge section elements (using REAL FUNCTION S3) 



£*■/« 






72yo 



M 

>J>cr//Oc, 



xtrhoi^ 



ii) Copy defined elements diagonally 





^K',3M 





30 







ij 


^«r; ^/y 


^K;ZJ 


^k; a 




Oir,Zn 




^k;JZ 


^K;3S • • . 


'7k-3Af 


"y 


t^;Arrr 


^k,ii 


^Ki2 


'7k-, 13 




^Xill 




'7k-, U 


^K’.in 




%sz 


^k^ii • • • 




t • t 


« • o. 

Ctk-HZ 


'^kjhS 


« • » 

^kj nn 



This subroutine that needs NPNTSl, NPNTS2, M and IPOWRl (== k-1) as 

"k 

input parameters, returns (SUBMl (7 , 7)) to SUBROUTINE MAT3, which 

presets the entire system matrix in the following way 

i) Preset edge section submatrices (using SUBROUTINE MAT31) 



K 



L-H 



5 



I 

A 



Ji-/ 



7^ 






Ah , , 



r 



'f 



vt^ 
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ii) Copy defined submatrices diagonally 



V y 



ZL 






A, A 



A 



i 

n 



A“ 












ZL 



Ku 



B.. 

K 



A, I 

Al 



/A 



/L 






/ 



/At 
At 
• • * 
A* 



A 



A 






'L 

A * 

Ail 



u 



A 



Ll 



The parameters to call this subroutine are NPNTSl and NPNTS2 (// of 

data points), L and M (<the order of the polynominal fit + 1> in XI 

and X2 - direction) and IPRINT. The coefficients C., are returned* 

11 
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3.4.2. Right hand side vector B 



The right hand side vector B will be divided up into L 
subvectors 1. class. 



& - 



B, 



Again, the asterisk (*) denotes a subvector 1. class. Now the 

•k 

subvectors 1. class B, can be written as 

k 



B. 









EEY„-xi„-A2 



Y .... I 



k 

Finally, the elements of the subvectors 1. class can be written as 

L j L- 1 ; ,M 

To compute b ^ > REAL FUNCTION S3 is used. Data XI, X2 and Y are 

available through a COMMON block named DATA3. NPNTSl, NPNTS2, IPOWRl 

(= k-1), IP0WR2 (= L-1) and lY (=1) have to be passed to the function 

and the value of b, is returned through the function name S3. 
k;i. 
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Using FORTRAN programming language, the allocation of b can be 

K. ^ Lt 

performed by two stacked DO loops. 



1=0 

DO 08 L1-1,L,1 
IP0WR1=L1-1 
DO 08 M1=1,M,1 
IP0WR2=M1-1 
1 = 1+1 

08 B(I) = S3(NPNTS1,NPNTS2,IP0WR1,IP0WR2,1) 
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3.5. Software: 



The software to compute the coefficients for a three dimensional 



approximation is described in Appendix A and is implemented in the 
TPL HP-21MX computer system. To work correctly with these program 



modules, the user has to conform to the following conventions: 



i) Provide the data in three arrays (Type: REAL) of 
16 * 16 elements through a COMMON block named 
DATA3. 

COMMON / DATA3 / XI, X2, Y 

REAL XI (16,16), X2(16,16), Y(16,16) 

ii) Dimension an array (Type: REAL) of 7*7 
elements to contain the coefficients. 



REAL COEF (7,7) 



iii) Define the parameters NPNTSl, NPNTS2, L, M and 
IPRINT. 

NPNTSl ... // of XI variations 

L ^ NPNTSl < 16 

NPNTS2 ... // of X2 variations 

M ^ NPNTS2 ^ 16 



L ... (desired order of approximation 

polynominal w.r.t. XI) + 1 

M ... (desired order of approximation 

polynominal w.r.t. X2) + 1 



IPRINT 



controls quantity of print out 
2 ... Print system matrix and 

right hand side vector 
before and after Gauss 
Jordan Elimination 

1 ... Print equation system after 

Gauss Jordan Elimination 
<0 ... No print out 

>0 ... Display equation (3.1.) 

with the actual parameters 
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iv) If a user program uses subroutine MATS, the software 
modules have to be loaded using the procedure out- 
lined in section 2.4. iv. 

If all these requirements are met, the correct call for the subroutine 
is : 

CALL MAT (NPNTS1,NPNTS2,L,M,C0EF,IPRINT) 

Upon completed execution of this approximation routine the array 
COEF contains the coefficients. Externals used by MATS are: ABS, 

DATAS, MATSl, lELS, SS. Under no circumstances may the user use 
any of these names for modules of his own program. 
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3-6 • Sample User Program 

3.6.1, PTN4 Compiler listing of sample user program 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 
0 029 

0030 

0031 

0032 

0033 



0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 
0 053 

0054 

0055 
0 056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0 067 
0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 



FTN4,L 

PROGRAM LAB (3,99) 

*,Plot calibration points and calculate c oef f ic ien t s . 

COMMON / AFLD / PLOTR 
COMMON / DATA3 / X1.X2,Y 
COMMON / AB3 / A,B 

REAL PL0TR(256) 

REAL XK16. 16) .X2(16,16) ,Y(16,16) 

REAL A<49,49) 49) 

REAL YT(16> 16) >R(16) .C0EFF(7 ,7) 

INTEGER II)CB(144) , IFiLE(3) , NOCR ( 2 > , ICLR ( 3 > ,YT1TLE(4) ,2TITLE(4) 

DATA PI /3. 141593/ 

DATA NOLF /006537B/ 

DATA NOCR /O 0 0 033B , 040433B/ 

DATA ICLR /01SS24B, OlSSiSB, 006537B/ 



a R M 

121 FORMAT ( 

122 FORMAT ( 

123 FORMAT ( 
t . . Ueloc 

124 FORMAT ( 



V <*:> l-.< 3 r A R *T* 

“ Enter raw data file nane 
3A2) 

Select 1 . . 

ity X “2A2> 



( 



* = 

*nd Delta 
125 FORMAT ( 
»X( , ) 1 

% C X=X( 

160 FORMAT ( 

131 FORMAT ( 

132 FORMAT ( 

133 FORMAT ( 

134 FORMAT ( 

127 FORMAT ( 

128 FORMAT ( 

129 FORMAT ( 

130 FORMAT ( 

150 FORMAT ( 

151 FORMAT ( 

152 FORMAT ( 

153 FORMAT ( 

154 FORMAT ( 

Max, 

155 FORMAT ( 

Max, 

156 FORMAT ( 
»x, ♦ par 

157 FORMAT ( 
»x, ♦ par 

158 FORMAT ( 

101 FORMAT ( 

102 FORMAT ( 
107 FORMAT ( 
105 FORMAT ( 

^ II o •• 

103 FORMAT ( 
% 

104 FORMAT ( 
149 FORMAT (•• 

611 FORMAT ( 

612 FORMAT ( 

613 FORMAT ( 

614 FORMAT ( 

615 FORMAT ( 

616 FORMAT ( 

617 FORMAT ( 

618 FORMAT ( 

602 FORMAT ( 

603 FORMAT ( 



Select 1 

I. 2 

C = ( , ) > 

•* Select t . . 

“/ 

" 2 . . 

) } -2A2) 



"2A2) 

Pitch angle or 2 . 

Pitch angle v/s Ganna and Beta C 
•• u/s Ganna a 

*2A2) 

Velocity u/s Ganna and Beta C X= 



", IH" , " 



u/s Ganna and Delta 



n in 



*• Enter ♦ of Mach Nunbers and ♦ of pitch angle 

" Enter 1 to print input data "2A2) 

/*• Calibrate = ( , )^V) 

/'• Calibrate = ( , )“/) 

/” Calibrate X=X( , ) ‘ V ) 

/*• Calibrate X=X( , > ”/) 

24X" "9X” ") 

24X" ”9X" “9X" “) 

24X" "9X" ••^’X^X") 

M^^ll 

” nininun ua lues , 3 ( IX , F9 . 3 ) ) 

'* naxinun values" , 3 ( IX , F9 . 3) ) 

" of user coordinates‘'3X"Xl •’8X"X2"9X’‘Y" ) 

'* Enter 1 to redefine these data "2A2) 

“ Enter nin, nax , nin, , 

♦ partitions , 6 ") 

“ Enter nin, nax, nin, ", 

♦ partitions , 6 ") 

" Enter nin, nax, nin, nax, X 

tit ions , 6 X" > 

" Enter nin, nax, nin, nax, X 

tit ions , 4 X" ) 

1X,3(F8.3, 1X,F8.3,5X) ,314) 

" Enter order M of XI - appr o x i na t i on 

" Enter oredr N of X2 - appr o x mat i on 

" Enter IPRINT "2A2) 

" Select 1 ... to calculate the absolute error"/8X 

,4H" ",9X'*the relative error ‘’2A2) 

** Change coefficients ... 1"/ 

" Change axes ... 2 "2A2) 

" Enter ALPHAX, ALPHAY, XO and YO “ ) 

"< (3A2) ) ) 

"1 "////" T* l-i %•> <=> i=> <:! Cl -V Cl r r' o ra 

/16X,I2". Mach nuntaer"/) 

6X"Ganna”6X*’Beta"7X"Ph i " / ) 

6X "Ganna "5 X " De 1 1 a "7X " P h i " / ) 

6X"Ganna"6X"Be ta"5X"X uel "/) 

6X "Ganna "5X"Delta"SX"X uel"/) 

1X,4( 1X,F9.3) ) 

/9X"X1"8X"X2"9X"Y") 

(3X, 16(6X, 12) ) ) 

lX,I2,i6( iX,F7 .3)/3<3X,16( 1X,F7.3)/) ) 



X 

nin , X 

"2A2) 

"2A2) 



na 



na 



3 A 
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0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 



0002 LAB 



3;03 PM FRI . , 26 



SEP., 1980 



604 FORMAT (///“ Absolute error at each point 
♦value - Measured value) }“/) 

644 FORMAT <///** Relative error at each point in X 



C error = (calculated 
C error = (calco 





♦ated value 


- neasur 


606 


FORMAT 


( IX, 


I2.18(iX 


1101 


FORMAT 


(8H 


Beta 


) 


1102 


FORMAT 


(8H 


Delta 


) 


1103 


FORMAT 


(8H 


Phi 


) 


1104 


FORMAT 


(8H 


X vel 


) 



LI = LOGLU(I) 
LO = 6 



C 

C 

c 

c 

c 



read raw data file. 



41 


WRITE 


(LI, 


121) 


NOCR 




READ 


(LI 


122) 


IFILE 




WRITE 


(LI, 


149) 


ICLR 


42 


WRITE 


(LI, 


123) 


NOCR 




READ 


(LI 


♦ ) 


IPX 




WRITE 


(LI 


149) 


ICLR 




IF (IPX.LT, 


. 1 .OR 


. IPX 


43 


CONTINUE 







IF (IPX.EQ.l) WRITE (LI, 124) NOCR 
IF (IPX.EQ.2) WRITE (LI, 125) NOCR 
READ (LI, ♦) IBD 
WRITE (LIl 149) (ICLR, 1 = 1,2, 1) 

IF (IBD.LT.l .OR. IBD.GT.2) GO TO 43 
IRECY=13 
IRECY=9 



IF (IPX.EQ.l) 
■■ :.EQ.2) 



IF (IPX . 

IF (IPX.EQ.l 
IF (IPX.EQ.l 
IF (IPX 



AND. IBD.EQ.l ) WRITE (LI, 131) 



_ .AND. 

_ .EQ.2 .AND. IBD.EQ.l ) WRITE (LI^ 133) 
IF (IPX.EQ.2 .AND. IBD. EQ.2 ) WRITE (Ll! 134) 



IBD. EQ.2 ) WRITE (LI, 132) 



CALL OPEN (IDCB,IERR,IFILE,IOPTN, 0,28, 1^4)' 
IF (lERR .LT. 0) GO TO 41 



CALL READF ( IDCB , lERR . XI .512 , LEN , 1 ) 

IF ( lERR .LT. 'O ) ^TOf^ 0062 
CALL READF ( IDCB , lERR . X2 ,512 . LEN , 5) 

IF ( lERR .LT. 0 ) STOP 0003 

CALL READF ( IDCB , lERR . YT . 5 12 , LEN , IRECY ) 

IF ( lERR .LT. 0 ) STO^> 0004 



CALL CLOSE ( IDCB . lERR , 0 ) 

44 WRITE (LI, 126) NOCR 

READ (LI, ♦) NMACH,NPITCH 
WRITE (LI, 149) ICLR 
IF (NMACH.GT . 16) GO TO 44 
IF (NPITCH.GT. 16) GO TO 44 



DO 801 I=l,NMACH.l 
DO 801 J=1>PITCH,1 
801 Y(I,J)=YT(I,J) 



C 

C 

C 

C 

C 

C 



Correct input data and optional output of data, coefficient 
are based upon. 



47 

48 



IF (IBD.EQ.2) GO TO 48 
DO 47 I=1,NMACH,1 
DO 47 J=1,NPITCH,1 



X2( I, J)=X^(I, J)/il(I, J) 
WRITE (LI, 160) NOCR 



READ (LI, ♦) IDUM 
WRITE (LI. 149) ICLR 
IF (IDUM.NE.l) GO TO 55 
WRITE (LO, 611) IFILE 
DO 54 I=1,NMACH,1 
WRITE (LO, 612) I 

IF (IPX.EQ.l .AND. IBD.EQ.l) WRITE (LO, 613) 



PAGE 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 



0003 LAB 3;03PM FRI., 26 SEP., 1980 



IF (IPX.EQ.l .AND, IBD.EQ.2) WRITE < LO , 614) 

IF (IPX.EQ.2 .AND. TBD.EQ.l) WRITE (LO, 615) 

IF (IPX.EQ.2 .AND. IBD.EQ.2) WRITE ( LO , 616) 

DO 54 J=l,NPITCH,i 

54 WRITE (LO, 617) XI ( I , J ) , X2 ( I , J ) , Y ( I , J ) 

WRITE (LO, 618) ' > > > 



C 

C 

C 

C 

C 

C 



55 



45 



.Calculate nininupi and naximm values for all coordinates and 
optional redefinition of ninirtun and naxifiuM values by user. 

xuMA)(=)<i(i *i) 

XUMIN=Xl(lU) 

YUMAX=X2(1,1) 

YUMIN=X2(1,1) 

ZUMAX=Y(1,1 ) 

ZUhIN=Y(l, 1 ) 

DO 45 I=i,NMACH,l 
45 J=i'nPITCH,1 

(X1(I,J) .GT.XUMAX) XUhAX=Xl(I,J) 

(X1(I,J) .LT.XUNIN) XUhlN=Xl(I, J) 

(X2(I,J) .GT.YUhAX) YUMAX=X2( I , J) 

(X2(I , J) .LT. YUMIN) YUMIN=X2 ( I , J ) 

(Y(I , J) .GT .ZUMAX) ZUMAX=Y(I,J) 

(Y(l' J) .LT.ZUMIN) ZUMIN=Y(lU) 



DO 
IF 
IF 
IF 
IF 
IF 

IF X . vx , 
CONTINUE 
IF (IPX.EQ.l 
IF (IPX.EQ.l 
IF (IPX.EQ.2 
IF (IPX.EQ.2 
WRITE (LI 
WRITE (LI 



WRITE (LI, 
INCX = 4 
INCY - 4 
INCZ = 4 
WRITE (LI, 
READ (LI, 
WRITE (LI, 
IF (I .NE. 1) 



.AND. IBD.EQ.l) WRITE 
.AND. IBD.EQ.2) WRITE 
.AND. IBD.EQ.l) WRITE 
.AND. TBD.£a.2) WRITE 
XUMIN,YUriIN,ZUHIN 
XUMAX,YUhAX, ZUMAX 



150) 
151 ) 
152) 



(LI, 
<LI 
(LI , 
<LI 



127) 

128) 

129) 

130) 



IF 
IF 
IF 
IF 
READ 
WRITE 



(IPX.EQ, 
(IPX.EQ. 
(IPX.EQ, 
(IPX.EQ. 
(LI, 
(Li; 



153) 

*) 

149) 

GO 



AND 

1 .AND 

2 .AND 
2 .AND 

158) 



NOCR 

I 

ICLR 
TO 49 



IBD.EQ.l) 
IBD.EQ.2) 
IBD.EQ. 1) 
IBD.EQ.2) 



WRITE 

WRITE 

WRITE 

WRITE 



(LI, 

(Li; 

(LI 

(LI 



154) 

155) 

156) 

157) 



XUMIN,XUMAX,YUMIN,YUMAX,ZUMIN, ZUMAX, INCX, INCY, INCZ 
XUMIN,XUMAX, YUMIN, YUMAX,ZUMIN, ZUMAX, INCX , INCY , INCZ 



C 

C 

c 

c 

c 

c 
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Initialize plotters define plotter area; calculate scaling 
coefficients XA, XB, YA, YB, ZA and ZB. 



CALL INITG (13) 



XPMIN 


= 


0 . 






XPMAX 




10 . 






YPMIN 


= 


0 . 






YPMAX 


s 


10 . 






ZPMIN 


= 


0 . 






ZPMAX 




10 . 






ALPHAX 


= 


330 . 




• 


ALP HAY 


ss 


30 . 






ALPHAZ 


2T 


90 . 


-xphIn 




XA 


S 


(XPMAX 


)/(XUMAX-XUMIN) 


XB 


= 


(XPMIN^XUMAX 


-XPMAXtXUMIN ) / ( XUMAX-XUMIN ) 


XL 


=: 


(XPMAX 


-XPHIN 


) 


YA 


ss 


(YPMAX 


-YPMIN 


)/(YUMAX-YUMIN) 


YB 


•= 


(YPMIN:<cYUMAX 


-YPMAXttYUMIN)/(YUMAX-YUMIN) 


YL 


s 


(YPMAX 


-YPMIN 


) 


ZA 


ss 


(ZPMAX 


-ZPMIN 


)/(ZUMAX-2UMIN) 
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PAGE 0004 



LAB 



3;03 PM FRI. > 26 SEP 



1980 



0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 
0281 
0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 
J300 

0301 

0302 

0303 

0304 



85 



86 



87 



88 

301 



302 



303 



C 

C 

C 

C 

C 



01 



C 

C 

C 

C 

C 



= (ZPMIN*2UMAX-2PHAX3KZUMIN)/(2UMAX-ZUMIN) 
= (ZPMAX -ZPMIN ) 

= 8. 



ZB 
ZL 

xo 

YO =1 

PL0TR(61) 

PL0TR<62) 

PL0TR(63) 

PL0TR<64) 

PL0TR(65) 

PLGTR(66) 

PL0TR<67) 

PL0TR<68) 

PL0TR(69) 

PLOTR(70) 

PL0TR(71) 

GO TO (86,87) 

CALL CODE 

WRITE ( YTITLE,il01) 

GO TO 88 
CALL CODE 

WRITE <YTITLE .1102) 

GO TO (301,302) IPX 
CALL CODE 

WRITE (ZTITLE, 1103) 

GO TO 303 
CALL CODE 

WRITE (ZTITLE, 1104) 

CALL SETSM (113.1 . ) 

CALL AXIS (XO,Y(D,XL, ALPHAX,6H Ganna. 8 .XUMIN .XUMAX . 4HF6 . 2 , 6 . INCX)^ 
CALL AXIS (XO,YO, YL,ALPHAY, YTITLE,-8,YUMIN,YUMAX,4RF6.2,6,INCY) 
CALL AXIS <X0,Y0,ZL,ALPHAZ,ZTITLE,-8,ZUMIN,ZUMAX,4HF6.2,6,INCZ) 



(ALPHAX*PI )/180 , 
( ALPHAY:<tPI )/180 
(ALPHAZ3KPI)/i80 . 
XO 
YO 
XA 
XB 
YA 
YB 
ZA 
ZB 
IBD 




CALL SETSM (113,2, ) 

DO 01 I=1,NMACH,1 
DO 01 J=1,NPITCH, 1 

CALL THRTU ( XPLOT , YPLOT , XI ( I , J ) . X2 ( I , J ) , Y ( I , J ) ) 
IF ( J .EQ. 1 ) CALL PLOT ( XPLOT , YPL0T,2 ) 

IF ( J ,GT. 1 ) CALL PLOT ( XP LOT , YPLOT , 3 ) 

CONTINUE 

DO 02 J=1 ,NPITCH,1 
DO 02 I=1,NMACH. 1 

CALL THRTU ( XPLOT , YPLOT , XI ( I , J ) . X2 ( I . J ) . Y ( I , J ) ) 
IF ( I .EQ. 1 ) CALL PLOT t XPLOT , YPLOT , 2 ) 

IF ( I ,GT. 1 ) CALL PLOT ( XPLOT , YPLOT , 3 ) 

CONTINUE 



91 




92 



91 



CALL SETSM 
WRITE (LI, 

READ (LI, 

WRITE (LI, 

IF (MORDER. 

WRITE (LI, 

READ " 

WRITE 
WRITE 
READ 
WRITE 

IF (NORDER, 

M=MORDER->*l 
N=NORDER+l 

CALL MAT3 ( NMACH ,NPITCH ,M ,N ,COEFF , IPRINT ) 



(Li; 
(LI, 
<LI, 
(LI , 
(Li; 



( 113,0 . ) 

101) NOCR 
«) MORDER 

149) ICLR 
GT.6) GO TO 

102) NOCR 
NORDER 

149) ICLR 
NOCR 
IPRINT 
ICLR 
GO TO 






107) 

t) 

149) 

GT.6) 
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40 



ii 
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0345 

0346 

0347 

0348 

0349 

0350 

0351 

0352 

0353 

0354 

0355 

0356 

0357 

0358 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 
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0380 



0005 LAB 



3:03 PM FRI., 26 SEP,, 1980 



C 

C 

C 

C 

C 

C 



Overwrite data array with calculated data, based on the 
just obtained coe-Pf icients . 

411 uRiii' ICiy iosi^mcR 

READ (LI, <c) lAR 

URITE <LI, 149) <ICLR,I=1.2) 

IF (lAR.LT.l .OR. IAR.GT.2) GO TO 411 
IF (lAR.EO.l) URITF. (LO, 604) 

IF (IAR.EQ.2) write (LOJ 644) 

WRITE (LO, 602) ( J , J=1 , NP ITCH , 1 ) 

DO 20 I=l,NMACH,i 
DO 19 J=1 ,NPITCH,i 
SUM=0. 

DO 96 11 = 1, M,i 
IEXPi=Il-i 

IF ( lEXPi .EQ. 0) GO TO 401 
X1EXP=X1( I, J)«¥IEXP1 
GO TO 402 

401 XlEXP=i. 

402 CONTINUE 

DO 96 Jl=i,N,i 
IEXP2=Ji-i 

IF (IEXP2.EQ.0) GO TO 403 
X2EXP=X2( I,J)*)riEXP2 
GO TO 404 

403 X2EXP=i. 

404 CONTINUE 

96 SUM=SUM+COEFF ( II , J 1 ) *X1EXP*X2EXP 
GO TO (412,413) lAR 

412 R<J)=SUM-Y<I,J) 

GO TO 19 

413 R( J)=< (SUM-Y(I,J))/Y(I,J) )«100 . 

19 Y(I,J)=SUM 

20 WRITE (6,603) I , ( R ( J ) , J=1 , NP ITCH , 1 ) 

WRITE (6,602) < J , J=1 ,NPITCH, 1 ) 



C 

C 

C 

C 




CALL SETSM (113,3. ) 

DO 81 I=i,NMACH.l 
DO 81 J=l,NPITCH,i 

CALL THRTU ( XPLOT > YPLOT .XI ( I , J ) > X2< I , J ) , Y ( I , J ) ) 
IF (J.EQ.i) CALL f^LOT (iPLOT, YPdOT,25 
IF (J.GT.l) CALL PLOT ( XPLOT , YPLOT ,3 ) 

81 CONTINUE 

DO 82 J=i,NPITCH,l 
DO 82 I=l,NHACH.i 

CALL THRTU ( XPLOT . YPLOT . XI ( I , J ) . X2( I , J ) , Y < I , J ) ) 
IF (I.EQ.l) CALL I^LQT ( )!:PL0T , YPdOT , 2 i 
IF (I.GT.i) CALL PLOT ( XPLOT , YPLOT ,3 ) 

82 CONTINUE 



C 

C 

C 

C 

C 




DO 802 I=l,NMACH,i 
DO 802 J=i'nPITCH 
802 Y(I, J)=YT(I, J) 
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0381 

0382 

0383 

0384 

0385 

0386 

0387 

0388 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 



C 

C 

C 

C 



Next step ? 



WRITE <LI,i03) NOCR 
READ (LI, t) IDUM 
WRITE (Ll! 149) (ICLR .1 = 1,2,!) 

IF (IDUM.EQ.l) GO TO 91 
IF (IDUM.EQ.2) GO TO 89 
CALL STOPG 
STOP 077 

89 WRITE (LI, 104) 

READ (LI, ALPHAX,ALPHAY,XO,YO 

WRITE (LI, 149) ICLR 
GO TO 85 
END 



FTN4 COMPILER: HP92060-16092 REU. 1926 (790430) 



tt NO WARNINGS tt NO ERRORS PROGRAM = 04465 



COMMON 



00000 
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3.6.2. Load map 



LAB 



10043 20633 Plot calibration ooinT<=> and cnlrulatp c f ic lents 



to load these program modules, enter (from LOADR): MS , XTPLBL 



DATA3 

MAT3 

MAT31 

S3 

AB3 



20623 23622 
23623 27455 
27456 27651 



27652 

30071 



30070 

41534 



3D-Appr 0 X inat ion' / DATA3 / 

3D-Appr o X inat 1 on arrange systen natrix and vector, 
3D-Appr 0 X i«at i on arrange subnatrix 1. clas*^, 
3D-Appr 0 X inat 1 on conpute sunnations. 
3D-Aoor 0 X i«at i on / AB3 / 



AFLD 41535 42534 GSP / AFLD / A<256) 

AXIS 42535 45054 GSP Draw and label axes. 

INITG 45055 46067 GSP read control array Aj initialiie plotter. 

PLOT 46070 46620 GSP nove pen to a defined point. 

SETSh 46621 47561 GSP change plotter nodes. 

STOPG 47562 47624 GSP terninate grafics. 

THRTW 47625 50111 GSP converts 3D user coordinates to plotter units, 
FCTR 50112 50115 GSP / FCTR / FX,FY 

FACTR 50116 50135 GSP vary size of the plot. 

LOGLU 50136 50213 92067-16268 REO.1903 790228 

READF 50214 51155 92067-16125 REO.1940 790719 

OPEN 51156 51453 92067-16125 REV. 1903 790215 

CLOSE 51454 51663 92067-16125 REV. 1903 791229 

CLRIO 51664 51672 750701 24998-16001 

OVRD. 51673 51673 92067-16125 REV. 1903 780526 

SSMVE 51674 51762 92067-16268 REV. 1903 790202 

LURQ 51763 52345 92067-16268 REV. 1903 790223 

.DADS 52346 52455 780818 24998-16001 

.DhP 52456 52623 780818 24998-16001 

.DDI 52624 53124 781021 24998-16001 

SESSN 53125 53142 92067-16125 REV. 1903 780413 

R/U% 53143 53301 92067-16125 REV. 1903 791214 

P.PAS 53302 53330 92067-16125 REV. 1903 740801 

.DNG 53331 53340 780818 24998-16001 

PAUSE 53341 53441 771122 24998-16001 

SALRN 53442 53557 92067-16268 REV. 1903 770715 

FMTIO 53560 55056 24998-16002 REV . 1926 790417 

ERRO 55057 55146 771122 24998-16001 

TAN 55147 55253 780424 24998-16001 

ABS 55254 55262 750701 24998-16001 

.SNCS 55263 55424 780424 24998-16001 

.DDE 55425 55436 780818 24998-16001 

.DIN 55437 55444 780818 24998-16001 

.RTOI 55445 55540 780921 24998-16001 

.FPUR 55541 55602 781106 24998-16001 

.SBT 55603 55643 770518 24998-16001 

..FCM 55644 55660 750701 24998-16001 

PAU.E 55661 55661 750701 24998-16001 

ERO.E 55662 55662 750701 24998-16001 

.CMRS 55663 55746 780424 24998-16001 

iOPEN 55747 56123 92067-16125 REV. 1903 790103 

RUSUB 56124 56471 92067-16125 REV. 1903 781003 

RWNDS 56472 56614 92067-16125 REV. 1903 780801 

FRMTR 56615 62252 24998-16002 REV. 1926 790503 

FMT.E 62253 62253 24998-16002 REV. 1901 781107 

IB4A2 62254 63554 59310-1X013 REV 1940 790802 1153 

REIO 63555 63701 92067-16268 REV. 1903 790316 

RMPAR 63702 63744 781106 24998-16001 

PNAHE 63745 64012 771121 24998-16001 

LUTRU 64013 64121 92067-16268 REV. 1903 790223 

IPUT 64122 64142 92067-16125 REV. 1903 740801 

SSETP 64143 64167 781106 24998-16001 

.CFER 64170 64245 750701 24998-16001 

.LBT 64246 64276 770518 24998-16001 
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/LOADR : SEND 



NO PAGES MSEG 



43 



3 . 6.3 



Results 



i) printed output 



L* 1;> kK> 


c:l cx V cx f 


' r' o 




1. Mach 


nunber 


Ganna 


Beta 


Phi 


-.810 


-.005 


-35.000 


-.519 


-.005 


-25.000 


-.280 


-.005 


-15.000 


-.182 


-.005 


-10 . 000 


-.103 


-.006 


-5.000 


-. 0S9 


-. 006 


-2.500 


-.018 


-.006 


0 .000 


.019 


-. 006 


2.500 


.054 


-.006 


5.000 


. 133 


-.006 


10 . 000 


. 197 


- .006 


15.000 


.360 


-. 006 


25. 000 


.605 


-.006 


35.000 




2. Mach 


nunber 


Ganna 


Beta 


Phi 


-.781 


-.016 


-35.000 


-.506 


-.016 


-25.000 


-.280 


-.016 


-15.000 


-.178 


-.017 


-10.000 


-.092 


-.018 


-5.000 


-.048 


-.018 


-2.500 


-. 004 


-. 018 


0 . 000 
2.500 


. 032 


-.018 


. 061 
. 130 


-.018 


5.000 


-.019 


10.000 


. 190 
.346 


-. 019 


IS. 000 


-.020 


25.000 


.583 


-.018 


35. 000 




3. Mach 


nunber 
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Beta 


Phi 


-.787 


-.024 


-35.000 


-.510 


-.024 


-25. 000 


-.276 


-.024 


-15.000 


-.172 


-. 026 


-10 . 000 


-.085 


-.027 


-5.000 


-.050 


-.027 


-2.500 


-. 003 


-.027 


0.000 


. 037 


-.027 


2.500 


. 066 


-.027 


5.000 


. 131 
. 191 
.345 


-.027 

-.029 


10 . 000 
15.000 


-.029 


25.000 


.592 


-.027 


35.000 




4. Mach 


nunber 


Ganna 


Beta 


Phi 


-.790 


-. 035 


-35. 000 


-.519 


-.035 


-25.000 


-.276 


-.035 


-15.000 


-.174 


-.037 


-10.000 


-.055 


-. 039 


-5.000 


-. 043 


-.039 


-2.500 


-.003 


-.038 


0.000 


. 037 


-.038 


2.500 


. 064 


-. 039 


5.000 


. 133 


-.039 


10.000 


. 196 


-. 041 


15. 000 


.349 


-.042 


25.000 


.596 


-. 039 


35. 000 




5. Mach 


n unber 



:i: A A . 
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Ganna 


Beta 


Phi 


-.795 


-,045 


-35.000 

-25.000 


-.521 


-.045 


- . 282 


-.045 


-15.000 


-.175 


-. 048 


-10.000 


-.082 


-.050 


-5.000 


-.040 


-, 050 


-2.500 


-.010 


-.050 


0.000 


. 034 


-.049 


2.500 


. 067 


-.049 


5.000 


. 134 


-.050 


10.000 


. 199 


-.053 


15.000 


.350 


-.055 


25. 000 


.593 


-.050 


35.000 
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-.065 
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-.068 


-10.000 


-.085 


-. 071 


-5.000 


-.032 


-.072 


-2.500 


-.004 


-.071 


0 . 000 


. 024 


-.071 


2.500 


. 058 


-. 071 


5. 000 


.131 


-.071 


10.000 


. 199 


-. 075 


15.000 


,349 


-.078 


25.000 

35. 000 


.603 


-.073 


XI 


X2 


Y 
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ii) graphic output 
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4. FOUR DIMENSIONAL APPROXIMATION 



4.1. Problem: 




A data set of NPNTSl ^ NPNTS2 * NPNTS3 data points is given, 
where Y depends on the parameters XI, X2 and X3, The data pattern is 
to be approximated by a function Y = f(Xl,X2,X3), so that the error 
between data points and analytically determined points is lower. 
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4.2. Approach 



Again, as in the two and three dimensional case, a polynoininal 
method is used. For 

V". fcM, i-Z, xz) 

we assume 

K- {[€,„ . ... ^ 

» • ♦ 

Af-/ j 

.... ^ 1-xz ]t 

/^A-/ 

I^'^ZZ, iC^XZi. .. ^ X3 ] -Xl * 



N-t 



* C„,XS- ... * C,js 1 

" i ^ 

"Ul ^ ^LZ2 



■( 
A/-/. 



-> C c,,, i C,„y3 ■! ... C^^t^XS J-xZ 



^ ^ ^LHi ••• C-i^f^XZ ] " XZ fx/ 



or 



/= E [E fE C;.«''}Aa''''} Xl'"' 

** 



(4.1.) 
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This data arrangement which discriminates between data points taken at 



constant XI, X2 and X3, respectively is needed in the application of 
this method to the description of pneumatic-velocity probe characteristics. 
The various surfaces in the above sketch for example can be seen as data 
taken at different Mach number (XI) . At each constant Mach number 
(XI = constant) yaw angle (X2) and pitch angle (X3) are varied, and the 
probe is measured to be Y. Even though probes are usually balanced to 
the average flow yaw angle in some cases this is not possible and the 
four dimensional approximation is then needed. A similar method for 
describing probe calibration data is used at the institute for Jet 
Engines and Turbomachines at the Technical University Aachen (Ref. 2). 

The least squares criterion is then 




where the indices rst denote data points. Using equation (4.1.), R 



becomes 








r«i 
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4 - 



V~f 

C Cm + ^iZZ^^rst* '••'*’ ^iZH ^ 

« « • 

The term in the [ ] - bracket shall be called B 

>^><<TS/ NTA/Ki 

^-E{2fE fi'i{ 

r»/ t*i t-*/ 
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4,3. Solution: 



To minimize the error, R is partially differentiated with respect 

to the coefficients C.., and the partial derivatives are set to zero, 

ijk 

thus ; 



bC-k 






/. / ... i . j.. I ... M ^ U- / ... ^ 



'^P NfNKi p 

EfLiE o 

r»i -»«/ Uf 






3 c 



Assuming, that the summations extend over all data points, ZZZ should 
be understood to be Xj Xj YJ ' even a very large sheet 



r*/ 






notation. 



can* t 
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system matrix. 


unless 


we sk: 
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... A* 








S* 
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« o « 














K 


K. 


. . . 




-Jt 1 




s! ' 



(4.3.) 



The asterisk indicates a submatrix 1. class. On the following pages 
the submatrices and subvectors 1. class are listed. 
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4.4. Structure of Equatiog System ; 



4.4.1. Svstem Matrix A 



It can be seen, that all submatrices 1. class on the same 



diagonal bands 2. order of the system matrix are identical and therefore 

* * 



will be renamed from to A^ 



A*’ A* 

4 ■ 4 ■■ K 

4;. < .4; 



A 



L*f 



■ 4.Z 



- ^ - 

- ^L-Z,Z - ■■ 


■ ■ 


- A* 


-K;z- ■■ 


^ A* 

L-t 


- 


■K,, ' ... 


-yA * 





K 

A-/ 



- 

'^L , 1 , 



The submatrices 1. class A, (k=l, • • . ,2L-1) themselves are divided 
up into submacrices 2. class. 





K-i,> 



A 






A 



kjlZ 



1,'* 



4 A 






k)m 



A 



w-r 

k;n,z 



/A^* 
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The double asterisk indicates a submatrix 2, class. All submatrices 

2. class on the same diagonal bands 2. order of the submatrices 1. class 
* 

(k = 1,...,2L) are identical and therefore will be renamed from 




to 






1 



A 



a 






= A 



* * 

k,i,i 





A** 

^ 2,1 


hZ 










A** 


- A** 


' u 








A** 


A** 

“ A 2 


••• 


‘ A**' 








-A** 


-A*" 


- A** 


- 

- 


K'h,^ 








-- A^^.. f fi.i 


A 

-- 



A 



« *■ 

k- 2M-! 



where 




can be written as 



-A 



it* 

k; !,fi 



A 



k;L 



EEExi^-xz':; 








xzM-xi:, ... 


EzaiX-x3i 


ZECxC«‘>: 


• « « « « • 

EEE)!1"x2;;|xJ^ ... 


• • • 

EEZxi^.xz^'C 


k- !,..., ZL-i 
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Note, that this submatrix 2, class is not only symmetrical, but also the 

elements on each diagonal line 2, order are identical. The elements 

a, , , . of Che submatrix 2. class A, therefore are renamed to a, _ ; 

k; 1 k;L;m’ 

L specifies the place of the submatrix 2. class A, ^ in the submatrix 

^ > 1 

* -k 

1. class A^, k specifies the place of the submatrix 1. class in the 

system matrix A and i,j specifies the place of the element a. . . in 

A > i * J 



the submatrix 2. class A, 









Cl ^ CL 



%l;Z. *' *' ^KU;I.3 



^y-Li3 ‘ 





















= a 









U-L-ZH-i 



-■ Q. 



k ;L; H.n 



where 






can be written as 









r-\ S-t t'l 



(4.4.) 



/, .. ,2L-( y j M. /, 
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This operation is programmed in REAL FUNCTION 34 (NPNTS1,NPNTS2 ,NPNTS3 , 
IPOWRl,IPOWR2,IPOWR3,IY) . The data XI, X2, X3 and Y are known to 
this function through a common block named DATA4. So only NPNTSl, 
NPNTS2, NPNTS3, IPOWRl (=» k-1), IP0WR2 (=» L-1), IP0WR3 (=» m-1) and 



lY (=0) have to be passed to the function and the value of 
returned through the function name S4. 

SUBROUTINE MAT42 presets the submatrices 2. class A, in the 
following way 

i) Preset edge section elements (using REAL FUNCTION S4) 









^X;ii W 

%L; 






^yecno^ 






ii) Copy defined elements diagonally 






I,M 



X}i;N,N 
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^k;l; V,A/ 





^k;L,l,Z 




^kXl.H 
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' ‘ ' 




%L;3,!^ 
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‘■kfhHN 
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%L,2.^ 


^;hH 


%l;2Z 
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9 • • 
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, • • • 

^:hu 


^k;L,i,N 




%1;12 




2,V 


%IA' 


^k-MZ 


%iiZt ••• 




• • • 


m • • 


• « • m • • 






°-k;l;N.Z 




^•4- 



71 



This subroutine that needs NPNTSl,>rPNTS2,OTNTS3,IPOWRl(=k-l) and 
IP0WR2(=L-1) as input parameters, returns the submatrix 2. class 
A, (SUBM2(4,4)) to SUBROUTINE MAT41, which foresets the submatrices 

ilC 9 JL 

* 

1. class in the following way 

i) Preset edge section submatrices (using SUBROUTINE MAT42) 






A. 

•O.- 












A 



** 






hi 









ii) Copy defined submatrices diagonally 



A 

A 



k; 2. k 

** 

k;lfi 



.44 
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A** 




A* * 




^U;22 


A:2rJ 










/ * 


• * * ^ 


/* 
m m ^ 










At* 

* k ,M,M 






“k;U 


• . • 


\.,r. 






A^* 




.4-# 




K%,z 


A** 




A* * 








. • • 

^‘‘.nn 






c 








A** 






KIm 


A** 


A*"* 

Ak.j.z 

/ 


^k;Jj 




J* 

nk;3,H 


A** 


4^* 


• • • 

A ^ 


• • • 





This subroutine, that needs NPNTS1,NPNTS2,NPNTS3,M,N,IP0WR1(-L-1) 

* 

as input parameters, returns the subroutine 1st class (SUBMl (16,16)) to 
SUBROUTINE MAT4, which presets the entire system matrix A in the 
following way. 



73 



i) Preset edge section submatrices (using SUBROUTINE MAT41) 



.. 



. =0 

A., 






& 

^iL-t 




ii) Copy defined submatrices diagonally 






ZL 



A 



LL 



K 



A 






Al 






^ ‘y ^ 



'JL 



/A 



ji. 
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4.4. 2, Right Hand Side Vector B ; 



ik 

The right hand side subvectors 1. class generally can be 



written as 




I S 

• • • 

Z S 

• • • 

EE2y„^xi„^;'x2„v<* 



E S EX.-X1HX2:: 
ESEXfX|H<:'x?^ 

I ESsr^-xi^'<-<' 

k- ... . i 
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The subvectors 1, class (k 
up into subvectors 2. class: 



= L) themselves will be divided 









*¥■ 



c 



** 



where subvectors 2. class B, can be written as 

iC 9 






H* 



ZSEVxi|;;'-x2‘^x? 






k’’ I J ..., L ; , M ’ rr,^ !j ... , N 



Finally, the elements of the subvectors 2. class B can be written as 

k ] 1 






k) Lj n% 



k--l, ... jM ; m-- k ••• y ^ 



To calculate b, , REAL FUNCTION S4 is used, again. Data XI, X2, X3 

k ^ X 

and Y are available through COMMON block DATA4. NPNTSl, NPNTS2, NPNTS3, 

IPOWRl (=• k-1), IP0WR2 (=• L-1) IP0WR3 (=• m-1) and lY (=1) have to be 

passed to the ftinction and the value of b, , . is returned through the 

k ^ X , m 

function name S4. Using FORTRAN programming language, the allocation 

of b, can be performed by three stacked DO-loops. 

k 5 1 ^m 
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1=0 

DO 07 L1=L,L,1 
IP0WR1=L1-1 
DO 08 M1=1,M,1 
IP0WR2=M1-1 
DO 08 N1=1,N,1 
IP0^^R3=N1-1 
1 = 1+1 

08 3 (I)=SA (NPNTSl ,NPNTS2 , NPNTS3 , IPOWRl , IPOWR2 , IPCWR3 , 1) 
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4,5* Software; 



The software to compute the coefficients for a two dimensional 
approximation is described in APPENDIX A and is implemented in the TPL 
HP*21MX computer system. To work correctly with these program modules. 



the user has to conform to the following conventions. 



i) Provide the data in four arrays 
(Type: REAL) of 5*5*5* elements 
through a COMMON block, named DATA4. 

COMMON / DATA4 / X1,X2,X3,Y 

REAL X1(5,5,5),X2(5,5,5),X3(5,5,5),Y(5,5,5) 

ii) Dimension an array (Type: REAL) of 

5*5*5 elements to contain the coefficients, 

REAL COEF (5,5,5) 

iii) Define the parameters NPNTSl, NPNTS2, 



NPNTS3,L,M,N and IPRINT 



NPNTSl. . . 
NPNTS2. . . 
NPNTS3... 



// of XI variations 
L ^ NPNTSl ^ 5 

f/ of X2 variations 
M < NPNTS2 ^ 5 

// of X3 variations 
N NPNTS3 ^ 5 



L 



M 



N 



(desired order of approximation 
polynominal w.r.t. XI) + 1 
1 < L < 4 

(desired order of approximation 
polynominal w.r.t. X2) + 1 
1 M ^ 4 

(desired order of approximation 
polynominal w.r.t. X3) + 1 
i ^ N ^ 4 
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IPRINT . . . Controls quantity of print out 

2 ... Print system matrix and right 
hand side vector before and 
after Gauss Jordan Elimination. 

1 ... Print equation system after 
Gauss Jordan Elimination 
<0... No print out 
>0... Print equation (4.1.) with 
the actual parameters. 

iv) If a user program uses subroutine MAT4, the software 
modules have to be loaded using the procedure in 
section 1.4. iv. 

If all these requirements are met, the correct call for the subroutine is 
CALL MAT4 (NPNTSl ,NPNTS2 ,NPNTS3 , L,M, N, COEF , IPRINT) 



Upon completed execution of this approximation routine array COEF 
contains the coefficients. Externals used by MAT4: AB4, DATA4, MAT41, 

MAT42, S4. Under no circumstances may the user use any of these names fo 
a modules of his own user program. 
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4.6, 

4.6. 
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0002 
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0007 
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0016 
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0019 

0020 
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0023 

0024 

0025 

0026 
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0028 

0029 
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0033 

0034 

0035 
0 036 
0037 
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0059 
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0065 
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0070 

0071 
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0074 

0075 

0076 



Sample User Program 

1 FTN4 Compiler of sample user program 

FTN4,L 



PROGRAM FITR (3,99) 

APPROXIMATION PROGRAM TO FIT A SURFACE THROUGH DATA POINTS Y 
WHICH DEPEND ON THREE INDEPENDENT PARAMETERS XI, X2 AND X3 . 

AUTHOR i HANS 2EBNER 
DATE : AUGUST 26, 1980 



Hans Zebneri 



40-Approxi«at ion j Y*f <X1 ,X2,X3) . 



C 

C 

C 



COMMON / DATA4 / Xi.X2,X3,Y 
COMMON / AB4 / A,B 

REAL Xl(5,5,5) ,X2<5,5,5) ,X3( 5 , 5 , 5) , Y (5 ,5 ,5) 

REAL A(64,6^'> 

REAL CQEF<S,5,S) 

INTEGER IOCB(U4) ,1F1LE<3) ,NOCR ( 2) , ICLR <35 

DATA NOCR /0000333, 040433B/, ICLR /015S24B, 01SS15B, 006S37B/ 

101 FORMAT (• ENTER NPNTSl NPNTS2 NPNTS3 *2A2 

«) 

102 FORMAT <“ ENTER L M N *2A2) 

103 FORMAT (* ENTER DATA FILE NAME *2A2) 

104 FORMAT <3A2) 

131 FORMAT <* ENTER IPRINT •2A2) 

132 FORMAT C OUTPUT RAU DATA YES OR NO “2A2) 

105 FORMAT </• XI*) 

106 FORMAT </lX,I2*. DATA SET"/) 

107 FORMAT <3X, 12110) 

108 FORMAT ( IX, 12, 14F10 . 3) 

109 FORMAT </“ X2*) 

110 FORMAT </“ X3“) 

111 FORMAT </" Y*) 

149 FORMAT <*“(<3A2))) 

1101 FORMAT ("HANNES") 

READ RAU DATA FROM FILE 

01 URITE (1,101) NOCR 

READ <1, ») NPNTSl, NPNTS2,NPNTS3 

URITE (1,149) ICLR 
NMAX=S 

IF (NPNTSl .GT .NMAX) GO TO 01 
IF (NPNTS2.GT.NMAX) GO TO 01 
IF (NPNTS3.GT.NMAX) GO TO 01 

02 URITE (1,102) NOCR 

READ (1, *) L.M.N 

URITE (1,149) ICLR 
IF (L»M*N.GT.S4) GO TO 02 

03 URITE (1,103) NOCR 
READ (1,104) IFILE 
URITE (1,149) ICLR 

31 URITE (1,131) NOCR 
(1, *) IPRINT 



TE (1,149) ICLR 
TE (1,132) NOCR 



READ 
UR IT 
URIT 

READ (1,104) IDUM 
URITE (1,149) ICLR 
CALL OPEN (IDC8.IERR, IFILE) 

IF (lERR.LT.O) GO TO 03 

IL=»2«NPNTS1*NPNTS2*NPNTS3 

CALL READF ( IDCS, lERR ,X1 , IL, LEN, 1 ) 

IF (lERR.LT.O) STOP OOOl 

CALL READF ( IDCB , lERR , X2, IL, LEN , 54) 

IF (lERR.LT.O) STOP 0002 

CALL READF ( IDCB ,IERR .X3., IL, LEN , 107) 

IF (lERR.LT.O) STOP 0d03 

CALL READF ( IDCB ,IERR,Y ,IL, LEN, 160) 

IF (lERR.LT.O) STOP O0O4 

CALL CLOSE (IDCB.IERR) 

IF (lERR.LT.O) STOP OOOS 
IF (IDUM.NE.3HYE) GO TO 17 
URITE (6,105) 
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PAGE 0002 FITR 3 ; 47 PM FRI., 26 SEP., 1980 

0077 DO 04 Kl=l .NPNTSl, 1 

0078 URITE (6>i66) K1 

0079 URITE (6,107) < K3 > K3=»l , NPNTS3 , 1) 

0080 DO 04 K2=1,NPNTS2, 1 

0081 04 URITE (6,108) K2 , ( XI ( K 1 ,K2 , K3 ) ,K3»1 , NPNTS3 , 

0082 URITE (6,109) 

0083 DO OS Kl=l .NPNTSl, 1 

0084 URITE (6,106) K1 

0085 URITE (6,107) ( K3 , K3=l , NPNTS3 , 1 ) 

0086 DO 05 K2=l .NPNTS2.1 

0087 05 URITE (6,1(58) K2 , i X2(K 1 ,K2 , K3) ,K3»1 ,NPNTS3 ^ 

0088 URITE (6,110) 

0089 DO 06 Kl=^l .NPNTSl,! 

0090 WRITE (6,106) K1 

0091 URITE (6,107) ( K3 , K3»l , NPNTS3 , 1 ) 

0092 DO 06 K2»1,NPNTS2,1 

0093 06 URITE (6,l08) K2, <X3(K 1 ,K2,K3) ,K3»1 ,NPNTS3 ; 

0094 URITE (6,111) 

0095 DO 07 Klal .NPNTSl,! 

0096 WRITE (6,106) K1 

0097 WRITE <6,107) (K3 ,K3«1 ,NPNTS3 , 1 ) 

0098 DO 07 K2=l .NPNTS2,1 

0099 07 WRITE (6,l08) K2, ( Y(K 1 ,K2 .K3) .KS^l .NPNTS3 , 1 

0100 17 CALL MAT4 (NPNTSi ,NPNTS2,NPNTS3,L,M,N,C0EF, 

0101 CALL CODE 

0102 WRITE (IFIL£>1101) 

0103 CALL OPEN ( IDCB > lERR . IFILE , lOPTN , 0 ,26 , 144 ) 

0104 IF (lERR.LT.O) STOP ()006 

0105 CALL URITF C IDCB , lERR . COEF , 250 ) 

0106 IF <IERR.LT,0) STOP 0007 

0107 CALL CLOSE (IDCB.IERR) 

0108 IF (lERR.LT.O) STOP OOli 

0109 STOP 0077 

0110 END 



FTN4 COMPILER: HP92060-16092 REV). 1926 (790430) 



1 ) 



1 ) 



1 ) 



) 

IPRINT) 



NO WARNINGS «« NO ERRORS 



PROGRAM a 01259 



COMMON 3 00000 



^,6.2* Load mao of sample user prog:ram FITH 



FITR 10042 12414 Hans Zebn^r: 40-Aoor o x i na 1 1 on ; Y=f < X 1 , X2 » X3 ) . 



to load these program nodules, enter (from LOADR): MS , ^TPL3L 



DATA4 1241S 14364 
MAT4 1436S 202S1 
MAT41 202S2 20666 
MAT42 20667 21066 
S4 21067 21370 
AB4 21371 41S70 



4D-Aopr ox i«at ion / DATA4 / 

4D-Appr 0 X ifiat 1 on arrange systen natrix and vector. 
4D-Approx i«at ion arrange suPwatrix 1. class. 
4D-Appr ox i^^at ion arrange submatrix 2. class. 
4D-Appr ox mat ion conpute sumnat i ons . 
4D-Appr o X inat ion / AB4 / 



LOGLU 

READF 

OPEN 

CLOSE 

CLRIO 

OURD. 

SSMME 

LURQ 

.DADS 

.DMP 

.DDI 

SESSN 

R/WS 

P.PAS 

.DNG 

PAUSE 

SALRN 

FMTIO 

ERRO 

ABS 

.DDE 

.DIN 

.RTOI 

.FPUR 

.SBT 

. .FCM 

PAU.e 

ERO .£ 

SOPEN 

RWSUB 

RUND4 

FRMTR 

FHT.E 

REIO 

RMPAR 

PNAME 

LUTRU 

SSETP 

.CFER 

.LBT 



41S71 

41647 

42611 

43107 

43317 

43326 

43327 
43416 
44001 
44111 
442S7 
44560 
44S76 
44735 
44764 
44774 
45075 
45213 
46512 
46602 
46611 
46623 
46631 
46725 
46767 
47030 

47045 

47046 

47047 
47224 
47572 
47715 

53353 

53354 
53501 
53544 
53612 
53721 
53746 
54024 



41646 

42610 

43106 

43316 

43325 

43326 
43415 
44000 
44110 
44256 
44557 
44575 
44734 
44763 
44773 
45074 
45212 
46511 
46601 
46610 

466i? 

46724 

46766 

47027 

47044 

47045 

47046 
47223 
47571 
47714 

53352 

53353 
53500 
53543 
53611 
53720 
53745 
54023 
54054 



92067-16268 REV. 1903 790228 
92067-16125 REV. 1940 790719 
92067-16125 REV. 1903 790215 
92067-16125 REV. 1903 781229 
750701 24998-16001 

92067-16125 REV. 1903 780526 
92067-16268 REV. 1903 790202 
92067-16260 REV. 1903 790223 
780810 24998-16001 

780818 24998-16001 

781021 24998-16001 

92067-16125 REV. 1903 780413 
92067-16125 REV. 1903 781214 
92067-16125 REV. 1903 740801 
780818 24998-16001 

771122 24998-16001 

92067-16268 REV. 1903 770715 
24998-16002 REV. 1926 79Q417 
771122 24998-16001 

24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
92067-16125 REV. 1903 790103 
92067-16125 REV. 1903 781003 
92067-16125 REV. 1903 780801 
24998-16002 REV. 1926 790503 
24998-16002 REV. 1901 781107 
92067-16268 REV. 1903 790316 
781106 24998-16001 

771121 24998-16001 

92067-16268 REV. 1903 790223 
731106 24990-16001 

750701 24999-16001 

770518 24998-16001 



750701 
780818 
780818 
780921 
781106 
770518 
750701 
750701 
750701 



20 PAGES RELOCATED 20 PAGES REQ^D NO PAGES EMA 
LINKS? BP PROGRAM: LB LOADrTE COMMON : NC 
/LOADR: FITR READY AT 3:53 PM FRl., 26 SEPT, 1980 

/LOADR : SEND 



NO PAGES MSEC 
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i/i>urOr^ 04 rv) ^ fo oi^oiroK^ ui ui<^04roM^ (/i4^oiroi^ 04 i/i>G 4T^»-*^ fu i/i^uron^ 



4.6*3. Results 

printed output 



DATA SET 










1 


2 


3 


4 


5 


.100 


.100 


. 100 


.100 


.100 


. 100 


. 100 


. 100 


. 100 


. 100 


.100 


. 100 


. 100 


.100 


. 100 


. 100 


.100 


.100 


. 100 


. 100 


.100 


.100 


. 100 


.100 


.100 


DATA SET 










1 


2 


3 


4 


5 


.200 


.200 


.200 


.200 


.200 


.200 

.200 


.20 0 
.200 


.200 

.200 


-.m 


.200 

.200 


.200 


.200 


.200 


.200 


.200 


.200 


.200 


.200 


.200 


.200 


DATA SET 










1 


2 


3 


4 


5 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


.300 


DATA SET 
1 


2 


3 


4 


5 


.400 


.40 0 


.400 


.400 


.400 


.400 


.400 


.400 


. 400 


.400 


.400 


. 400 


.400 


.400 


.400 


.400 

.400 


. 40 Q 
.400 




-.m 


•M 


DATA SET 










1 


2 


3 


4 


5 


.SCO 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 


.500 



DATA SET 


1 


2 


3 


4 


5 


- 20.000 


- 20 . 000 


- 20.000 


- 20.000 


- 20.000 


- 10.000 


- 10.000 


- 10 .000 


- 10 . 000 


- 10 . 000 


0.000 


0.000 


0. 000 


0. 000 


0. 000 


10 . 000 


10.000 


10.000 


10 . 000 


10.000 


20.000 


20 . 000 


20.000 


20.000 


20.000 


DATA SET 


1 


2 


3 


4 


5 


- 20.000 


- 20.000 


- 20.000 


- 20.000 


- 20.000 


- 10.000 


- 10.000 


- 10.000 


- 10.000 


- 10.000 


0.000 


0.000 


0.000 


0 . 000 


0 .000 


10.000 


10 .000 


10.000 


10 . 000 


10 . 000 


20.000 
DATA SET 


20.000 


20.000 


20.000 


20 . 000 


1 


2 


3 


4 


5 


- 20.000 


- 20 . 000 


- 20.000 


- 20.000 


- 20.000 


- 10 .000 


- 10 . 000 


- 10 ,000 


- 10.000 


- 10 . 000 


0 . 000 


0 . 000 


0.000 


0 . 000 


0 . 000 


10.000 


10.000 

20.000 


10.000 


10.000 


10.000 


20.000 


20.000 


20.000 


20.000 


DATA SET 


1 


2 


3 


4 


5 


. . 


— 
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— 


. 



in ^ 'H(\jM^in cu ro ^ in 'HpjM<ru> ^ (m 



- 20 . 000 


- 20 . 000 


- 10 .000 


- 10.000 


0 . 000 


0 . 000 


10 . 000 


10.000 


20 . 000 


20 . 000 


DATA SET 


1 


2 


- 20 .000 


- 20 . 000 


- 10 .000 


- 10 .000 


0 . 000 


0 . 000 


10 .000 


10 . 000 


20.000 


20 .000 



DATA SET 


1 


2 


• 20 . 000 


- 10 . 000 


- 20 . 000 


- 10.000 


- 20 . 000 


-10 . 000 


- 20.000 


- 10 .000 


- 20 . 000 


- 10 .000 


DATA SET 


1 


2 


- 20.000 


-10 . 000 


- 20 , 000 


- 10.000 


- 20 . 000 


- 10 . 000 


- 20.000 

- 20,000 


- 10 .000 


- 10.000 


DATA SET 


i 


2 


- 20 . 000 


-10 . 000 


- 20 . 000 


- 10 .000 


- 20 . 000 


- 10.000 


- 20 . 000 


- 10 .000 


- 20.000 


- 10 . 000 


DATA SET 


1 


2 


- 20 . 000 


-10 . 000 


- 20 .000 


- 10 . 000 


- 20 . 000 


-10 . 000 


- 20 .000 


-10 . 000 


- 20.000 


- 10.000 


DATA SET 


1 


2 


- 20 . 000 


- 10 . 000 


- 20.000 


- 10 .000 


- 20 , 000 


- 10 . 000 


- 20.000 


- 10 .000 


- 20 . 000 


- 10 . 000 



DATA SET 




1 


2 


.400 


.700 


. 4 S 0 


.750 


.500 


.300 


.450 


.750 


.400 


.700 


DATA SET 




1 


2 


.463 


.763 


,513 


.913 


.563 


.963 


.513 


.813 



000 


- 20 . 000 


- 20 . 000 


000 


-10 . 000 


- 10.000 


000 


0 . 000 


0 . 000 


000 


10.000 


10.000 


000 


20 . 000 


20 . 000 



3 


4 


5 


000 


- 20. 000 


- 20 .000 


000 


- 10 .000 


- 10 . 000 


000 


0 . 000 


0 . 000 


000 


10 .000 


10 .000 


000 


20.000 

• 


20 . 000 



3 


4 


5 


000 


10 . 000 


20 . 000 


000 


10 .000 


20 . 000 


000 


10 . 000 


20 .000 


000 


10.000 


20 . 000 
20 . 000 


000 


10. 000 



3 


4 


5 


000 


10.000 


20 . 000 


000 


10 . 000 


20 .000 


0 00 


10 . 000 


20 . 000 


000 


10.000 


20 .000 
20.000 


000 


10. 000 



3 


4 


5 


000 


10.000 


20 . 000 


000 


10 . 000 


20.000 


000 


10. 000 


20 . 000 


000 


10.000 


20 .000 


000 


10.000 


20.000 



3 

000 


4 

10 . 000 


20 . ool 


000 


10 . 000 


20.000 


000 


10 . 000 


20 . 000 


000 


10 . 000 


20 .000 


000 


10.000 


20 .000 



3 


4 


5 


000 

000 


10 . 000 


20.000 


10.000 


20.000 


000 


10 . 000 


20 . 000 


000 


10.000 


20 .000 


000 


10 . 000 


20 . 000 



3 


4 


5 


900 


.700 


.400 


850 


.750 


.450 


900 


.800 


.500 


950 


.750 


.450 


900 


.700 


.400 



3 


4 


5 


863 


,763 


.463 


913 


.813 


.513 


963 


.863 


.563 


913 


.813 


.513 



-20 

-10 

0 

10 

20 

-20 

-10 

0 

10 

20 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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5 


.463 


.763 


.863 


.763 


.463 


3. 


DATA SET 












1 


2 


3 


4 


5 


1 


.650 


.950 


1.050 


.950 


.650 


2 


.700 


1 .000 


1 .100 


1 . 000 


.700 


3 


.750 


1 . 050 


1 . 150 


1 . 050 


.750 


4 


.700 


1 .000 


1.100 


1 . 000 


.700 


5 


.650 


,950 


1 . 050 


.950 


.600 


4. 


DATA SET 












1 


2 


3 


4 


5 


1 


,963 


1.263 


1.363 


1.263 


.963 


2 


1,013 


1.313 


1.413 


1.313 


1.013 


3 


1 . 063 


1.363 


1.463 


1.363 


1.063 


4 


1.013 


1.313 


1.413 


1.313 


1.013 


5 


.963 


1,263 


1.363 


1.263 


.963 


S. 


DATA SET 












1 


2 


3 


4 


5 


1 


1.400 


1.700 


1.800 


1 .700 


1.400 


2 


1,450 


1.750 


* 1.850 


1,750 


1.450 


3 


1.500 


1.800 


1.900 


1.800 


1.500 


4 


1.450 


1.750 


1.850 


1.750 


1.450 


5 


1.400 


1,700 


1.800 


1,700 


1.400 



Coefficients C0EF(1,J,K) 



I » 1 



i 

J 



h 

3 



I » 2 



J 

J 

J 



1 

2 

3 



I » 3 

J a 1 
J * 2 
J a 3 



K a 1 
. 4iag-^04 
, 134E-^01 
. 121E+02 

K a 1 



K =» 1 
.458E-^05 
, 143E+02 
. 133E+03 
K a 1 



K a i 

.994E+05 

.238E+02 

,289E-»‘03 

K a 1 



2 

.965E+03 

.390E-01 

.284E+01 



2 

. 106E+0S 
, 417E-»-00 
,312E^02 
2 



.230E-^05 

.696E+00 

.678E-f02 

2 



3 

.691E^02 

.282E-02 

.203E-^00 

3 



3 

.757E-»-03 

.300E-01 

.222E+01 

3 



3 

. 165E-»-04 
.499E-01 
. 484E-^01 
3 
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5. CONCLUSIONS AND RECOMMENDATIONS 



Compared to Che probe calibration approach described in Ref 3. 
the method presented here provides a much simpler way to both calculate 
the probe calibration surfaces from measured data and to apply Che 
calibration to on-line data reduction. The iteration required in the 
method of Ref 3 is completely omitted. It is noted however, that the 
Gauss Jordan Elimination lacks some sophistication. The system matrix 
for example undergoes the elimination without prior conditioning. 
Particularly in the case of the 3-D and the 4-D approximations, numerical 
round off errors have an influence on the accurracy of the coefficients. 
Consequently it is recommended, that Che Gauss Jordan Elimination 
routine be revised to use double precision constants and a routine that 
conditions the system matrix be added. Time constraints made these 
steps impossible for the author. 

Following these refinements, a calibration procedure, similar Co 
McGuire’s (ref 3) should be formalised. Since the 4-D- approximation 
requires large arrays and therefore extensive C? - memory, extended memory 
access (EMA) is necessary if the HP 21MX computer is used to perform 
the calculations. 
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APPENDIX A: 


Some 


Useful Matrix 


Conventions and < 


Operations 






Al. 


Submatrix 


Notation 














0. 


































... 


A* 

fra 




% 










0 ^ 


< 


oL - 




















Ki 






CLf 



















defined through elements . or submatrices A.. 

ij 13 

A2, Diagonal Lines and Diagonal Bands 1» order 



Diagonal lines 1. order and Main diagonal line 1. order 
Diagonal bands 1. order and Main diagonal band 1. order 




V 






A 




* Diagonal bands 1. order 
w,r,t. submatrices 



X 

% V 



% * 

A3 . Diagonal Lines and Diagonal Bands 2. order ^< 9 ^ 

Diagonal lines 2. order and Main diagonal line 2. order 

Diagonal bands 2. order and Main diagonal band 2. order 
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* Diagonal bands 2, order 
w.r.t. submatrices 
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APPENDIX B: Software Description: Flow Charts 



The following are given: 



Flow Chart 


Page 


MAT2 


93 


S2 


97 


MAT3 


98 


MAT31 


106 


S3 


109 


MAT4 


111 


MAT41 


119 


MAT42 


124 


S4 


127 
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F/ouv c^or6 /^/^7~2 



(entry ) 






r/0 
LI • Lo< 
L0= ^ 


f(otr>r»ievi'f^ • 

3LU (Z) 





93 



o 




^ , / 





96 



F/aw chori nSoZ 




97 



/* /o vv/ ^ari MArj 







98 



JJO 0^ A': U, / 





o 3 

Cxp<} 



? 

yj 

^u^ma^Ttx / chss 

C SUSM/ ) C*Sit^j 

SUS’^Ol/rfNB. MAVS> 



i 

D- 0 
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:x) 0 *- ^ 









1 


'Df+t 


\ 

< 


— — 



Suhry<or^iJC 
SUSfH cv den \T(/ 

fro/rt ^ocjthfouth kdt 

Mats/ hnh cv>" 
rtTc/ ploaP 
Ad,:/)‘-^Ci3M/cp,7/) 




HO 



NO 



100 




101 






£io X'j^A^r,ji:TOP, ! 







102 




I 



S^ri 3> /ore^^ 
hah) hoiryci ^idc 
'i^chr S’ 

I-Q 




I.L,I \ 










iPow/d = u~ / 




l.^.l ] 




L ..._.^ 


1 


^WRZ ~ Ml- 1 

r.-i-^i 





oe 



(^fyyou^ <S^U<9 Uon 

3. S’. (A'Siir^ 

rUNCT7QN 
anal ^a/ue 

ihfo ^ri) 




0 



1 



/^hi equciiioin 
^ysi£yn ctfler (T'7. 
s/irr)ityCKt>o 





L 





COQfPiCiey^^ 


\ 




Creturnj ^ 
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^ow chor'i 'H-) 7-J/ 





106 



/ 



Pf^oc^ed io ncxi, 

TnCf^^CL^ po^er . 
IPowRZ / 




I 

Copy dinned 
e'-SmQr^it dta^ 
QOnsdy' 

)^/x pf.l 



L 



V 

Copy 

SOSMIKC* /, J-/ J--€oSHirT.d 
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108 



Flow chart 'S3 




109 



o 

L_ 

jr- Yda.KZl 




1 





\ 






( ggruRH ') 



no 



i-hw chart h^AT~4 




/ 



111 




.112 






3o o7 




114 




117 




IS 

_i_ 



/Ae 



V 

(RLTUPN ') 
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Flow chaA MAT4I 



(entry 






'^i.ori io 
^dpe sech'om 
SuLrr\crfl-iCSf . 
Z. c/ass 
IP0 Wj? 2- 0 
Ml -2 *M-f 




02 

c>o <3 

4 

Sa^tna/rf ’^ 2. dasT 

^SUSM2 ) 

^eRQUr!N£ HaT4Z 
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00 asi'is77iKr,jsroi*,t U 



W -0 



00 oz u--Jstart,-^op,i 




V 

7 /-- 7 /*/ 



V 

12 . cfoM SU^MZ 
OS 

Sohrouhnf HAr 42 $ 
irth correct ph>ce 
iin suhmaf'r'iie 

/. class ^IaSHI . 
SuiW ft. 3h ^SMZ(a;si) 




Yez 



120 



3 

\y 




CO 0( f 



Co 06 UI-Z.M. / 



h 

u 



s 

V 




121 




1 



Define siari ond 
j^op demem-h 
for 

ISTOPm T!*N 
JSrWTt (Jt-I)MA/*/ 



Do os i’2sTMr,2snP, / 



Do CS D-^'*Rr,DSTVP, i 



os- 



r 

r 

i 




I 




123 



r/^n/ cl^ctri MA^^Z 





? 

Cofnf»(Ai* <y«a ^on 
4 . CtSin^ 
Fu/^cr/ON Sf- 
arid wri^f 
mio SueMS^r,7j) 



/ 

V 



124 




125 




126 



ric\A/ <dnari ^4 




127 




128 



3 

\/ 






1 











(return ) 
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APPENDIX C: Software Description: Listings 



The following are given: 



Listing Pag 

Two-Dimensional Approximation 132 

BLOCK DATA AB2 132 

BLOCK DATA DTA2 133 

PROGRAM DEM02 134 

SUBROUTINE MAT2 137 

REAL FUNCTION S2 140 

INTEGER FTJNCTION IEL22 141 

SUBROUTINE IEL2 142 

REAL FUNCTION F2 145 

Three-Dimensional Approximation 146 

BLOCK DATA AB3 146 

BLOCK DATA DATA3 147 

PROGRAM DEM03 148 

SUBROUTINE MAT3 151 

SUBROUTINE MAT31 155 

REAL FUNCTION S3 156 

INTEGER FUNCTION IEL33 157 

SUBROUTINE IEL3 158 

REAL FUNCTION F3 161 

Four-Dimensional Approximation 162 

BLOCK DATA AB4 162 

BLOCK DATA DATA4 163 
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APPENDIX C (Cont'd) 



Four-Dimensional Approximation (Cont'd) 



PROGRAM DEM04 


164 


SUBROUTINE MAT4 


168 


SUBROUTINE MAT41 


172 


SUBROUTINE MAT42 


173 


REAL FUNCTION S4 


174 


INTEGER FUNCTION IEL44 


175 


SUBROUTINE IEL4 


176 


REAL FLT^CTION F4 


179 
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BLOCK OAfA AS2 






FTH 



cw/ 12/ 



I 



BLOCK OAfA aB 2 
cannON / AB2 / 
A&AL AC 7 # 71 «B( 7 ) 
LNO 



&«3 






SLOCX OAfi 






3COCX OAfA 3TA2 
COBROn / OTki / ti,t 
A€At XLC2t6i*T(2d&> 
(NO 



FfN <••«*««! 
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I 

3 

10 

13 

20 

23 

30 

33 

40 

41 

30 

33 

40 

63 

70 

73 

90 

83 



3k10Z 74/175 iPT^^ «Tn 4,o»-.-#« 

AAOctAH dEf^oz <iN^uT»aur^ur* rA*€x* r4pEo#rAPE4«(N»un 



THIS IS A OEnQNSTSAnON ^tOOftAH AMO IT :hJwS Tni C3AAECT dSc 
Qf The r^t. 6IMAAV Lis«Aar using EXAnpie:* 



AUfHOt • 
OAtEi 



nans :E3N€i 

OkcenaU 24» i9to 



CoimaN / 0 TA 2 /xi«r 

8EA(. U(236>*Y(236} 

seal C0EP(7> 
lacICAI. IPIQT 



PT /J.i4l392634/ 

iHor /.rfloE./ 



lOl 



103 

m 

ioS 

601 



data 

OAlA 

FQRNATS 0EPO2 START 

PQRNAT inniti N6LL3* TNIS IS P80GPAN *>En02 4N0 d^*lL ScE# HQk T 
•0 USE the NAiVEtLJUS/ 70N TPL dlNAftf LiaPAPf. rOU8 InTEAcaT klLL 
•9E GREATtr AMAROeO if cASIeR/ZSN AHO FASTER PR0GRAPnr<6«///l 
FQRHAf UZH plotter ruffALlIsO.) 
format fZ4N OATA POINTS SENERaTEO.) 

FORMA T U4X#3HU1 IN«0X* $HX IN AX# 7X«AHrNlM» 7 X# 4HTNA X / 6X* 7NPL OTTER # A ( 
•IX«F 10.3)/4X#4NUSER#A< IX#F10«3) ) 

FORMAT (47H PLOTTER ANQ JSER ARfs OEPI'ICO ANO AXES ORAUN* I 
format iZOH OaTA points 39AUN.) 

FORNaT (34H enter IROER JP <*OLTHONIAL (HaRORl) * IPRINT anO ISTaT) 
FORMAT <2ZH CURVi FIT CONPLETEO*) 
format UdH CURVc PIT ORAwH. > 



format UMl/lNTl 
FORMATS OEMQZ STOP 



1/0 ASSlGAiMENTS. 



Ll*LOCLUf ISESSN) 
LP*13 

:n\i !hS: iiii 



INITIALIZE THE PVOTTfR* 



II* <IpLOf » *^ALL*lNlf5* (LP>* 
kRlTE Ul# IC^I 



YARIaTIQN op POITNONIAL DEGREE. LOOP STARTS NeRE. 



mi^trT.i 
M iSTQP*6 
00 06 lOROi 



Rl«MISTRT#niST3P*l 



GENERATE OATA SCATTERED AROUNO A PQtTNQnlAL. 



iiniv 

XSTOP 

hpnTsi 
ox 



•i;r 

♦2. 



<XST0P-XSTA«T1/(NPNTSI-1) 
NgPORl * 4 

caEP<I>* -i. 
mm\: ::3 oo 34 

C0EPf4)« .i 
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90 

93 

100 

103 

no 

113 

120 

123 

130 

133 

140 

143 

130 

133 

160 

165 

170 






r4/i/5 a?r-2 
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FORNATS IEL2 STOF 



L0*6 

Rwifl 

REARRANGE VECTOR 3 IN ARRAY V# IF N IS LCSS THAN 7. 

IF (N.£0,7) GO TO 02 
50 OL 

01 V(t#N)*V( I#6) 

FAINT SYSTEM MATRIX A ANO RIGHT HAND StCC VECTOR a 3EF0RE THE 
ELININATIQN4 IF IFRINT IS SET .TRUE. . 

02 IF CIFRINT) VRIT? CLO# 601) 

IF IlFRINT) «RIT: (LO# 602) (J»J*l»H»L) 

00 03 I«I»N#I 

03 IP nPRlHT) yRTTS tLO 
IF UFRIHT) VRITE (LJ 



603) I#(va#J)#J*l»M#l) 
602) (4*0*L#H#L) 



? 8 s: 



linear EaUATIQN SYSTEM V(N#N )«G(N )«V( N#HaU 

liMli 

04 A(t» J)«V( 1 ^ 0 ) 

OQ 09 

05 3ir)*.0 

CHECK FOR stability. 



/ 42 /. 
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suiftourtnc 









FT< <*»o*«*i 



oo / 1 i/ 



^5 



160 



16S 



170 



06 

tl 



54*jS 



r-i#«a 



100 


09 




10 


lOS 


1 

c 


no 




ns 


s 

c 


120 


11 




12 


12S 


13 


130 


16 


13S 


IS 


160 


16 


14S 




150 


17 

il 


15S 


1 



DO 07 

0(1 l*0( (U J)*«Z 

0(ii •$oftr(G(l 1 1 

30 10 {•i#A*A 

30 09 X3«I*N,I 
SK«$)(«y(i»i3i«y( j#ui 
f •$!(/ (0(1 i«o(J n 
If (Aas(F).or.As$(Fn 1 ri*f 
co«r iNue 

FI-A8S(FI| 

FI»ACGS(FII 
F|*(Fl«iaO.)/FX 

ST*«T XTE«ATXO.i TO FINO AiO XnF»aV£ lOLUTIOH* 

00 I9 

FIVOTISATXON ASO OAJSSlAf^ STEF F0« COLUW J. 

XX»J 

E-y ( J#4i 
00 ii 

IF (Aas(v( X* Ji }.ir.Aa$(EM oo ro ii 
OO TO 12 
Xa« X 

E«y(X#ji 

{ OMTIyOE 

F (lA.Io.JI O') 73 IS 

y(lx* (}*y (4«xi 
y(j» Xi*0( XI 
0(r<l*.0 

J ♦ I 
00 15 

M*y(X»Jiyy(<i«y I 
00 Is 

y^^AA^jy ( x#AA i-y (j#AAi«y 

XF^(V<X#l l•<a••0) 00 TO 32 

soiurioN. 

00 16 X«X*M*l 
0 (X 1 ••o 

0(MI -ycHa r<i/v( y«y I 

DO 19 lX*lpA»l 

I»y-xi 

IF ( v(X#X I .feO,,Cl 00 TO 32 

lyiWxU.i 

JM-JJ 

IF (J.eO.Xl 00 TO la 
iili •oT I i- 0 (j I #y( C»j I 

CO>*TlMUfc 
0(X 1 -odl /y( Ip X) 

COaSECT S0tUTX3N 3f AJ13U9T OF t?»0B* rfHlCX WAS ItyOLVEO I.y T^is 

MSaArrCN STEF. take r^XS wEClOa as start vector FCR Trie iEAT 
ItgRArioh. 

20 2?if2aiH2;i. 

COflFUTE THE aErtAllIlC ERR3R. 

OQ 22 I«1jN,1 

y(XpA i*«o 

00 21 J*IpMp1 

21 yTx#^i«y( w^i^OTJ iwAapji 
22 y(i#/ii«AU#iii-y(Xp<ii 
a*. a 

00 23 l*l#M.l 

23 R*R«>y( t#i5 (•yap/n 
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SUSRJUTIng (EL^ 



74/A75 P10*1<» 



f TH « *» 



tw/ u/ 






173 

180 

185 

150 

155 

200 

205 

210 

215 

220 

225 



•>SQ«T(M) 

(a.cr.^^SLutu ;u r:3 2d 
If CLA.&r • iTA'iui ca T3 33 
du zs 
aa»*6< I) 

!»«INT iTSTEW rtAmx 4 4M0 atSHf -»4 nO StOC ^gCTlR 3 IfTSR Tuf 

f ivoriSArioN* 



LX 



25 

26 

27 

20 

25 

30 



IF UMINT) M«:TE (LO/ 6041 

lo 

If tiavlNTi bRtri CL<1» 6031 I / I V ( I « J 1 # J«l/ H» I ) » fi ( H 
if (IFkInM «aITC IIU# 6021 U#0*I#ct#13 



00 26 
OO 2 6 
VU#0)«A( 

STOf LOOF TO fINO SOLUTION* 

courtNUf 

ffCARilANOi SOLuTtHH IN THf ARRAt 
CONTINUC 

00 29 I*lfN,l 
Vll»d>«0( X) 



fUlNT SOLUTlOMf If IfRlNT IS SET «TfUE. • 

(f (iffiNn bftirf CL’)* 6o>» 4 
If (ifaiNri (to* o06i c i*cc 1 1» 

If clfAlNF) wAITe (LO* 607) PI 
If UfRINT) MaiTc (LO* 60d) R*LX 

OfUCAiN OCR HAfaiX 4 JNO 3&S SECTORS 3 MACH ELtFINAflOH. 

IF IlftlNT) wans (L0» 6IJI 
If (IfAlNM WAITq tLJ* 602 ) 

31 If <}faiNT)*iaiT| (L0» 603) - . 

|fr(|faxHT) naiTE clj# 602) 

EtaOl LKOUNTEaEO* ST3F. 

32 XTta — I 

«fire (L0» 6C9) I*tWlX*X) 

GU)«GCI)/Vfl»i) 

aEruRN 

NO CONVEROiNCE Cf TH€ SOL'jriaN. 






33 



ITER«-ITER 

airuRN 

ENO 



V 
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\ 









FTpI 



et / 



fUMcrnu f2 



r*/i7> 



1 

J 

10 

13 

20 

23 



»€*t fUNCTIOH fZ (CJgP. 



ccfl»uT6 PUHcrnN v*tjf r*<soci» gp x :?-• »oltmq»<ul usi'*;; 

THe piQfiN£P a«Oe« r3 lfiC?.£Alu SPttO 4N0 4CCU-i*Cir. 



4UTHQA I 
04Ul 

\r4»i4aL£ 

tmuv 

N^aoai 

P2 



^kn% :E3N£a 
)£C£^8i:a 24» 1980 



exPLA»iArioN 



i«?NaAUC*!5ef?!cHNfl 

ti-y4LUfe4 #H£ag f ^(LL ac cn.ipurco. 

laoe* :p th« »gltno*^ui. :h the 



<l v)iacCTlQ»< 
jftr 



oiN4re Jf fH6 20««»attNor^:4c. 



Trp£ 

••••••• 

>)£4L 
*e >L 
iMTcOea 

a£4t 



31 



R£4L Ca£P<7) 

T-c2£?(HCQ£pI) 

DO Cl II«l#NCSOAl#I 

l*HCOZ^X»\\ 

r«rlxi*-C3£r <n 

F2*r 

9£TUBN 

£N0 



f 



145 



S(.ocic OAfA Aai 



?A/175 






I 



HOCK data A83 

ConnoN / A 81 / 

i€AL A(A4»A9l»a(A 

fMO 



A«a 
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StXK 


OATi OATAi /A/175 JPT•^ ^ ''it*' 


1 


SLOCK OATA OAfAj 

COHHON / OArAJ / Xl#X2*T 

«€AL U(lA#ie)*A2Udf iol* r(la«16l 

£NO 



PHUCRAH 9t10i 






FTn 






1 

3 

10 

13 

20 

23 

SO 

33 



«3 



fRQORAH 3Eno3 ( iMi*ur»our^ur» r ar€i# rAPcc*r AFEe# TJi^fc')«(N**ur I 



THIS 1$ A aEnjHSTRAriC^ ^RQSSAn AMO IT SHOWS fHe CJRrfcCr USE 
OP THE IPU 3XH*Rt Wl3AA«r USI'tfi £XA«Pl.6S. 



AUrHOR I 

QATEI 



»AMS 253N|R 
OECeHB^A U# L930 



CDIHOM / AflO 
CQAHON / OATAJ 



/ PtOFR 
/ U»X2#T 



101 



If! 



RCAL PLQTR(2d6) 

REAk Xl(l««lo}»X2U6*16)#r(16*lA) 

real coif (A91* SU5Ml(7»n 
InrEGER l»iU(:i 
LOGICAL {PLOT 

OAT A IC«T /lO-«» / 

DATA {PLOT /.TRUE./ 

OATa pi /4.1A159263A/ 

forhats 0EnO3 start 

FORMAT (///72M ^CLL^• T-4t 

• U USE The fiAPVELLOOS/ 73H r _ . . , 

• SE SREATLT AwAROeJ 4T cAStER/EjH AMO FASTER PP wG» AMM£'4U«7 / / 1 
FORMAT (22M PLOfTeR IN IT I AL I 2|0 . ) 
format (31H OAIA POINTS REAO ?R3M TAPS 3.» 

FOMMaT a<»«»5H<l.'irN»'>«**-<XlMA<»6A»5HX2*i:‘!»<»I*;HX2*A«»7<»*HrM£N»7X» 

• A«rMAl/6A» 7 nPL1TTcR* 1( 1 WFl0.3l/9I,4,MiJSi?,o(H»FlC. 3) » 

^^3TT5R asO_JS£R_*P|S OErlrjCO ANO AXe3 334#^.! 

XI- AMO x2-3{aewTIGH» IP 



15 PROGRAM PEM03 AMO #e^LL i?*» -fCW T 
--part* rOuR INTEREST WitL 



\l\ 

Hi 



format C30H NEASUREa DaTa SJRFaCc OPAHM.I 

>3^ poLTNaMiAc :ri xi- 

OSUR^Ice’^oiAWrt.J 



I 2 LH surface FtT 
C29H APPR ixZ.iaTE!: 



• RIMT AMO ISTATI 
FORMAT ‘ 

Format _ 
format C1H1/1HT> 

format (Al#<»Aia*A3l 

format C2I1S} 
format (SFIS.oI 
FORMATS 0EM03 STOP 





c 


I I/O ASSIGMMeMTS* 




* 


30 




* 




• 




LI«L0GLIJ( tksriT 
L0*6 










L0«6 

LP*13 






53 




•RITE (LI# lOlJ 
wRiTC CLQ# 601) 






60 


c 








1 

1 


1 IHXTlALaE the plotter* 
« 




• 


63 




!? ciPLifJ !all Up5 

WRITE (LX* loll 






70 


c 










I TARIaTIOM of POITMOMIAL 3EGREE. LOOP STARTS nER£. 

• 






73 










to 




N2STRT-1 

M2ST0P-6 

00 09 mOROR1»M1STRT*m 1STOP*1 
00 09 MOaOR2«M2STRT*M2STOP»l 










7^/175 lt>T «2 9 H 0 rti» 






100 

lOS 

110 

115 

120 

125 

130 

135 

140 

1^5 

150 

155 

100 

105 



01 



• 3EA0 1£aS0«£0 OATi. TiP6 ^ 4HP $ IrtULTiNEOUSC T Se*^C-» T-tt 

• fllNIflOrt ANO lAXMiJi yALJS QP c4CH t', (> 4Nf) f (^cJJHEO Tj 

• G€T THfe PLaT 4 iMf 4 iHc iou'<Oi«ie: V X r-*e:i op papc-^i. 

* 6 SxnO*LO * * 

• »9999<». 

Xl«*X • -99999, 

I2«IH • #99999, 

I2?1AX • -99999, 

TfllN • »99999. 

T>^AX • -99999, 
lViL»lOM 

9U0 (LO* dOl) iyAL«IRA4 
IF ( IVAL.eO.XCiTl iT T1 CL 

9EC000 (30#a02#UAK» r4P«i rU #.1 MTS2 
00 Oi li-UNPursui 
00 03 {2«i«^t^*<rsl#i 
IVAL«lOH 



IVAL«lOH 

02 A£AO (LO# Sul) t^AL#PA 9 
IF « IVAL. 60 .IC ID GJ n 02 

OECOOi (<»0 *oC3#I9amI 1 1 ( a • 1 2 ) • X2 ( I It 1 2 ). f 1 1 1 * 1 2 ) 
IF (xli h# 12 >.LT.Aifii*n ( 1 1 # ) 



IF (X 2 (ii!u) 

I? ■■ 



;#i 2 i .Li.Ti 

Ii#l 2 ).GT, XI 



I^n 

LT.XZTU 
■ TIM » 
’lAXJ 



. - _T 

A 2 '^li»i 2 ( i;# :zi 

T'ltK •T(!i,ri:) 

ii-iAx«xiai#i 2 ) 



03 



I *2iu< U>*‘5I**2’1AX) t2lAX-<2( II, !2) 

IF (t7u#I 2T .ar.T'iAx ) MAX .TiiaUt 



CONUHui 

vAire (Li# 103 ) 



OftFtNt FLOTTE) ANO JSE4 AftEAS: Or.AW AXES. 



xo 




TO 


3.0 


XiPMXN 


.0 


ZIPMAX 


1 C. a 


X 2 PMIH 


•c 


X 2 PMAX 


10. 0 


TPMM 


•9 


TPMAX 


ic.2 


XIUMIN 

xIjmax 


■hS 


X 2 UMIN 


-.1 


X 2 UMAX 


.0 


TUMIH 


- 35.0 


TUN AX 


3J.0 


ALFAXI 


30.0 


ALP AX 2 


15U.0 


ALPAT 


9C.Q 



FL 0 T«( 6 U 

pLgr«(o 2 ) 

FL 0 r 4 ( 63 ) 

FtOT* (ftA) 
PL 0 Tai 65 » 
FL 0 T«(O 6 > 
FLOra ( 57 ) 

?9fSUS! 

Pt. 0 T«( 70 ) 
f»L 0 TH( 71 ) 
XIL 






^aiTE (LI# 
• MAX# X2UMIM 
IP.( IFLOn 
•1#7#a) 

IF (IFLOn 

*IF^a^LOT) 

• 3 # 7 #d) 

• RirE (U# 



( ALFAX1«*I ) /ISO. 

( alFax2*fI>/1»0. 

(ALPAT ••I)/H0. 

13 

TO 

(XiF'lAX J / (XIUMAX-<1UMINI 

( x!i>'UN*XlUMAX-<lF1AX*XlU.1in) / ( XiUlAX-tlJIlNl 

( X 2 PMl,S«X 2 ’JMAt-C 2 P' 1 AI*X 2 UMrH| / ( X 2 UflA<-» 2 U^Nl 
( TP'^AX - TP*»l'< )/( T*J-#AX-TLPM » 

( ''»ilN*ruJiAx - tpiax*tu«IK j/I TUrtAX-ru^M i 

(XlP'iAX -<IP*I9 \ 

|X2P'<AX -«2P'«IH I 

( rPMAX - fPMN ) 

109 1 XIP^IN# IlPMAX#X2PMlNtX2PnAX# f PM IN# T PMAX # XI U«I N . XlU 
#x 2 UMax# rg'iM, rjiAx 

CALL AXci {X3#r3#XlL#ALFAXl#2HXl,^2#XlUMlN#XlJMAX»9f4F7. 
CALL AXES (X0#ra#X2L#ALFAX2#2’iX2#«2#X2UMlN#X2UMAX.>.fiF7, 
CALL AXES (X0#T0#TL #ALFAT # CHf #-2#rU«lN #TUNAX #%MP7. 
IQS) 



170 



ORAri MEASURinENT 3 a7A SURFACE. 






149 



OEiua 



jpr*2 






179 



lao 



o5 0<» II«WnPhT51#1 
00 Ok u»i*NPNr$2#i 

IF (iFtgn CiLL rH«rw t iptor* fFtor#<ur;# i2)» <2( n# I2»* r( ii#i2 ) J 

IF (Iplot ,Ano, I2,ei.i) CALL ^Lor (7H^^#rHaT*n 

04 IF (iPLor .ANO. li.ir.i) :aul puot cxpchT/ rPiJT»2 j 

09 IF llPCOr .ANO. li.ir.i) :all PCOT (XPLOT#TPtOT»2» 

imitE ttl# 106) 



tl9 



140 



149 



200 



209 



00 ^PxrE (lia ;c7) 

RfAO ( 9. •) '«aRO*l»MO*0»2AlPaif<TAlSTAT 

h3i*0«- 



NOi<dfll*f*QXOPI 
N3«0«2 

r ‘ ‘ 



STAfZ 



L»4Q*0P1»1 
/i«Naao«2tI 
IP (l.lt.i 

IF (rt.c r.i . 

IF a*«.LT,l .JP. (.•*<. uT.<»9» ,, 

CALL "AT3 (HPNTSlANP4rS2AL».1*CaeF»IPPIMTArSTATASU04l,I6iiPO6) 
IF (|EPPOP«^r. 0) oa TO 09 



.Qp. 

• Q». 



L,CT.7) 

7) 



S3 73 
S3 T3 



76 
36 

S3 73 06 



210 



219 



229 



03 07 II* l»NPNrsiAl 
00 w7 |2«I#HPNT|2»i 

7APP«Fi(CQEFAXlTli, U)«ir2( U» 12 ) .HOPnitt^MOPOAl ) 

or IF (IPL07 .AND. U.37*i) ;alL PL3T (XPLUTa 7PL07a2 ) 
00 ud j2*lANPs7S2# I 
00 C6 a* lA«PN7Sj.#i 

IF (IPL07 .ANO* I1#£0*1» call Pl 07 { X PLQ^# r plO 7# i » 
00 IF (1PL07 .ANO, a. 37a) CALL »t37 (XPL JT # 7PL0T# 2 ) 
•«IT6 (lIa ;C9l 
04 IF (IPLOf) call sgPAOe 



230 



239 



IF (IPLOT) CALL STOFG 



2A0 



STOP 0077 
ENO 



150 






♦“I* 



ci ! lit t ;♦*.<♦• 



I 

5 

10 

13 

20 

23 

30 

33 



Si 



33 



63 



70 



T3 



10 



33 



SJ8«0UTtN6 rtAT3 ( N^MTS iTS 7# L# •« »Ca£F # ! P« IN T» I S T* T* 5 I , is- kU-* I 



«»«6S£T SYSreH iAT^IX X iitj ^ICHf H4MD SlOe /feCTUa 1 -3^ A 
3U*»*QcyN0«I Ai»a«0XIf1Ari3N. T^6 uINCAP sQUAriQN ifSIei frt6N 
IS S0LV£0 »r USING A Gauss J3ftOAT4 cCl'^lNAT IGN. 



AuTMQAi 

0AT6t 



-♦AMS z«iN5a 

u» i960 



vxaxASLS 2 i^v*matiom 



NPNTSl 

Hanrsl 

CQ£f 169) 



TTfg 



IPRZNT 



rSTAT 



sua«uM#M) 

I£8«0« 



U(16» 16) 
S7(16* Id) 
T( 16«16I 
A(^9»*9) 
3(69) 



Nun<4€« JF OATA POinrs at c instant <i 
Nunaea GP gata joints at cimstant tl 
oULTnINIaL ;Ge?FICI£NTS 
Nunaga OF coePFicigNTs for th6 polt- 

NO-llAt IN Tl GiarCnaM, 

NUNaea of cieFFicuNTs pop tn£ polt- 

NONIAL IN tz GIRtiCTiaM, 

COMTajL OUANTirr GF -»9IMT GUT 
•2 palMT euuAfiaN stst^n G-EPoas 
ANO AFTEa The Gauss jnpiGAn 
SlIMlNAfiaN. 

•l >aiMf gaUATION STSTf.i AFTCP 
THE G4USS ijaOAN tL;.’*:r;4TiON 
>0 PaiNT iASIC EGUAflOflS -- 
ACTUAL PA^AiETiaS. 

<0 NO lUrPUT AT ALL. 

coNraGL STATisTsiwS ra it 

/lSTAT/*2 PaiNT GaUlNATC: 

tJ3 5l5!S?[i5fi5' 

3NLT paiMT caana:. 

S(?r A3S0LJTE 
itT acLATive 2?r.ijP, 

1. CLASS IN rne :rsT6P 



stat/.i 






P9 integ 
GF oat A 



iNTsGea 
I M Ts Gc a 
if Ac 
(MTcGca 

iMTEGea 

iNTcGga 



liTcGSa 



/rsi 
rsTAT>2 
lsTAr<o 
suaNATa IX 
NAfaiX. 

saa3a flag aETuaNco m the calling 
PaGGMAN. 

••t NO jaaoa getecteo. 0 no joa, 

•-1 AM ?aa3a GCCuaeO :’i t-»u GAUiiiA« 
ELiNiNAruNf :5EFFic::.MTS ahe not 
AC^ utATS. aEvi&M pafai;t!ps« 

Xl-COOaOlNATSS 3F- The ?A*’i POINTS 
Xi-CJOaOlNATCS GF The ’'aTa pGIMTS 
r-c33aoi'urs5 3 ^ the oata »gimts 

STST 5 N MATRIX, 
aunf HAMO ilOf YeCTOa. 



aEAL 

INTrGrt 



Hi Ac 

--4L 
He AL 
*?AL 
HcAC 



IBPOITANTI TH« OATA POINTS XI# *2 Y haVE TO 96 »a3YlO€3 

TwSquGH Th| LAaELeO CUN.H»1P 3L3CA GaTaJ A.nG ThE 

StsTen HATaix A ANO thc ’ic-tr iano side - = :t3h a 

AIF AlTAILAdLS THROUGH ThF 'ABtLtO CO^NON iLlZ'i 
SAVES A lot of 7haT PaeCiOUS Thing 
cIllEO iEMOftr. 



COMMON / OATAi / Xl#X2#r 
COMMON / A33 / A# ) 

«UL Xi(l6#l6)#X2(l6f :6)*r(l6#l6) 
aEAL A| ^9#49|# 3U9) 

REAL C0EF(69)#SU8MlIN#M)#STATa6»16) 



»I3#nIH STEPS NO SOLUTION VITHIN THE aCCURaCT 

if./* 



formats mat3 start 

loi FORMAT (///7H after., 

• OF. elQ. J# UN { I GUir./) 

601 format {///51H CaOER GF A P»« Ol I MA TCON >nirNOMtAL IN U GIP6CTI3N. 

• 8X#I3/4IN OiOER JF APPROX INA TION ><5 lTMGNIAL IN U C I R ?C f I 3N. 5 X » I 5 
*//60H suMiEa :f gata points ^rrn comstjuit xi»i6x# i N uMata 3 

6O2*F0a*I} ?///62H^^ jTSf SN^MATRlX^A^4.N6‘yStrGR 3 3EF0RE GAUSS JOROAN E 

• linihatiom/i -n. 

603 format UAX.I3MX# 12) )) 

606 FQ«maT (2X#I2# I3( IX»lPE9.2)#U/< (6X*13 (:x«IPE9.2)) ) ) 

605 format (///6«N FGUAfrON SiSTEl AFTE® GAUSS JOROAN ELIMINATION/) 

606 FdaMAf itttZSH iJcFMcIENtS C3£r(I,J)/j 

607 format (12X»3H4 • . 12# lU a X > ( 2 ) # IX / 3( 7X # IZ t 3X# 12 ) / )) 

60a format (2x#3nI •#I2#l2UX«e9.3)»iX/3(6X#L2(U*c9.3) /)) 

609 format <///3RH OaaiNATES OF GATA POJN^S • AS GIVEN/) 

-* - * AS APPHOXINATFO/I 



610 FQRM4T 

6 11 format 



w / / J '•H Jmi ^ ''VJ * ^ ^ * g 

(///66 H jaGlNATcS GF^OaTa points - AS 

(///30rt AiSOLUTE saa3a aT ;AcH pqint/i 



151 



SUBRJUflNe ^kTl 



7^/175 JPT•^ PnO^-» 






«</:2/w<>. it. 



95 

100 

lOS 

no 

119 

120 

129 

130 
U9 
1^0 
149 
190 
199 
160 

169 

170 



ti: 



612 POPflAT (///35‘» ^tLATlVk Af cACrf »TI*tT IM t/J 

613 FOPNiT (///2*r» »r*risTics i j«*«4aTt// 7r»;iHArft«Av»t ^«»o9»i 

•2/7* . i9Mrt*xinu*t e»«n4 \T 

614 FOilnAr <///2Ir» STAfiSTICS Sun'i a« t t / / 7r , . irtAVLP. aGc c?- i3*.2PP9. 
•2»2M ?/7*»19Hn*x[.'iUi c»6a« AT t • • I 3# :Xi ; IJ • » i 3# a a» 2 c ?) 

619 Puawir APT£9»£J»AIH STEPS •<'7 iU'lTiaN Mr aCCU^ACT 

•0F» clO.3* llH i : QJU./I 
fOSiVATS 1A?3 3rOP 

^LQCLUC ISESSil 

IER6Q«*1 
ITANAX«60 
EPScaN«.lG0E^O7 

PRIMT APPflOZINATl3N PaRANETSRS* IP IPRIMT IS GREATER TKAH w 

HC*0RI*<.-1 

MUR0R2-W-I 

IP ( IPRIMT.IT.OJ G3 TO 01 

MRITE UO# oCli .S0R0Rl»NQR0R2/riPnTSI«MOMTS2 

START TO preset EDGE ScCTrQ^t SJAMATRIOcS 1. CLASS. 

91 IPQUR1«0 
Ll*2*L-l 
00 09 

IP IK.GT.L 1 GO T3 02 

TOR E06E SECTIJM SU3iAT«ICES I« CLASS. 



JSTART«<K-l)*Hn 



JSTA 

J5Tgp-<*« 
g5 Tq 03 

right HANO £06£ SECTIQH SUSHaTRICES 1. CLASS. 



02 ISTarT»<a-LI91*1 
ISTQP«< 

JSTAaT*(L-ll*H*l 
iSTOP«C«N 

PRESET EOGE SeCTIQM SUSMATRICES 

03 CALL RAT31 (l«PNT>l«NPNTS2»n» I PQwRi . SUS51) 



04 



1. CLASS IM STSTEH HATRU JSIHG 4aT3a. 



I»ISTAar#ISTOP»l 



Ih'l. _ 

00 04 J*3SrA6T*4STQP*i 

a7i*Ii« 



SUSoIIIIaJU 






09 IPQUR1*IP0«R1»L 

COPY OEPINEO SUSNATRICES 1. CLASS OlCOHALt^. 

00 07 Ji*2»L#I 

I5TART«<U-H*'’*l 

I3T0P-I1*h 

j!tart«ui-i)*h»i 

JSTQP«Jl*« 

oo oo i«inARr.isrop» L 
30 06 /•iSTA«T»jSTOPf i 

A( |«A< Ij,l) 

CONTINUE 

PRESET RIGHT NANO S106 '/ECTOR I. 

I-O 

90 06 Ll*l#L/l 
IPOWRl.Ll-l 

OO G4 -K 

|pQWR2«Hl*l 

08 9(n*S3(NPNTSl9NPHTS2»lPaRRl# IPQVR2«1) 

Nl«t*N 

PRINT SYSTEft NATRU ANO RIGHT HANO SlOE VECTaa S» IP IPRINT 
IS SET TO 2* 

IP IIPRINT,hE.2I go to 10 
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175 

110 

115 

190 

199 

200 
209 
2 X 0 
2X9 
220 

229 

230 
239 
290 
295 
250 
299 



MtlTe U 3 « {) 02 ) 

«AIT€ (L3« 'y*i3l U»4»i»Hl$U 
OQ 09 

09 wA{u UQ# [« (A( r* n 

w 9 (re ua« 6 on ( u 

I?t IHT»1 

CAUXS J090AN HI'lliUriQM. 

10 6 ^SLJH-€^SL 0 N *‘<1 

CALL I&L3 (NI# lfEtl« IT99AX«£^SLaN) 

If (irEfl*Lr*4^} oj ra 27 

FtiNT cOUATION STST69, IF [FtllT tX :6T T*n 1. 

IF i IF9lNT.Nt*I) 30 TO 12 
V917E (to# 90S) 

««1T€ UO# 9V3) (J*A*I#Nl«U 



IX w 9 



00 U 

ifc ?L0» 9v») 

«RIT€ (LJ# 6C3) 



M.%) !• (A(lW)»i*l#n«X)#d( n 



OEFIflk C06FFICt€flTS AnO (JAI'IT T^Ci, IF IPOINT IS SaEAfE^ f^AH 0* 

12 CCWTimJE 

30 U Li«l#Nl#I 

13 



tlpirSf I lIJo) 31 TO IS 



wairg (LO# 9 C 6 ) 

WAITE (L0« oon ( 12# C2*I»^#I) 
“ 19 Il*l#L# i 

TAT«<I 1 -U* 9 #l 



OQ 



ICS TAT _ _ 

ICSTOF* (1*9 

19 WAITS (lo# 9Ci) ri#(CQ£F(ic)#r:*tcsTRr# icsTiF#:) 

WAITE UO# 907) ( r2#IZ*l#l»U 

statistics# if ISTAT A.wrTmi3 else aUT o* 

IF (ISTaT.OT.O) IEAA3A-2 
xsrAr*iA8S(isrAT) 

COnFUTE 49SOLUT6 SAAOA* 

C •(AFFRQtlhAnO ikVJi • 3IVE9 VALUE) } 

* 00 16 I1»1#MF*ITU-, I 

16 ?TAf?Il5li|JFJ(cii^ * XICI* 12)#X2( Il#:2)#^iQP0Al#NaA0A2)-r( U«I2) 

CCilFUTF RELATIVE SA^QA# fF ISTAT WAS \,Z:Z THAf» li»0, I,c. IcAAQA*!* 
C •(AFMOxI/IArSO VALJE - UvEW VALUL}/(GIVE)4 VALUE) ) 

SS 

00 17 r2*l#9AHTS2,i 

17 srAf(U,i2J*STAf7ll#I2)/rl ru I2) 

‘ SF*-Ai!JuS*^lAli* A‘/€Aa;E EAAOR AnO FIp»0 lccatiom 

18 sAAAVE-.J 
FAAnAI*AaS(STArU»l) ) 

I9A*»1 

jnAX« I 

00 19 ri»l#NfNTSl#l 
00 19 I2«1#.WAN rS2# 1 

^F ( AsUsTAfC I1^;2)).LE.€AA9AX) 30 TO 19 

AlHAAaS “ 

I'^Al.t! 
jWA X* f 2 

19 eAAAVE«€AAAVE#A9S( STAT( :i# 12) ) 



ir i • 9 j j > • I I i 6 

gAANAX-STArdi^U) 

If*Al “ 

^WAX 

cAAa _ . _ ^ ^ 

E«aav«-e»a AvE/ (wAnrsi*N»wrS 2 ) 

If aiTAr.^E.2i 30 22 

faint 3A0IHATES If OATA FOINTS ANO AFFROIMaTEO VALUES# 
IF /ISTAT/ IS SaUAL 2. 

wa[?| IlS# 6 C 3 ) Cr 2 /I 2 *l#NFNTS 2 #l) 

00 20 r i«i#NFNrsi#i 

20 WAlTe (LO# 909 ) 1 1 > ( f ( 1 1 # 1 2 » • 1 2 • 1 # NPNTS 2 # 1 ) 

• AITS (LJ# 603 ) ( :S r 2 *l»w#NTS 2 »l ) 
wAlTE iLO# 910 ) 

WAITE (LO# 0W3) ( 12 # 12 *l#NFNTS 2 #l) 
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260 

26S 

270 

273 

210 

283 

200 

203 

300 



3U8ftOur(MS 1AT3 






21 

22 

23 

20 

23 

26 

27 

28 

20 



00 21 ii*i»NPHrsi«i 

waU£ (LO* 6C01 {l#(F3(C0EP>Xl(lI#X2)«XCtIl*I2)»HQt0Rl,N04a.<21*l4* 
•l»#«PNTJ2#U 

tiriirg (L3# 603) (I2»I2«l*^PNrS2*U 



rilNf SRaO«5* If /IST4T/ «&S I 0* 2. 

IF { iSMr.Ni.D so n 26 

IF (IERR0I.£9.U iO TO 20 
POINT AiSOLUri E080RS. 
wRlTf aO« 611) 

MRlrl (LO* 603) (I2*I2«l*'(PNrS2*l) 

00 li tl«l*HPNf3if i 

«aire (Co* 604) ri«(srAr((i#( 2 )*( 2 «L*Npfirs 2 *i) 

WRITE (LO* 603) ( l2*l2*l*'«PNr32*U 
00 TO 26 

PRINT RELATIVE ERRORS. 

wrHI IlS* ( I2*I2*1*WPNTS2*1) 

OO 23 I l*l*N«'HrSlil 

:iin !t3: tti\ 

PRINT ERROR STATISTICS. 

IF (1EPR0R.E0.2) WRITE (LO* 613) Eff*AVE#:RAX*gNAl«ERRAAZ 
Is VRIfc (L3* 614) E4RAV£*!'1AX*WN4l*e4A^4X 

FAfAL ERROR 0ETE:TS9. SUrPUT ERROR RESSAOC AN PRESET COcFFICIeNrS, 



ITEfi»-(Tc« 

IER«OR*-(ERROR 

!4|ipi:iiNwi 

Syfriti: IIW 

GO TO 12 
RETURN, 

RETURN 

ENO 
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1 

5 

10 

IS 



Z3 



30 



3S 






50 



SUM3UTI(ie lAf3i 7^/175 lPT«i PTi 

SUitOUTIhfc 'lATll S2>^* IPn^Pl/ JUftiiJ 



P«6S6T SU8n*T«IX l, CLASS SJB^l THl STSTCn ^*T9U 4 PQIt 4 

SO^acTNC/ltAL 4PP40U.i4n3N« 


AUT«409 I 
0AT6I 


^A^«S 

OfcC6Mdcit 31# 1990 




V4AUaL6 


EXPLAN4TIQN 


TYP6 


KP#iTsi 

NPNTS2 

N 

IPOVAl 

SUI*3X(M#rt) 


iU^96ft IF 3474 A0i?<fs AT C^MSTANT «l 

*iUfl964 IF 3*f4 PgiMTS AT giMSTANT <2 

'*U.Y9£R JF COfcFFIClcNTS tmC PQL/- 

MQ.1I4L IH t> OIRSCTtOM. 

powfip Fa* T-»g a V4LU6i IM THE 

CU'iiAriQl. 

SUBYATRU 1. CLASS IN T^g CYrTgN 


i'ir=‘.ep 

iNftGe^ 

LYTCvtC* 

iNTeSe* 

46al 



01 



a&Al. SUBMl(rt»/i} 

!P0WP2*0 

OU 03 

If (A*iT.H» SO TO 01 
TOP 6006 SfeCnON gLElCNTS. 

I-l 

GO TO 02 

tXGHT HkHO HOGS S6CTI3N 616^6NTS« 



PP6S6T 60C6 SEGTlOff cL6n€iTS I"* SUBHATPIX I* ClAiS JSIMG S3. 

02 SUdMUI«<i »*S3( '^P^TSl»riP*«TS2*tP0yAL,IP3tf((2*O) 

03 IPQMK2* IPQwP2^; 

COPY 06PIN60 £L£«6NTS OUlOtlAUf. 

00 QA 
00 cA 

agfuPN 

6N0 
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1 

) 

10 

15 

20 

25 

10 

35 

«0 

45 

5C 

53 

• 



74/1?; jpr-i ?iofi4 FTfi 

RIAL FUNCTION n ( NPN T S U 4 »»«T S2 » I FQWR 1 , I? 3WP2# I T I 






PIRFORN SUN^AflQNS FJ »ReSkT FH6 ELTMCNTi 4(1*4) M F-ti; 
STSriN 4ATRU 4 A40 r-»t lltHeNTi 9(!| :N :H*; ^Ur»f -*141 

Sloe vICrCR i for r-tl lIicAR C4U4TMM SYST^i, 1^,41 iuST 
9fc SOLVIO t-UR 4 ^aLFNjNUL 30-APPP.ajanA T ION, 



AUTHQAl 
OAFE I 



U.'<5 :k3NER 
OECEiSER U* l^eo 



vARiAiifc Explanation 



NPNFSl 

NPNFS2 

IPORRI 

IP0NR2 



OAFA POINTS /*T CINSTAHT XL 
OATA POINTS at Constant <2 
FOR TnE XI xaluE: in The 



ir 



S3 



rlUN8£« Oi 
MUN86R 3 
»Q.XeR «0l 

:unnation. 

P04tR for TNC t2 /ALUCr :M TH6 
sunnation, 

CONTROL •aranETER 

•0 r values are not irCLUOEO Ih TnE 



•I 



XU16* le) 
X2U0#16) 
ri id»iN) 

INPQRTANTi 



SUNNATION, 

(-» STSJEN 4*|Rrx JLC^ENTp 
T VALUfei ARE INCLUOcO In THE 
SUNNATION, 

(-> RIGHT NANO :iOC /CCTOR ELEN,) 
RESULT OF SJNNaTION. 

Xl-COOffOlNAT£S OF the 04TA *0INTS 
<2-C30ROlNATtS GF The AaFa POINTS 
r-CJOiO MATES OF THE DATA POINTS 



^N F ^ (i e R 
iNFe^eP 

integer 

INTEGER 



INTiGsR 



REAL 

R-aw 

AL 



The data points <1# *c A.‘;o T have fq a p«jvioto 

THROUGH the LAdELcO COnnON 3L0CK GATA3, 






03 

34 



05 



CONNQN / QaTaS / <l»X2 , t 

REAL xm6#l6)«X2(l3#ld)>m4Ald) 

s*o. 

00 C5 N1«1#npnTSI#1 
00 03 <2«1#NPnTsI#L 
IF (IPQHR 4 .E 0 .O) GO TO n 

51- 1. 

IF < lP0«R2.eO«0) GO TO 03 

52- X2 (AlfAA)** (P0^R2 
GO TO 04 

S2*U 

ST-T|Kl4<2) 

IF (ir.E0«0) ST«1« 

S*S*Si«k*ST 

S3«S 

RETURN 

ENO 
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1 

S 

10 

15 

20 

25 

30 

35 

40 



55 



60 



65 



74/175 

XMTESEA FUNCTION tgl33 tHl 



FTN 



GAUSS JGROAN gLlNlNATTON diTH NaTFIX A ANO /ECTOR 9« TH£ 
/ARUSte Ut3 IS SET TO I, IF T-E ELINInATION IS COwPUsTfeO* 
Uti LESS than 0 INOICATES an ERROR. 

AUTHOR I 
OATi< 



VARIABLE 



IEL3 



A(49«44) 
SI 44 I 



'^ANS ZBSHSR 
0ECEN6ER il* i960 

ElFLANAriON 

wsirsrnuUioHs 

caHTaoL FARANsrea# returneq to the 









^ -AH cLININATinH VAS 

FcRFdRHEO correctly. 

<0 ERROR EHCQUHTERfcO OURIMG GAUSS 
Iqaoah ELiniNAriQN. 

STS TEH HA rail A 

RIGHT HAHO SIOE VECTOR 3 



TTR6 

!;f£Gi5 

IHTcGcR 



»5AL 

real 



IHRQRTANTl ’•ATRIX A ANO VECTOR S *UVE TO 3E RROVIOEO THROUGH 
The LASELtO CaHHOH SLQCA a33. 90Th a kHO i .ILL 
RILL UHOERGO CHANGES OURING EXECUTION# I.E. ThE 
EtlfilATiaH. 



CONNON / A83 / A#6 

REAL AI49#49j»3f44T 

ilH3«IEL3 

cRS-l.GCe-<3C 
00 G6 l«i#N#l 

?? ?K.E«i!l'io TO 02 
CONST«-A|K.l)/A(I,II 
SO 01 

a7a#r i^Ai a#ji*cqnST«a(i#a) 

IF U.EO.n A|A,JJ-.0 
01 CONTINUE 

3UI«9IA}*CaNST«9m 

03 A?l5a)iA<njt/CQNST 
All# I )«t. 

8(t)«si iwcoNsr 

00 06 J«l#N#l 
lERRQKO 

00 04 

IF ( ASStAl J«<) I.GT.ERS) S3 T3 04 
lERROR* Ic9RaR4l 
COHTIftUE 

IF (A9S <9Un.GT.EPS) GO TO 35 
UaAQa«lERRGR*l 
- le I ItRRJR.Gl. H*U GO TO 07 
06 CONTINUE 
. S3 TO 08 

” ihv.hL, 

08 RETURN 
ENO 



06 



05 
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ec/42/»o. 



I 

5 

1C 

13 

20 

23 

30 

33 

VO 

V3 

30 

33 

60 

63 

70 

73 

50 

53 



J?T«2 



M 'I 



SUtMQUTlhi liLi (.'4«lUfl#irSrt4l*MSLa.1) 



04USSMM 6Ll!’IH4nj.'4 ^IT‘< Ll^gAfl J3U4‘»’nM STSTe* ^ (Me ^IC-iT 
H4NI> SlQh /£Cr38 ^ IS l'4T6a?4T60 WII-H fHt ifSTcV <4- 

T4U 4 iH n£ iraix V). iica is lie muh^ck ir irei4ria>iS -t- 

TuaHfo fo Tie CkLLiHQ aa^ca4M. me :T§?4riQN rj :n(>-iv6i r^e 

SQLuriOM sa ^kHt Tties. j'«riL th£ ae^Atniii taaoa is withcn 

TMt «4C^*lTU06 Jp :?SL3'«. S-iJUtO ThC MU««c« a? ireaAfirjNS zX- 
?iSQ ifflH4x# rse suaaojTiife rca^uurcr rne C4LCuL4fiJ.4 4.'<o 
KlS iteA M 4 ^cGaYIVc tf4UU^« 



4UTHQRSI 

04 re I 

vittme 



ireii 



ITaTt4X 

ePSLON 



41 49#43I 
3(491 



IHPOtfAHTl 1AT«IX 
Th6 



X0NR4Q yo6eL£«* 44NS :c9nex 
oecenaea 3i* i9so 



eXPLANATlQM 

MuA4ea*:?*rjuATfaNS 

Nu.iaea JF iTeaATiQNS rn as peTuasea 

*0 G4JSSI4H SLIfimATIOM PCa.fORrttO 



fYFS 



_ 'JiZlkH _ 
CQaaeCTUT. 

<0 HO cawycaieiji 



nf th; sjLjTi'iN 

44xf«(ji r 



I"i£ 06SI9 
treaAil'jM 
irayAx 



iTepArroH sic? 



oeiUeo aesQtuno^ i 



ST3T|.3 14T«IX 4 

aiGHf ^4N0 Sloe yecTQR s 



c?sioN. 

AFTea 

an.*' 1 

£M4rJj»t STfePS 

F saturiQN 



FOM 



iNTt^ea 

mTiCga 



iNTSOeA 

«cAt 



a^AL 

^e4t 



4 iHo yecTo* 3 HAVC to se p-toyioea TnaouGrt 

tA 3 €l.eQ COrtHON aloe’s 433. SOT*^ a h'*0 ^ -ItL 

- . 



FL^HlAri^N^^ :^4^G£S QURIMG cXCCUlIQM* l.c 



COnnOM / 463 / / 

8c 4L y(4q»6Ol«G(30)#4(49#50l*5(5QI 

logTcai ip&Int 
QAT 4 IPRIfit /.FALSe./ 

OATA FI /3.l4139i304/ 



FO««ATS iei3 START 

601 Fg»H4r (///39H i4Tar< 4 jho vexTo* 3 voR =li*um4T10n/i 

602 FOariAT ( ( 4X* 13 (8 ( • 1 2 M I 

601 Fg»«AT (2XW2# L3( U. LFf9.2l/(C4t#13llX* :FS6,2I ») I 

604^FaanAr (///23-t »ivor iSiea res srsfert itcfatiomsscm 

603 F0HH4T (///2X# I2»24-» <aeFFin EiTSi 4US GAUSS/I 

606 POFMAr H2x«I2#iX#I<>E<).22) 

607 F^nAT I///39-* SLcI <ST?a Ze I Le Hy c <TaR? C I TTwl N<e L ••eX2.3l 
605 FCF«4r (27 x, 12-» acSIOIUt •#6l>,3#5f4 </.Cl#n»L2*< 1 f FF aT 1 

609 FQRrtAr <///2CX«U4 -e-^tca ^ • f2 • 1 h, [ 2 , J-i | 

610 Fnfl«*T i/i/ziH irsTs^ 4 a:h eti5iNAi:jM/i 
FORAAii m3 srsF 

L0*6 

X«M*l 

aeARRANce vecroR 9 ih array y* if h is less ihan %9« 



IF (N»e0.49) G1 10 32 
00 QI I«l#N»t 
01 y(i*^i*v( i«3C) 

FtlMT SYSTtR MATRIX A AHO aiGMY H4MQ SIOC VeCTOR 9 3eF0Re Twe 

eu.ninATiON# If iftntj IS SET .raue. . 



ao* 

(LJ« 



601) 

602) 



02 IF (IFRINT) yRire 
IF (IFRINT) «R!Te 
OO Cl t*WN#I 

03 IF (iFRiNT) «RlTe (L3# 6031 I # (V ( I • J I » J«l# ) 
IF (IFRINT) yRITc (LO* 602) 

Lt^eAR COUATI09 SYsreft v(M,ij*«on.Y(«,rt*u 
00 04 I«l«M,; 

00 04 ««U5*I 

04 Aii#4 )«y( i«4) 

90 99 

03 i(I)-.0 

CH6CX FOR STA3IUTY, 
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90 

95 

100 

105 

110 

115 

120 

125 

110 

113 

190 

195 

130 

135 

160 

163 

170 



suaRjurifie icci 



/9/i/S 2P1*i 



F • a*«» -*3 



0^ &6 

06 6I1I-.U 
00 06 
00 C7 

3 ! 

00 10 
JJ» l*x 

00 U9 13«1#N,1 

09 S/l«SK»vaf (3I*V(J*13) 

F«$K/(GU )*GU II 

[f (A8S(F I.GT.AaSIfXn 

10 dONflr«UE 
F1«a3S(FI 1 

START ITERAriQ:! T3 FX^<0 i^O X.IRtOVE rOLUTION* 

00 27 LX«l#999»l 

IUr«IX 

90 13 

RIVOTXSATXCN ANO GAUSSXAM STcR FOR COLUnh J. 
ri«j 

E«V< J«J I 
00 12 

IF (ASi tV< X«4I I.GT.A^StEll 33 T3 11 
GO TO 12 

12 coi«rlNU6 

IF IXX.cJ.J) GO T3 19 
00 il 

G4X l-yl XI* XI 
y<xx> xi-y<j«X) 

13 V(J« X )•«( X) 

S(>4I«.0 

19 

00 13 X««'**»H*1 
R«V( X« J )7V ( J* J ) 

OO IS 

VrX*KKI«vI X«i4K l-V( J«KKI *6 
13 cofirxwui 

X«XI*ca**OI GO TO 32 



IX 



16 






SOLUTION. 

G(Ni*y(n»Ai/vlN»xi 
0 19 IX-l#K,l 
•N-IX 

IF (V <X«X )•£Q••9l 33 TO 12 
G(X )•¥< I*N| 

00 17 4J«I#fC»l 

IF (J.EQ.X) G3 TO 19 
G(X l*Gt X)-G(4) «V| X«J I 
ccNTiHue 

Gai«G( ti/v<x#xi 



T SOLUTION 9T itOUNT OF ?Rt0t» y»JTCN yAj I>iy«LV60 iN THfi 
:a< 6 r-»XS VECTOR AS START VECTOR ’Zn T^E lEXf 



. 20 

1. 

C 



UERiTiON’STEF 
XTERArXON, 

00 20 X«X«N*1 
G(X I'SI XI*G(XI 

CONFUTE THE RENAXNING ERROR. 

00 22 I-l#H*l 
V(X*H|«.0 

00 ll 



n V<I»HI«yiI#9|«Gf JiAAf X#J) 

22 v<X«N)« A( I.n|oy(I«H| 

R«.0 

OO 23 X*1#H#1 

23 R«R«V ( X#H I«VC X#N) 
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su6«auri'«€ tcLJ 






175 



ISO 



US 



1^0 



195 



200 



205 



210 



213 



220 



225 



29 

25 

26 
27 

20 

29 

30 

31 
12 

11 



f)«SQAr(R) 

CF (i)«t.r*£FstaiM sj ri 2s 

IF ( tx,ar* rj 3) 
3a 29 i*I*f«*i 



FtiMT SYST|« 

FivariSAf iON. 



H«r«(X k ANO iUCHT <tANO SIOC YECTIt S AFTe9 rH£ 



F (iPR|Nri ^A||| 



t ipaInh I 



(L0« 6091 LI 

(1.3* 6021 U*J*L*N»II 






OO 23 

IF (iFt*!**?! -«IT« (LO# 6Q3I 

IF ciFtiNTi •air§ aa* 602 i 

§8 n 

Va»Jl*A( 1#4) 

STQF L0OF TO FIHO SOCUTtOd, 

COMTiHUi 

agAOAAMCE SOLUTIGM li THf kHHkt 

CCHTIMU6 
00 29 

v(X*ioi«<< n 

Foxnr SQiurraM* if i»9int is set .tfue. • 

IP (IP9INTI WBIT6 (to# 6051 i 
IF (ifuimti »«it§ (l::, b06» ( 

IF IXFAINTI mAITc Iti* 6073 FI 
XF IXFtXNTi 49XTe ItO* 6081 ^»LX 

OKUCXEN OEfl nAffllX A 0*10 3ES VECTOAS 1 ^UCH ELXR XNATI3r«. 
IF IXFtiNTi VAXTe U3# 6101 



F (XFAXmTI vAIT: UO* 6021 

a ii ‘ 



_ _ __ i«i/f«*i 

n ishist'i :tm itii 

•etuan 

EAAQt EnCOUNTEREI}. STQF. 

AETUAN 

HQ CONVEAGEMCE OF- T'lE S3CUTI3M.. 

ITEA--XTEI 

aetuan 

ENO 
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fUNCTUM 



1 

3 

10 



20 



25 



30 



35 



45 



Si 



03 



/4/175 


lPT-2 ^lOrtP FTn <« 


C 


REAL FUKCTIQN 


f3 IC06F» <l» 42>N0R0Rl#fUPCP.2) 




COnFUTE FUNCTION 1/4LJE r*FUl»X2) 3f \ 30*<»0 l TN*j«l *L JS£«G 
Tn« nONNEk SCIcNE in order O InCP.c/.CT SPEED iHO ACOUriACT. 


iUTHOB 1 
OaTEi 


NANS ISiiNeR 
3eCE*15cR Hi 1990 




VAilASLE 


eXPLANAriON 


TYPE 


u 

nSaoii 

noroRE 

f3 


WfWUL*!^fUI{I5Af? 

Xl-V4LJt# anERE f AILL 3C ^ONPUTtO, 
*2-V4LUt* aNERE f AlLL it CQrtPUTSO. 
lajER Jf Ths POtfNONUL l:i Tnc 
XI DIRECTION. 

DR3ER Jf TnE pOLTNONUl I*i T-t€ 

*2 DIRECTION. 

OROlNAft OF INE 30*P04.YN*VHAL* 


• • * • . • . 

AEAL 
RE4L 
Re 4t 
iNfsv^ER 

iNrSGER 

RE4L 



CW/U/C'#. !•». 



86*L Clgr t4Qj,0<4<J» 

WC3EP l-HORO«l*l 
NCQgF2«MaAQA2«X 

CO^FUTf COEPFICIEMTS 0 3*Se3 QH COEFFICICHTS COEF *N0 <2. 

30 02 1«1*mC 0£FI#1 
IC- I*»*CCgf 2 
r«C3€P4 IC» 

30 Ul ii«l#K0fl3fi2#l 

J-NCa£F2-Jl 

lC*U-l)*NCGeF2*J 

T-CO£f C IC )♦<2•r 

COftFore ir 34S£3 C3cFFi:i£MT$ 3 kHO XI 

r« 0 (HCQ£f 1) 

00 Ui JI«l*NO03Al#l 

J»NCOtf l-Jl 

T-OU )»U*T 

F3«t 

06TU4M 

ENO 
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8C0C^ 2kU 



I 







1/5 J^»r•^ 


SLOCK 


OkT 


K AB^ 


COnnON / 


AB<» / A»d 


HcAL 


kits 


» 0^ 1 • 9( 94 1 


eNo 







Ft'i 



»*ic 



( 
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SIOCA DATA OATAA 7A/175 JPf-2 FTm aOii/i'J. >•;?.;> <*4CL l 



1 ItOCX oata oataa 

C3nnON / OaTaa / <l>iZ#x3#r 

SEAL XU9i 7»3»»A2H#;«;)»ii(9«7*3)*r<9# 7/3) 

ENO 
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^aoaRAfi 0 £.i.u 74/179 



rjn 4«4t4</9 4^*v^* i ''44C 



I 



9 



10 



19 



20 



? 



^•QGiAII Ol«04 UMfUT.OOTPUr# 74^61# TAP€6#’’*?t7» r***fc?*l‘^*»ur I 



THts ts A oEnoNsraAMON aaca^A.i amo :? sriJus rne c^asEcr <jse 

Of TMfc TPL 3MAnr LlsaAar using sxArtPtfis. 



AUTHQRI 

OATEl 



HA**s iEiNfa 
3 ECeMlE«i n# 1983 



CMNQK / AFLO / FtOTR 

COiinoM / oaTaa / ii#A2«x3*r 

alit X^? 97 ^fitix 2 f 9 » 7 * 9 )*l 3 ( 9 # 7 » 5 )#ri 9 * 7 # 9 ) 

tut CQCF164WSU9N:149*49I*S(JBN2(7«7I 
INTEGSR taAw( 6 l 
IGC1CA4 IPiOT 



29 



30 



39 



M> 



A 9 



90 



99 



OATA icht /:on« 
OATA IFLOT /.rauE./ 



OAlA FI 



/ 3 . 141592694 / 



XOl 



FORNATS 3 ENQ 4 STAar 
FOFNAT C/// 77 H r-tis 



•0 USE THc WAAVELLOUS/ 7UH fFt 3MAaT LlSFAFf* 
• 9E JaEAfLY awaAQcD if sA5IEa729H * 



FQAftAr 

FQanAT 



IS »«OGFAN f^ENHA ANO -E»tL SEE# T 

ar LiSF.AFY* TOua iNTc-tbSr auu 

AHO ^ASTt/t FACGaAfi^lM(i«<//i 



(2ZH FlOTTcA iMITULUEa.) 

<31M 3»rA »3tsTi U*D Fao.i fAFC ?•! 

FOaWAT (19X»9MU 1lH«4f, 9NXl4iJ|,oX, iMX2N(M#dX# :4X2NAX# 7X*«HTnlN« 7X# 



• 4 rtT«AX/qX# 7 flFt 3 Tri)#»(.X#FL 3 , J}/ 9 X* 4 MU;CF« 6 (;X#F u.m 
9 FOitNAr 7 a 7 h i>L 3 TrEA 4 .i 0 JScA AAtS JErt »:£0 AMO AXtS Oaa«M.} 
b format < 31 m mSaSjaeO 3 aTA Suacicsj o^Ayn.i 

7 format (74m E-4TjA OaOEa JF fovtnUMIAL :*1 xl-#. X2- 4N0 U-Ol9SCrX0 
•H, laacNT 4N0 is7aTi 
a format (29H SuaFACE FITS :OMFt6TEO.J 

H Oh xffaOxImaTEO SUaFACeS 3F.AVN* ) 

Ml/IMTI 
|A1#6A13| 



FOaMAt 
FoaMAT 
format 
format 

FOaMAf . 

FOftMArS OEM 04 



|A 1 # 6 A 10 



TQF 



I/O aSSICMMEMTS* 



l!*L 0 GLuI xsinro* 

tp 

iaiU Ui 4 ICII 
vaiTE (L 3 * 9 G 1 ) 



XMXTXAtXZe Tm€ flotter. 



IF «IFVCTI CALL 

waxTC Tti* 1C2) 



IMI7G (LFl 



70 



79 



90 



99 



• VAaiATXOH OF FOCTHQ.MXAL SEGSEE. LQQF STARTS 4EtE* 

MiJTafSr 

Mj|TOF«b 



?8 IS ? 8 ;SSK 5 i!?l?:?HI 3 J:t 

30 IQ *iQapa 3 «M 35 raT£ i 3 ST 0 F#i 
iTM«<Moao 5 i*ii*<-»oR 5 a 2 *ii*<Mo*i 3 a 3 *iJ 
XF (XTM.GT.A 4 ) GO T 3 10 
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90 

90 

100 

lOS 

110 

119 

UO 

129 

130 

139 

190 

199 

190 

199 

160 

169 

170 



06^09 



79/175 JPT-e 



rTi «C/i2/iV. 



9EA0 «lASU«cO 04f4 pa'll T4pe7 kf*0 3 1 -lUL T ANt *lUSl T Sc 4<CH fHg 
AHO /ALJ6 OF 54C*< ANO f (aeJUlRcO TO 

»iT THE PtQT rfXfHLN Th§ aJUNOAaigr Qf k SHEET Of ?*Peaw 



*!hIN 
xinAX 
x2h 



01 



lOM 

«g9999« 
•99999 • 
*99999 • 
•99999* 
*99999* 
-99999* 



XZHAX 
T«IH 
TflAX 

OECOOfc <45,702WaA«l HfHTSl#HfHT52*HfHTS3 

88 31 

OQ <33 I2*i>HPHTS2>l 



02 aEAO (L0« 7C11 lVAt.»ta«M 
If « iVAt.EO.ICMTl G3 rf) 

OECOaE (oO>7C3jiaA^I XillU 

•2*131 

IF (Tui*u*i3i .LT.THtfi I 

If ixiiiuiz# i3»*iT.xiiAxi xi'^Ai*xu:i> 131 

If <x|tII*i2*m.>;T.x29Aii x4^tx.f2l!** !:*hi 

- if *iT*T9Ax I r9Ax *Tai*:2*m 

93 CONriNUE 

«R17E ai* 1031 



*i2*nitx2n; 



I3l*xi( tl*I2*I3l*r( 11,1 



OEfIMc PLOTTER A90 USER AREAS) OIUi/ AXCS. 



XQ 

TO 

xif mn 

IIP MAX 

X2P9IH 

X2PHAX 

TPHifl 

TfrtAX 

XIUHIN 

XigftAX 

X2UMIn 

xlunAX 

TUHIH 

TUflAX 

ALf AXl 

ALf AX2 

ALf AT 

fh8?5!311 

PLQTa(<»3) • 
PL0Ta(6AI • 
fLorai65i • 



• lA. 

• 3. 

• «0 

• IC *i, 

• .0 
• 10*0 
• *0 

• lO.O 

• XI^AX 

• X2*<19 

• I 2 HAX 

• r^iN 

• thax 

• 60*0 
•120. 0 
•UO.O 






CALFAT «*n/I39« 
XO 
TO 



37X19? 
ftara (6di 
PL3Ta(69t 

at 

TL 



5 XIP9AX 
i xl* 



X- a5;i!.uu,tM;aa 



uiix-x’u^lNt 

U.“1A X-X I 'J^I i » 



I*HIm*XIUHAX- 

(XZPHAX ) / 1 »’‘J'1AX-<2JH1N| 

lX2P1lH«*2U9AX-t2fHAX*X2U«:M) / 1 X 2 U« a X - X 2UH t Hi 

( rP«*4X • ro<nT>< j/( f«j.^Ax-r u*»fH i 

I t»-*[m«tuhax • fP'^iXPTUniM )/{ ruiAx-rgHlM ) 

cxipiAX } 

(X2f^AX -tSPHlH ) 

I rfiAX - rf9[H > .> 

WRITE Cir, ICA) XIP9IH, IlPH4f,X2f1IN,X2f'^AA»Tf'<IN, Tf9AX,XlUHlN,XlU 
•fX4X,X2UHIH,X2UHAX, TUHIS, TUHAX 

IF (Ifi.ari CALL AXES UO,Ta,XIL, ALfAXl,Z4Xl,*2,XlUniH,XlJH4X,4Hf 7* 
•3*7,0) 

If UfLOTJ CALL AXES (XO,rO,X 2L, ALf AX2, 2HX2,*2^ X2UnlN, X2UHAX, 4Hf 7* 
*|*^tljLOT) CALL AXES tXO,TO,TL ,AtfAY ,2MT ,-2^fJHIN ,TUHAX ,4HF7* 
•RITE (LX, 1091 
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175 

180 

185 

190 

195 

200 

205 

210 

215 

220 

225 

250 

235 

2^0 

245 

250 

255 



PflOGRAn OE304 



24/ 1 /i 1^T*2 ?ti<ynw 



F(N 



4i/4^/2'*a 






0R4M 4e4SUAgHcNT OAfi lUVFiCES. 



o5*C4*!J2I^«FNf!3ll - 

88 8: [i:i;KS3fh;t 

IF Tinon C4tc rn^rw (xFi.3T»rPtor*xir:i# Z2#i3i#x2( Ii#i2#i3)*r< n# 



• ! 2 #tan 
IF (iPLor .AMO. i2*c0.i) :all 
IF ( iPLQT »ANa.^:2,iT.i) ;all 
OO 05 l3«l#NPNfS3#I 
00 05 I2«I#NFNri2# 1 
30 05 Il-l»NFNrSl# 1. 

(IFlOT) call Th«TJ 



»LJT 

FLOr 



(XPtnr, rPL0T»3l 
(XFLaT#rpL0r#2l 



^IF (IFCaT) CALL ThSTJ (IFLaT#TPL0T#XIi:i#l2#I3).X2CXl.I2»I3)»rUX, 

!?5lot ,A«0. U.cO.l) CALL ?L3r (XFL0r,TPLaT#3l 
95 If (IFlor «AKO« li.iT.i) CALL ’COT (XFuaT » f OLOT • > ) 

MRITE (LI# 106) 



CALCULATE SURFACE FITS THRQua^ OATA POINTS* 



06 «!f2 (lI# 107) 
«EAO ( 4* •) 

NQaOfll*N3flOAl 
NOfl082*Na«OR2 
NOROfll*NQROR3 
IFRlNf-O 
lSTAT-1 

iirssssSKl 

N«NOflOR4»l 
NN«N*N 
IF <L«Lr*l 

I? 



•4a{ORl«HQROR2»NOaOI(3»IPFZMr#ISrAT 



Ki. 



:8|: 



_.0T. 7) 



IF (L*nN*lT*1 .UR* 



T3 36 

f 8 n 

1 6 ^ ) *0 









ir.o) 

1C8) 



50 ro 39 



* OftAW affroxinateo surfaces* 



oy5r!s:t;A«^!5;i 

00 07 Il*L>NFNrsi#i 
OO 07^12*1. NPNTS2# i, 
rAFF«FA<C0EF#XUll#I2 
FN0RQR3) 



I2«13)*K2(Zl#12»I3)«X3rZI«I2#I3)*NORORl>NOflOR2* 



(XFLQT#rFLaT*Xl< 1 1» 1 2# I 3 1 » 1 2 C U# 1 2# I 3 ) • T AF F ) 
eO*i) CALL ^LOr ( X PtUT, fPLOT* 3 J 
all plot (XPLOTj TFL0T»2I 



IF a/Ldr .ANO. I2 

30 09 I 2 * 1 #nPnTS 2#1 
OO 09 Ii-l#NPNTSl#l 

rAPP*FA(C 0 SF*Xl(Il* 12 « 13 )*( 2 ( I 1 #I 2 »Z 3 )#X 3 ( Z 1 » 12 #I 3 ) #HORORl#NOflOR 2 » 

CALL TMRTJ (XFLOr>fPLOT*Xlf ll,l 2 #r 3 )#X 2 (Il#Z 2 # I 3 )»rAPF) 
If (iPLor *ANo* ii.£ 3 *ii call plot (xPL 2 r#?PL 3 r* 3 ) 
iRlH^uF* *1091 (XFL 3 T#rFL 0 T» 2 ) 

99 IF UFLOn CALL iUFAC* .x 

10 co«rxfHi8 



I rERNlNATE GRAFZCS* 
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2*0 



0e^0% 7<»/l73 JPT«2 PTm ‘•.•♦-.ia 

[p (IPLOn CALL STQI»« 

STOP 0077 
SHO 



1 

3 

10 

15 

20 

25 

30 

33 

45 

30 

35 

60 

65 

70 

?5 

80 

•5 



f4/i7i lPT-2 



PIN 



£W/ i 



su8«0uriM€ ^^474 i'«p'irsA#^PHrs2#rtPnTS3#L/i#M#;a£PiiP«isr# isrAT» iu^rt 

•l#/iN#6U6n2/;£ARQiii 



P463ET STSrg5 '1AT9IX * *N0 ^ICMf H4Hn 5 106 V6CT0* i P*3< 4 
♦ 0-P0LThCWl4L 4i»P«‘3Xl.i4rnN, Tng €QU4rijH ifSTci TrttN 

ts SQLV60 ir JStHO 4 54USS iai04N £L:'nM4tlQN. 



4U THOAl 
04T6I 



H4MS Z69N6P 
a£C6'186il n$ 1980 



y4iuii6 6XPL4N4ria:< 



MPNTSl 



Hi»MrS2 

MPNrSS 

CQ6f<44» 

L 



IMINT 



HU486* :P 3474 PQI4TS 47 CnHST4M7 ll 

NUH5|» 2f : “ " “ “ 

su-ia?* OP 04 7 _ _ . 

pQtrsartUL :j 6 ffici 6 n 7 s 



2f 3474 PQIS7S 47 ^0NST4N7 tZ 

-kf4 ?oIs7S AT innS74M7 ll 



SU486« 



. OF C35FF7(;16 m 7S 
M01I4L ;m tl 3I«cCTI0M, 

•“Sea iF :36 ffi:ic 47S 



Fnn 7Me PQtr- 



, . PnP. TH6 FOLT- 

NOrtUL 72 3r«6C77a«. 

SU48e« OF C06FF lClc*«7S FQF, TH6 FOtf- 
fiait4L I.i <3 3IR6C710M. 

Ca.'47ft3L JU4y7lTr Q5 3u7 

•i PRINT ^QU47iaN STSfEi 3 lFUF. 6 
4N0 4F7ea 7Hg 3nP.04a 

|LI5IN4n3N, 

•I PrZnT 63U4710N STSTEN aFTCF 
7H6 G4USS J3»34n fL:ilM/,7!0N. 

>0 P8iN7 64SIC £3U47I0MS aiTH 
4C7U4L parameters. 

<0 NO 3UTFU7 47 *Lt. 

CONTROL OTATISTSICS TO 3E p*IMT 60 
/IS7A7/-2 PRINT OBOINAT^; OF ^474 
4N0 APPRQfr-4T50 POINTS 
4N0 3EVI47I0NS /ERP.aa). 



/iST4ry*i ON 

iST4r>o 



ONLT PRINT |P«a?' . 
i|7 49S0LU75 EPP/IP, 



SUINUNM*NN 



NN 

SU8N2(fl*N) 



I6M0R 



UI9*7»3» 
X2<9# 7*3) 

4(64*64) 

8(64) 



IST47<0 GET RELArUE sPPJR. 

) 5 u 3 .NA 7 RII U CLASi IN T §TST 6 fl 
M 47 R IX. 

3 IN 6 NSI 3 NIN 5 y 40 I 48 LE FOP SU 8 «l 
:UdN 4 f?lX 2 . CLASS IN SJdWATRIX 

1 . CLA>>. 

•RROR FL 4 G 467 UBNE 0 70 THc CALLING 

^ROQRAN. 

.♦L NO ERROR 06 T 6 CTE 0 . GOOO J 05 . 

•-V AN ERROR 3 CCURE 0 IM T^i« GAuSSlAN 
cLlNINATIONi C 06 FF!i:£NfS ARE mi 
ACCURATE. REVIEr PAPahETcRS. 

tl-C 3 JR 0 lMA 7 «S OF 7 N 6 OaTA POINTS 
- -- .. POINTS 



X 2 -C 00 R 0 lNAr£S OF 7 
<3»CuOagiNAr5 S cf- 
T-:ooR3t uTEi :f tn, 

irSTEi MAfRIt. 

RIGHT nano n06 VECTOR* 



-HE 3474 P9IN7S 

7Nc 3474 PllArs 



7TFc 



INfeGER 

[n7: ^cR 
Be AL 
IN7EGER 

INTEGER 

INTEGER 

Integer 



INTEJER 



INTEGER 
Re At 



iNTEGeR 



mOT 



BcAL 

real 

BcAL 

ws At 

Re At 
RE At 



inPORTANri TMf 0474 »OTNTS <1* »2/ IC *^0 T HAVE TO RE PRO- 
VI3eo Trough The CAiclEO COHNCN iLZi< 34744 AMO 
THE STSTEN matrix a anO TMO RIGhT hamQ 5I0£ VECTOR 
9 are 4V4tLA9LE TNRauC.f •Nt LASELEJ TUNnJh mLOCA 

A84* TmIS^SavES a tar of that precious thing 
called NENQRT* 



CQNNQN / 04T44 / XI*X2*X3*r 
CQ.NNON / 494 / 4*4 

REAL n(9>7»5)*l2(9*7*5)*X319*7*3)*r(9*7*9) 

real 4(64*64)*a<64) 

• CAL COCf 164)* SUaN101N»NN)*SU8N2(N».N}»ST4T(9*7*5) 



rOtNATS rata START 

XOl FQRNaT (///7M after* I3*41N STEPS NO SOLUTION VITHIN THE ACCURACY 
•QF*E10*3* IIH ) I 4UIT*/) 

601 FQRNAT 1///31H OROER 3 A P PROI I -A TION PCLTNONTAL IN t\ OIRECTION. 
•5x, I3/5IH ORDER OF aPpROx INA f ifiN )<5LrNQNUL IN <2 OlBeCTiQN.Son 



OPOER OF APPROX IMATION POLYNONlAL 



•/51H _ _ , . . . , . . _ . 

P/40H NUN9ER OF Data POINTS VlTN CONSTANT 
• 3474 POINTS mITh constant X2 *i 6X * ij740H 
•7N constant X3*X6X*I3/) 

602 FORNAT (///62H iTSTeN NATRIX 4 AnO VECTOR 3 9EFQRE GAUSS JORDAN E 



. , , <2 312 . 

IN X3 3IRECTI3N.9X* U 

XI, 14X,I3/<.3h mijm4cR of 

NUH9ER OF Data points *1 



168 



susKQurzNi 



74/17J 3P^•^ 



PfN 



9Q 

95 

10« 

lOS 

no 

lis 

120 

125 

130 

139 

140 
149 
190 
199 
UO 

169 

170 



• L IfllPIArZON^ I 

603 FQiMAr H 4X.13(8t» 12) n 

604 Fg»fiAT (2X,I2#13(lX#li»e9.4).IX/U«»x.;3(iy»IPe4.2) n ) 

609 FO«»n*r tftfsdri c3jaTijh jfsTen xPiCP dAHS? jopoak Um<ATl3M/» 

606 FO^nAT I///26H CJtFFlCUHTS C3g F ( I# J # X ) / ) 

607 FQRnAT t2X»3M( •* 12) 

605 FORKAr (X2X*3H< «i :2» I1(9X#12)» U/3(7 Xa:3(9X#12)/)I 

609 FGRKAr I2X.3HJ •• 1 2 » 1 2 ( 1 X* c9, 3 I • I X / 3 ( ♦ X * ; 2 ( IX* £ 9 . 3 ) / I ) 

610 FQR«Ar UJti^'4 OROInaFIs 3F DATA OOlNTS - A$ llaN) 

611 FORMAT (/dH T(l» J*« {2* M) I 

612 PDMMAI C///A6f4 JRQI'iArES DF ;aTA FQ'ST; - AS APPROU MA T 60) 

613 FCRMAf (///30H AdSdturl rPR3R AT cA^-i ^tlNT) 

614 format (/UH SRROR a*3»»(2»IH)) 

m format I///31H •^cLATWE jRRJf AT ?ACH PHTf^T «) 

format (///21H irArrSrrC* 3JMMARTi//7X*13HAV6RA0fe 6aA3»*26X«lP69. 
«2/7(* 19HMAXlntJM ERROR aT I • » 1 3* i X* 2MJ • » X 3* 2 X* 2 • » 1 3* ix , i R E 9. 2 ) 

617 format (///2lH STATISTTCS 3JMM4Rri//7X*.3HAV6?.AD6 6t4R j« * iaX* 2 - F9, 
• 2#2H 1/7X#19HMAX1MU'4 error aT 1 • • 1 3#2 X# 2 IJ* • i 3# 3A • 2 -(K«» 1 3 * U* 2 PP 'v 

•|*2H :} 

618 FORMAr (///7H after* l3*4iH STEPS HQ SOUJTION rillHlN THi aCCJRaCT 
•0F.E10«3*UH ; I DUIT./I 
formats maT2 stop 

tX«C00tUI ZSESSM) 

L0«6 

lERROR*! 

rTRMAX«60 

€P$LQM«.1COE-007 

PRINT APPROlinATION PARAMETERS* IF ZPRtMT IS 3REATER Than J« 

NO*ORl«L-l 

NQR0R2«m*I 

N0Rp«3«N-I 

IF TiPRImT.iT.O) so T3 Cl 

tfRXie CC9* 601) NariORl*MOROR2*MQROR3*NpnTSl*NPNTi2* MPNTS3 
START TO PRESET eOOE SECT13N SOftMATRICES l« CLASS. 

lP0riRl«0 
U«2«L-1 
MN«rt*N 
0 “ 



01 



?? U.VxUhHh ra o^ 



CLASS* 



TOP eOGI SECriOM SUdiATRICES !• CLaSS* 

istart-1 

lST0P«Mf« 

JSTART* IK-1)*MNa 1 
JSTQP«<*MN 
so TO 03 

right MAMO edge section SUSMATRICcS 1* 

02 ISTarT« <K-W l«MN*l 
XSTQP-( <4-1 1*1 )«HH 
JSTART*(L-ik)«Mr4Al 
3ST0 P«L»mn 

PRESET feOGE SECTION SUiMATRICeS 1. CLAJS IN STSTEM MATRIX USING haT41. 

03 CALL MAT4I (NPNrSl*NPNTS2*NPNTS3*M*M«ZPOMRl/SU8Mi*MN* SU3M2) 

I1«0 

DO 04 i«isTAar«isroPf 1 

ii-ii*i 

oJ'Sa j«jstart*jstop*i 

?9 



ltfRl«lROuRl*] 

COPY OSfINEO SuaMATRICSS 
Ll«L-l 
88 ?? 

TART« lll-l)*MN#l 
rop*ii«MN 
,_TART-(Jl-l)«MN»l 
JST0P*J1*MN 
DO 06 l«ISTART*ZSrOP»L 
DO 06 J-JSTARr*JSTD»* I 

5 6 Ail ».HN* a* ^ I 

7 CONTlNUfe 



I* CLASS OXGOMALLr. 



Ill 
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SU8K0UTM6 



7A/175 



FTN 



175 



(10 



IIS 



lio 



115 



20 C 



205 



210 



215 



22C 



225 



210 



235 



^les^r fttsHT ^am 5 sto£ vEcn^ s. 



r-o 

OQ os Ul«i#L«t 



08 



rwKii.ti-i 

go os 

OQ 08 

i(iUsA(NFNrsi*H^*irs 2 *'«^Nrs 3 »z^a 4 «i#rpaup. 2 * t^atfi 3 #u 

N 1 «L«NN 

P«(NT SrSTt/l naT«(X 4N0 S(G>4r '4kH0 SZOC VECrOR 3# IP IRitt^r 
IS SET TO 2. 



tP f (RKNT.NE.Z) sa TO (0 
VRtll (ca* 6021 

58^41 kuin' 



09 8RtT€ ILQ» tiOA) L#(«(r#J)#J-l#il#i)*eC) 
wRltE <C(j» 603) 

tn (NT»i 

OAUSS JORDAN ELININaTION. 

10 ERSlON«(RSCaN«Nl 

CALt IfUA (Hl,ITSR»IT«1AI#fe8$L3N) 

(f UTfcR.LT.C) OJ TO 32 

RKNT EOUATtQN SrSTE/1# IP I?R(NT tS SET TQ I. 

(I (iRtXNr.NE.;) SJ TO 12 
WRlt| (LO# 0C6) 
yflrl (LO# 6C3) (J«J«(#NI»U 
00 II (•L#Nl«i 

tl WRITE (LO* oOA) I#(A( t* J)»J«l#NI«U#8a} 

WRITE (CO* 603) 

OERIHE COEfFlClEHTS ANO RRINT THEN# IP :*RlHT IS ORfeATEii T»<an 

12 CONTINUE 

00 13 L1«UHI>1 

13 COEMaU-SUI} 






IP ( iparNT.Lr.-)) oa ro lo 

mRITE (LO# 60b) 

00 16 ll«l#c#i 

WR|T| UO# 11 

wRlti UO* 608} (l3» (3«Ui*l) 

I 2 -l)*H*l 

iN*r 2 «H 

WRITE (LO# 6UW) [2j(C3E^aC)#i:-ICSTTr«ICSTQR#l) 
■ Rite UO# 608 ) (l 3 #U*t«H»U 

statistics# if ISTaT ANTTHliO ELSE OUT o. 



0 IS 

5 HI 

ICSTOF«m»U* 



16 IF (ISTat.eq.O) so rn sa 
IF (ISTaT. u,0» IcRROR'I 
JF^( ISTAT. GT.OI LcRR0R«2 
iSTANlAiSlisr AT) 

CONFUTE AdSOLUTE ERROi;. 

C •(AFMOAlftAlcO VALUE • 

2 S n 



GIVEN VALUE) I 



17 STAT (U#Z2#i3)*s<»<ClEF#n( II# 
*Rl#NOROR2»f«ORORi)«V(ll# 12# (3) 



)«^A(COEF#n( II# 12# 13 )#X2( II# 12# (3 ) * X3( Ii#I2# 13) #n090 



2A5 



250 



255 



CONFUTE relative ERROR# IF ISTAT WAS LES; TNAN ZERO# I.e. lERRJR-I. 
C • ( AFFROX IRATEO value * GIVEN V ALUC )/( GIVEN VALUE) I 

IF (lERROR.NE.t) GO TQ (9 
00 18 ti-l#NFNTSI#L 

IS STAT7u#I2#I3)«STAr(Il#I2#I1)/T( Il#(0#:3) 

ERROR STATISTICS) CONFUTE AVERAGE ERROR AND PiNO LOCATION 
OF naiinun error. 
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2M 

265 

27C 

275 

280 

2SS 

290 

295 

300 

305 

310 

315 

320 

325 

320 

335 



SUBROUrUE 



74/175 IRT-2 



PTH 



2W/L2/2'*. 1^* 



19 6 ««AVE «.0 
gflRnAx«Aas(STA r(i»x«;n 
INAX«l 

KHA 1 «L 

DO 20 tl«i#NANTSl#l 
00 20 I 2 «L#HPNrS 2 #i 
DO 20 I 3 «UNPHTS 3#1 

If (A»S(3rAr(ii«:2»t3)).L£.E9RiAX) oo ro :o 
tflBnAX«STAT(U# 12 f 13 } 

I«AX*(I 
JHAX * 12 

XHAX * 13 

20 Ei(0AVt«£*4BAtfE*ABS( STATC U»l2»m) 
PAAAVE«EhaAVE/<N0MT31«H^NTS2«*4PNrs3) 

If (lSTAr.r«E. 2 ) ^0 ro 29 

PRINT OBOtNATES OP OATA POINTS ANO APPPOXINATtO /AtUfiS# 
If /ISTAT/ IS iOJAL 2 . 

00*22 t32l,Nl2}s3*l 

SiiH Il8» oud! i u* i2«;«'«»Nrs2*i} 

21 2«l?i ftoi'ScJ J^t0(r(tI,l2#r3}*t2-:#NRf|TS2>U 

22 WAITE <L 0 » f>C 3 l ( 12 * I 2 «i*NPNrS 2 #I) 



23 



BITE _ 

WAITE (L0» 612) 

00 24 13«1»NPUTS3* I 
WRITE (L0« 611) 13 
WBlTf (L3f 6C3) (I2/12- 
00 23 Il*UNfNTS;#r 

vAlfE (Co* 6 C 4 ) 



ISTAT-l 

PRINT ERRORS# If /ISTAT/ WAS 1 OR 2 « 

25 If (ISTAT. NE.l) DO TO 31 
If ( lERkOR.cg.U DO TO 29 



l#NPNTS 2 #l) 

4(C0Ef »XI( I1.12#I3)#X2(U#12#13 )#Xj( U# 12# 131 



8 S 0 L 0 TE ERRORS. 



wRlTfe (LO# 61J) 

00 27 13 «l#NP«TS 3. 1 






WRITE (LO# 61A) M 

S3*l! 

WRlTbCLO# 604 ) ti#(SrAT(I 1 . 12 #I 3 )*I 2 «I#^)fNTS 2 #l) 
-BITE (LO# 6 C 3 ) ( X 2 # 12 «I»NPNrS 2 »i) 

00 to 31 

PRINT RELATIVE ERRORS. 



28 WRITE (LO# 619) 

00 30 I 3 « 1 #npnTS 3 » I 
WRITE (LO# 614) [3 



WRITE (LO* 6 ^, 3 ) ( 12 # I 2 «l#wPNTS 2 #l) 

00 29 :i-i#NpJirsi.i 

WRITE (LO* 604 ) !l> (S 7 AT( Il»: 2 # 13 )#:z«l 

a# oo3) 



30 WRITE ua# o 03 ) 4 X 2 #i 2 «l#iPNfS 2 ;ii 
PRINT ERROR statistics 



#MPNTS2/X) 



31 If ( 1 ERR 0 R.E 0 . 2 I wRITc 
If (IcRROR.cO.i) wRXTe 
0<3 TO 34 



(LO# 

(LO# 



416) £»WAye#l4AX. J4AX#ArtAX. 53 PNaX 
6 17 ) ERRAVE#!’ 1 AX#DNAX#R.' 1 AX#cRRNAX 



fATAL ERROR OETECTEO. OUTPUT ERROR flESCAGE ANO PRESET COE f FI C 1 £•« TS • 

32 ITER— ITER 

33 

WRITE (LI# 101 ) rrER#EPSLON 
return* 

34 RETURN 
£N0 



XTER#EPSLaN 
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74/l7i ftr* 

ustounne .1*7 41 (H^irsi«H»irs2«N9Nr:3«i#N» i?ay0i»suiii*ii«^ud4*) 



P6ESET suan*ratz 

40-50LTNU11. ■ 






iUTHOAi 

OAfgi 

y*AUSt6 

;«i5r* 

•«^irs2 

N^«rs3 

n 



lAL *)*^ 

^4N$ ZiSM€«l 
?6C£«66« 31. 



SJ6H IM fHf* SYITCn i*T*H 4 (JJk * 



1563 



tl^LANATlQH 

.*u.;3ir5r5urjQ!;f§-u-i,v.sfiAf-ri 

1UH6R 3F 04TA POINTS /.r CONSriNf tl 
NU.ldg? OF 3*r* POINTS Ar«CU‘J-T*NT 13 
N013ca OF C06FPICIS1TS rOR fNg pOUT- 
N01UL IN * ' 

H NUiagR ■ 

101UL , 

IFQURl PQygR FOR 

JUNIATION, 

SUtfll(1M»1Nl SJeiATRlX 
1ATAIX, 

IN OIlgNSIONlMO yAflA«LS £0R ;ilONl.-k«N 

SU6N2(NfN) SUaiAT^U 2. CLASS Tn ThE ClldlATRlX 

1* CLASS* 



OF COEFFICIENTS FOR T Hi FQLT- 
IN <2 OIRECnON. 

OF CJeFFlCrSNTS FOR THg POLT* 
i.i ti OlFgCriOM. 

FOR TN€ <l y*LU6C :M the 

1. CLASS II TIC STSTCN 



TTF5 

ww 

r.MTcS 

INT! 

iNTe 

IHTe 

INT6 

RbAC 

lNTz« 

XcAL 



cR 

3SR 

•SeR 

icR 

^;cR 



75 



EAL SU8N1<NN.1NI#SU512(N»NI 
F0UR2«0 

??‘fni.‘A7»Si ra 01 

TQF goes SgCTtON SU9NATRt:gS 2. CLASS. 

ISTART*! 

istof«n 

JSTART«(<-ll*H»l 

j5tof«r*n 
60 TO 02 

RI6HT HANO goes SECTION SUSNATRlCtS Z. CLASS* 

OX ISTaRT»CR-1I*N*1 ‘ ‘ 

I§TO?"t (RPU-DPN \ 

JSTART»<N-l»*N*l 
* JST0F-R4N 

PRESET E06E SECTION SUONATRICES 2* CUSS IN SUONATRIX 1. CLASS 
USIN6 1AT42* 

92 ^ALL 1AT42 (NPNT$l#N*NTS2«NPNTS3#N«IPawP^^ IPQwR2*SU6N2) 
fii'oS :*tSTART#ISTOP#l 

ofl 03 4«4START.45TaP#l 
41«41»1 

03 suam 1.4 I«SU8l2i I1«4U 

04 IP0WR2«IP0UR2«1 

CCPY OEFINEO SU8NATRICES 2. CLASS OlACQNALLr* 

30 C6 tl«l.Nl.l 
00 06 41-2#fl#l 

istart* ni-u»H*i 
IsroPMiiN 
4START«(41*II*N«L 
4ST0P«4l*N 

30 05 I-IST48T.IST0P.1 
00 05 J-JSTAar.JjTOF. I 

05 SUSNK 1*n.4-n)*S03NK I#4I 
96 CONTINUE 

RETURN 

ENO .X 
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: 
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30 
19 



7^/i79 PTn 

suanuuuNe «4T«.2 (Npf<«fsi#*4i»*«TS2#«^‘<T?2 , m# ipqw«;# 



v-/l2/2>. t-*. : .. -1 



PRESET $U8fl 
SU3A1 FOR 4 


lATHfX 2. CLASS SJd-iZ l'\ rif. SU9»ATaiA 
. 40-«»OL 7N04I AL AFFBOaIHaT :;1N. 


L. CLASS 


AUrnORt 


<445 CEasER 




OATEi 


lfcCE49Eii 3i# I960 




variable 


EIFLA.VAT104 


TYP? 


liFNT SI 


NUiU« 3F OAfA Fofifs t} constant <I 
4U43PR OF data F3I4TS \r C1MSTA4T <2 


Nf cOeP 


4F4rS2 


£4Te",:^ 


NFNTS3 


40i3£R 3F OaTA J>0i.4rS AT C'INSrAMT <3 


tN fcOe 9 


H 


)F coeF«r:u4rj rop. thE polt* 
SOCIAL li A3 oiaEcrion. 


ISTci=R 


IFQWRl 


POV£R FOB T-<E <l VACUtS !M ’’HE 
s’jfiiAri 34. 


tHTEOER 


XF0WR2 


POvER FOR THE tz VALUES IN THE 


IHTESER 


SUBn2(N«fO 


sui.iAr xoH. 

SUdHArPll 2.. CLASS IN T-fC SU9 «aTRIX 


real 




1. CLASS. 





90 



AEAI. SUB^2(N#N) 

IP0W»3«O 

00 03 K«L«^L#1 
IF $0 T3 01 

TO^ SOOfc SiCTITW Etc^ENTS 

!•! 

flO TO 02 

K16HX HAHO EOOS SECriQM eUEiIcNrS. 

01 

F«£$cT SiCriQN ELE16NT9 IH SiJS^ATRlX 2. CLASS 'JSlwO So. 

02 sumzdfj ) •S4<MFirsi,^F*frs2.*4r'<TS3#iPawRi# xfoofi# ji 

Q3 I»0)<k3«LF0m43»L 

oo 04 

00 04 

04 SUdn2(£»L,w*U*SJd42( U4I 
«ETU«H 
ENO 
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7s/iJu <>rtQ.1» 









<* 



I 



9 



10 



13 



0 



5 



0 



3 



0 



3 



0 



5 



60 



iitAt FUNCTION :a {N^^rsi*ip'<rs2»Nf»<r;3* ipQi«3«iT I 



P6«^QH« SUrtfl4TI3SS T3 ^t6S|T fH6 CLtMiNTS A(I.J) I'i f-»2 
STSTg« ,i*rail A AMO fH6 iLE'^ENTS ?<:) !M THG JlCirtf -»»M0 
SIQe VECrC* 3 (=0« THE LI'<EA« cQUATriM irsrti, THAI <1Ui f 
it SOLvSo P3R 4 ;»3LT'«0M4L ^3-APP?.OxIflAriOH, 



AUTHOai 

OATEt 

VAilAUE 

nJ A f si** 

M^Mrs2 

M^MTSJ 

IPQtAi 

irOW«2 

IP0UB3 

IT 



S6 



HAMS 2c 3N€B 
0EC£rt9EM 31* 1990 

EX^LAMATIQM 

NuAiE«**jf 'alfi* pQ!ATS*M*$0NsfAMK5i 

NOMiga JF OAfA PQIirS \T CTfISTAMT <2 

SUMMaTIOm. 

POWER POR THE 12 WAIUCS lit THt 
SUiMATIOH* 

POWER FOR THE X3 W41UE5 IH THE 

^HmtJjl^paremeter 

•O T WALUES are hot inCUUOeO IH fHE 



•I 



SUMMATIOM. 

(-* STSTEfl MATRIX CLCffEHTSi 



T yAtUES are 
SUMMA flOM, 



IMCLUOeO IM TME 



IMT2SER 

IMlfeiER 

IMTlStW 

IMTEGER 

INTEGER 

[mteg:r 

iMTcGER 



iMTeGeR 



tlp«7* 3) 
X2Ufhsi 
t3(9#7#5I 
T(9# 7# 51 



T-> RIGHT HA10 SIOC VCCTCR ELEM.I 
result Of SUMMATION, 

XI-COOROIMATES Of THE OA^A POINTS 
r2-COJROlMATE5 Of The OAfi POINTS 
X3-C00ajlNAf £S Of THE 04TA POINTS 
T-COOROIMATES Of ThE OAfA pOImTS 



real 
real 

i<fe AC 

iiAL 

real 

IMfORTAMTi THE OaTA POINTS U# x2j Xj tNO f haWE TO dc PRO- 
yioEO Through the laacleo common slocr oataa. 



ccnnON / oaian / ii#a2,a3,t 

real X1I9* T*5WX2l9,7fOl*A3(9*7*3l#n9#7*$) 

S ■ • G 

00 07 R1*I#NPNTSI#1 

oq OT < 2 »WNPNrs 2 ,i 
00 g7 R3«I*NPNrs3#I 
If (IPOWRI.cO.U) so TO 01 
5l*XUM# A2iA3l«*IP0wRl 
G3 TO 02 

Si |}*hpawR2,£Q*OI GO TO 03 
S2«A2I<l*i<2#Ai)**IfOWR2 

03 t?.I! 

06 If ( lP0wR3.EQ.0l GO TO 05 
|j«xj(Jl#A2jK3l**IP0WR3 

05 $3«l. 

06 ST«r(KI.<2»K3) 

* If (iT.lo.CJ ST*l. 

07 S-S*USiPS2l*S3MSr 
56* S 

RETURN 

EnO 



r 
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35 



«0 



65 



':a 


• 




/6/179 


aRT«2 *»HonR 


FT/a ; 


XNrso£R fuNcrnN islaa 




OAUSS JORDAN SIXIINATIJ'I AXfH MATRIX A *N0 VtClOR 3. TM 
^AAlAiLc l|L6 IS 3cr rj 1, IP rnt ttlNlNAflON IS CaARCfilcO. 
X6L6 1£SS T^N 0 INOICAreS AN SRROR, 


AUTHOR i 
OATct 


HANS 2£3N«6 
OfeCgNSfeR n, 1680 




fARuaiE 


sxrcanatxon 


TTR6 


1SL6 


NUH36R hf SOuifiaNS 
CaNr«:L <»4RA«6r6Ri 6eTURMC0 
Calling, RauriNS • 

•1 Gauss joroan eLiNiHArrai 
ASRFORNeo COPRcCTLr. 

<0 6RR0R SnCOUNTbrco OURt'^G 
JORDAN SLXlXNATlON* 


rO TH6 [ntISeR 

4Ai 

GAUSS 


A(66f06 ) 
9<6A> 


STSTPN HArtll t 

ilSHf HANO SlOe JECTCIt 3 


R6AL 

RcAL 


INRORTANT* 


HATRIX a 4N0 JECrOR 3 HA^C TO 36 ••OVI069 THROUG-l 
fMe LASSLeO C06N3N 3Lf1Cx ..rt6, dOH A a.h 3 i ^ILL 
¥XLL UNOtRGO changes OURIMG cXECUnON# The 

ELlilUriON. 






COnnQN / Ad6 / 6«d 
R£Al A(6«#66»«6(64) 



fu:;iv 

£^S«UO0e-O3O 
3Q 06 
90 U2 

If <*.c0, I) GO ro 02 
C0NSr«-A( A* l)/4( X# II 
00 01 

A(K« J)*A(K«jMC0:4ST«a(I#J) 

If A(<#i)«*0 

01 CONflNUt 

iu )• A(AI *CQhSrM( u 
32 CQ?tTiMU6 
C3.^ST«aU« XI 
00 03 

03 Acx« j i«Ai i«i) /co*isr 

aU»II*1. 

S(X)«6(n /CONST 
90 36 
X£6a0R«0 

ir tAaSiiulil )•3T•^PS) 03 ro 96 
Xe4tOA* (SRROR^I 
96 CONrihJ^ 

(A9S(a(jn.or.«PS) so ro 39 

* ifcR«g»« tfeRROP*; 

05 If (Icaso«.»ir.N*u 03 ro 37 

06 cONrisue 
GO TO US 
I£L6 — L 

Itt<*6*L 

08 «6TURH 
£N0 



07 
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TO 

T5 

90 

95 

% 



QJ»r«2 PrtOrt* 






9C/i2/2T> 



SUtiQuriNf IEL6 (N> ITEIflTt'IAXfETSLan) 



Gaussian ELin(H4TiCf« ^irn linp 4« 19U4nQ>* Syste»i y ithe »icht 
H4N0 SIQE YECroa 3 is INTEG34rE0 TIOCTAep aXTl T»iE $rsf5>< ^4- 
TRIX 4 IH T «£ IfRIX VI. Ifga IS THC -1U19ER QP IT6R4T[Jf«S a«- 

THE «4GHITU06 :F gPSLUH. jHaULO The MUH6t« UF ITER4TI3NS cX« 



4UTHQRS* 

OATEi 



XQHR43 VQGELERf H4H$ ZC3NCR 
3ECtH9ER 3t« 1990 



V4RI49CE EXn4H4TIQN 



TTFe 



ITER 



ITRA4X 

EFSION 

4(6A»66I 

9(661 



NU19ER 3F cGUXriQNS 
NU(19ER 3F irSR4TI3NS TO 8C FETUfiNEO 
*0 ^4ussiA^( EtxniMArioN PtaFORKta 
CORRECTLY . 

<0 MO CQHVEROEiCf OF TH6 SOLUTION FOR 
the oESlaEO aESJLurroM tFSiaM. 
ITER4TIQi 310 TSRHt.MATE 4FTER 

lTRi4x Iteration step. 

NAXliUH iU«9ER OF ITERaTJTM STEFS 
OESIREO RESOLUTION OF TMC SOLUTION 



iNTj-GER 

lNTe<iER 



iNTEGcR 

real 



STST6N MATRIX 4 
RIGHT H4N0 SIOE VECTOR 



real 

real 



IIIFORTANTi matrix 4 4M0 VECTOR 9 H4VE T3 3E PROTIOEO THROUGH 
The L48ELE3 COMMON 3L0CX A03. aOTH 4 4MJ i «1LL 
NILL UnSsRGJ ChamGES OORlMC EXECUTIONi* 1.6. THc 
ECINIATIQN. 



SSmN^N f 496 / V 

RE4L V(66«o9l«G(69)«A(66*65l#9(65l 

L0GZC4L IFRImT 

0AT4 IFRINT /.F4LSE./ 

0AT4 FI / 3*161 392636/ 

FORMATS IEL6 ST4RT 

691 format (///39H matrix 4 JMQ VIKTOR S VOR ELIMINATION/} 

692 FOFMAI ((6l*l3(aXtl2MI 

til lassiF 

• I 

605 fORMAT (///2Xt I2#26H «3EFein£HTEN 4US GAUSS/I 
ORMaT ((2X,I2,H,;PBq,2M 



TERATIOMSSCHRITT/ 



606 

nl 

604 

610 



format C//V39H KLtiHSTER ZcILcNVEKTnprCHNrTTUTNKEL ->612.31 
FUAMaT (27X>I2H RESIOIUN •*6I2.JfiH HACHtlO^liH iTERAriONEHi 
<"^20XfllH ;§hl 6R 3H> -feia.i/l 

FORMAT (///25H !tST£A N4CH 6LlMZM4T!of4/) 

FORMATS ZEL6 STOP 

L0*6 

M«M*1 

K-N-l 

REARRANGE VECTOR 9 IN 4RRAT V* IF H IS LCSS THAN 66* 



rIJMT^HAmO SIOE VECTOR 3 3EF0RE THf 



It ^<3 TO 02 

00 3l I«UM*1 
01 V(I#MI«V( 14631 

ifixwnjs: uf 

« I? (inisfi :sHi it3; tiw 

93 ?? n?i iNtl'itlTE (LO* 6011 I«(V(If Jl* J«L4M«1] 
IF (IFRINTI VRITE IL0« 6021 (RfM-L^NfU 

LINEAR EQUATION STSTEM V( N» M l*G(N )«V( NfNRI ) 

§8 IX 

06 4(l4J|«V( I«VI 
00 05 I«l#N#i 
05 9m«.0 

CHICK FOR STA91LITT. 
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90 

95 

100 

105 

no 

ns 

120 

123 

130 

135 

140 

145 

130 

155 

160 

165 

170 



su86aurx;i£ 



fhtm 



FT* -» .a»«* »3 



# Jl**2 
>) 



o4*s2 

06 G(n«.c 

2? 

81 Siiiniiiirj 

00 :o 

00 10 J«JJ*M»1 
Sk« .0 

00 09 13*1«N»I 

09 S4«s<*y <1# m»y«J#i3i 

F«$K/ nt( M 

If (AssjF>.sT.Aas(Ft>> 

10 COiitlNyfc 
Ai«Aas(Fi> 



Fr«ACQS(F u 

p£.| - • 



•(Fl*iOO.)/Fl 



STA6T llfflATIOl r3 FINO *10 I/1F«0V£ CJtUTinH, 

00 27 LX*l#999#l 
IriR«i.x 
oa 13 

FIVOTISATION AHO SAUSSIAN STEF F06 CHLUnH J. 
1X«J 

t«y( J«J > 

00 12 

\f (Aa3uu#jn*iT.A3S{£M 00 ro u 
lo ro iz 

n 

12 



£«V< t*J ) 
COMf INU6 



If (IX.fcO.J) GO TO 1% 

00 U 

G(I>«yUX#l) 

y(ix*n-yu*i> 

13 y<j*i>«G(i> 

Gfi >«.0 

14 

00 13 l«yrt*Mfl 
00 la 

y4Z*AK)«y (i#AA)-y( j»<iO*9 

15 COMTlNUi 

iF'*(y(X#n«£o«.G) GO ro sz 

SOLUTIOM, 

16 

G(N)«yfN*n)/VlH*1) 

00 19 11«1#A#1 
I*l-ll 

. IF (y(l«I>«£0«.0> GO TO 32 
Gcl)*yc I#rt» 

00 17 JJ«1#A»1 
4*»4-J4 

IF (4.E0.1) GO TO 15 
17 G<II«G(I>-G(4l«y(X*4> 

19 S?l)iG(i»/y(l# X) 

lT£RATlOrt, 

00 20 

20 GIX )«d( I) »G(I) 

CQAFUTC TH£ 0C9AINX16 eRAOl* 

00 22 I«1«M»1 
y(i#rt».,a 
00 21 

ya«A)«y(i#'i»*^(4 )«*(i«4) 

<x#^»*A(i#rt»-y(t#rt) 






a .. . 

*••0 

00 23 :*1#«#1 
23 ««A^y(I#rt»»y ti#4» 
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SUSRaUftNfc lEL^ 



7Mfl7b l»»T-a 



FTh 






175 



110 



150 



15 * 



200 



205 



210 



215 



220 



225 



25 

25 

26 
27 

29 
25 

30 

31 

32 

33 



••S0«TI«) 

IF (it.cr*€PSlQ'4) <30 13 29 
IF aa.Gi , 31 ra 33 
00 S 
S(X)«<(n 



H4T81X 4 4N0 41SMr H4N0 SIDC VCCTOit ^ AFTER THE 



LX 



FRIHT STiTE*! 

FIVQtlSAf ICN. 

IF WRIT® (L3# i04l 

IF (IPrUtI mRITc (ta# 9021 
OO 23 I«UN»i 
IF UPftlNTI ^RITE (LJ# 603) I#( VI I#J) * 1 ) 

IF (iPtlNT! ^RITE kj» 5021 
30 26 
D3 26 

V(I#J)*4I1#4> 

STOF L30F TO FIMO S3LUTI05. 

COnriNUf 

REARRANGE SOlUTUN IN THE 4RR4T 

conriHUE 
OQ 25 I«i»N»I 
VII#63)«G(I) 

PRINT SOCUTIOH# IF IFRIHT IS SET .TRUE. . 

IF (IPRInT) write IL3# 503) N 

IF IIPRINM WRITE (LJ» 505) t I#Gl I )» I«l#N»:) 

IF (IprInT) wRITe UIa 507) 

IF IIPRINT) mRITE (LJ* 50()) «#LX 

ORUCFEN OER N4TRIK 4 UNO 3ES VECTORS 3 rtACH ELININ4TI3N. 

IF IIPRINT) WRITE (L3# 6l3) 

IF (IPRINT) WRITE (LJ* 902) U*J«I#N»i) 

00 U 

\i n;S{S{l :Slfl Iti: ‘55! 



return 
ERROR ENCOUNTERED. 



932) U»R*UM»I) 

ST3*. 



ITER**! 

WRiTfc (LO. 6G5) I*I»W(I«I) 

G(! )»GIi)/vn» I) 
return 

NO CONVERGENCE OF THE SOLUTnN* 

IT|R**IT£R 

RETURN 

END 
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74/173 lPr«2 PTi 

ft£AL function f4 (COEf* xi» x 2» o#NoaoRi*/nKOP.<2«aaaoiia) 



CONPUTfe FUNCTION /ALUe r-FUl#Xi#/3) OF A 40— 'UL TNJNIAL 'JilNd 
Trt* NOBNfeR :CN6N6 OBOES TO InCPCACE :?fctU amO 4CCJSACr. 

AUTHOBl NANS 

OATEi JECEA9EA U# l^flo 

VAlUacE EXPVANAriON TYPE 

isiMSU 34L?wfXt*?5i?^ia!f;f5 

XI X1'-VALJ€* srtc^E / SILL 3C COAPUTtO* 46 al 

X^ X2*/alJE» SNtSg r 4ILL dg C0**Pur?0* 4tA4. 

Xi X3-VALJE# SNE^E Y rflLL 3g COrtPUTSO. ic4L 

NOaOBl OBOEB 3F TNg polynONUL !m TNt iNTtOcS 

XI OIseCTlON. 

N0B0B2 OaOeS JF fn« pOLTNOnUl X-i The InTEJEs 

X2 3IRECTI3N, 

Nosoas oaoes of TnE POLTNONUL :« the iNTeieS 

X3 OIRECTIQY. 

FA IROINATE OF THE SO-PQLTHO.OI AC. ««AC 



Seal CJEF <&%)# 0(A4)f c(64) 

NCObf i«NOaosi *i 
NCakFj.NoaoaZst 
NCQEF 3«NOaOR3A; 

CONFUTE COEFFICIENTS 0 iASEO ON COEFFtCICMTS COEF am© X3 

§8 S! 

T-COEFUC ) 

00 C4 uI«i«Hoaoai«i 
I3 «nCOEF3-(31 

IC«1I1-;1*nCOEF2pnC0£F3»U2-IMnC0EF2p: 3 
01 f«CoiP< Ic I sx^PY 
92 OIlO)*r 

CUNFUTE COEFFICIENTS E 4A$tO ON COEFFICICMTS 0 ANO x2, 

90 OA Il*ltNC0EFl#i 
IO«I1*NCOEF2 
T *9 U 0 7 

00 03 Ui«i>Noaoa2#i 
I2*NC0«F2-121 ^ , 

10« i I :-ll •NC0EF2*I2 
r*9( lOI ♦X2*Y 
Elli) •T 

CONFUTE r SASSO ON COEFFIIIENTS E ANO XI. 

Y«6INC0EF1) 

00 w 3 lll«l«NOaORl*l 

Il*NCpEFi-IH 
05 r.EaI)*xl*T 

CA* T 

SETUBN 

ENO 



U 



tutil/ 
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